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An Indian and a white man fell to talking. The white 
man drew a circle in the sand and a larger circle to 
surround it. 

"See", he said, pointing to the small circle "that is 
what the Indian knows, and that", pointing to the 
larger circle, "is what the white man knows." 
The Indian was silent for many minutes, and then 
slowly he pointed with his arm to the east and turned 
and waved to the west. 

"And that, white man", he said, "is what neither of 
us knows". 
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Abstract 



In this thesis wc consider primarily the dynamics of quantum systems subjected to continuous 
observation. In the Schrodinger picture the evolution of a continuously monitored quantum system, 
referred to as a 'quantum trajectory', may be described by a stochastic equation for the state vector. 
We present a method of deriving explicit evolution operators for linear quantum trajectories, and 
apply this to a number of physical examples of varying mathematical complexity. 

In the Heisenberg picture evolution resulting from continuous observation may be described by 
quantum Langevin equations. We use this method to calculate the noise spectrum that results from 
a continuous observation of the position of a moving mirror, and examine the possibility of detecting 
the noise resulting from the quantum back-action of the measurement. 

In addition to the work on continuous measurement theory, we also consider the problem of 
reconstructing the state of a quantum system from a set of measurements. We present a scheme for 
determining the state of a single cavity mode from the photon statistics measured both before and 
after an interaction with one or two two-level atoms. 
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Chapter 1 

Introduction 



In this chapter we introduce quantum measurement theory and explain why it is considered to be 
philosophically unsound even though it is extremely successful. We note that realistic measurements 
often consist of observing a system indirectly by performing a measurement on a second system 
with which the first interacts. We introduce the resulting formalism of generalised measurements, 
consisting of operations and effects. If the second system remains unobserved, then the first system 
undergoes decoherence, and we examine the relationship of this decoherence to the measurement 
problem. We go on to explain how dissipation and continuous measurement result when a small 
system interacts with a very large system, or reservoir. If we choose to monitor the reservoir, and 
hence the system, then the system undergoes an evolution referred to as a quantum trajectory. We 
show that the nature of the trajectory is determined by the way we choose to monitor the reservoir. 

1.1 Overview of the Thesis 

This thesis is concerned with the measurement of quantum systems. We focus primarily on the 
dynamics of quantum systems undergoing continuous measurement processes, and the resulting 
noise on the measurement. In addition to this main theme, we consider the problem of 'state 
reconstruction', which involves determining a quantum state from a set of measurements, where 
each is performed on the state in question. In the first two chapters we introduce the aspects of 
quantum measurement theory, and the theory of open quantum systems that provide the background. 
The original work is then presented in Chapters 3, 4 and 5. Some supplementary material is also 
included in the three appendices. 

In Chapter [l] we introduce quantum measurement theory and the theory of open quantum sys- 
tems. The purpose of the first part of this introduction is to present the key peculiarities of quantum 
mechanics, that of the measurement problem and non-locality, and to sec how they are related to 
entanglement. The second part of the introduction presents the quantum trajectory formulation of 
the theory of open quantum systems in a physically motivated manner. That is, by considering a 
system undergoing a real continuous measurement process. As the theory of quantum trajectories 
developed in a number of separate ways, we also include a brief overview of the various approaches. 

In Chapter || we are concerned with the mathematics required to describe quantum trajectories. 
First we introduce stochastic calculus, which is required for the treatment of stochastic differential 
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equations, and explain why it obeys rules that are different from those for standard calculus. We then 
show that quantum trajectories may be formulated in a simple manner in terms of the generalised 
measurements introduced in Chapter ^, and show how quantum trajectories may be written in 
equivalent ways as either linear or non-linear stochastic equations. In the last part of this chapter 
we consider the Heisenberg picture formulation of open quantum systems, which results in quantum 
Langevin equations. We use these to derive the master equation which was introduced in Chapter [l]. 

In Chapter |j| we present a method for deriving evolution operators for linear quantum trajectories, 
and apply this to a number of physical examples. We are primarily concerned with trajectories 
involving the Wiener process, and which describe continuous projection measurements of physical 
observables. However, we also show that we may extend the method to provide a unified approach 
to quantum trajectories driven by both Wiener and Poisson processes. 

In performing a measurement of a physical observable, the uncertainty in that observable is 
reduced. Due to the uncertainty principle the uncertainty in the conjugate observable is increased. 
If the dynamics couples the two observables together, then the increase in uncertainty feeds back 
into the measured observable, and the resulting noise is seen in the measurement. In Chapter || 
we examine an experimental realisation of a scheme for measuring this quantum back-action noise 
in a continuous measurement of the position of a mechanical harmonic oscillator. We consider 
various sources of experimental noise and examine how the back-action noise appears among these 
contributions. 

In Chapter || we turn to the problem of the reconstruction, or measurement, of the state of a 
quantum system. This is possible if many copies of the state are available, so that independent 
measurements may be performed on each one. We present various schemes for determining the 
state of a single cavity mode, in which the photon statistics are measured both before and after the 
interaction of the mode with one or two two-level atoms. 

In Chapter ^| we conclude with a brief summary of the work presented in the thesis, and comment 
on the current state of quantum measurement theory, and directions for future work. 

We include, in addition, three appendices which contain material that we felt, while useful, was 
best separated from the main text. As we use the input-output formulation of open quantum systems, 
and as there do not seem to be any really detailed introductions to this topic in the literature, we 
have included an introduction in Appendix In Appendix [b| we include a calculation of the action 
of exponentials linear and quadratic in the single particle position and momentum operators, a 
result that we require in Chapter [| Finally, in Appendix ^|, we present some expressions which are 
required for the discussion in Chapter [si. 



1.2 Quantum Mechanics and Measurement 
1.2.1 Breaking the Rules of Probability Theory 

In the early part of this century a new theory was developed to model the behaviour of atomic 
systems Q. Referred to as 'The Quantum Theory' it broke from the inadequate classical theories 
which had preceded it by describing the world as existing in states that are sharply at variance 
with the states in which we perceive the world to exist. The relationship between quantum theory 
and measurement, or equivalently between quantum theory and our perception of the world, is still 
regarded today as a philosophical, if not a practical, problem. 
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A central peculiarity of quantum mechanics is that it breaks the rules of classical probability 
theory Q . A simple example of this is given by the electron interference experiment involving two 
slits ||. In this experiment an electron is fired at a barrier containing two adjacent slits which we 
will denote by slit 1 and slit 2. If the electron passes through the slits it falls upon a screen where 
its position is recorded. First of all we repeat the experiment many times with slit 2 blocked to 
determine the position distribution of the electron resulting from passing through slit 1. We may 
denote this distribution by P(x\l), which should be read as 'the probability distribution for the 
electron to be at position x on the screen given that it passed through slit 1'. Similarly we may 
block slit 1 and determine P{x\2). Using the rules of probability theory, we may now determine the 
distribution of the electron on the screen when both slits are uncovered. Denoting the probability 
that the electron passes through slit 1 and slit 2 by -P(l) and P{2) respectively probability theory 
then tells us that the final probability distribution for the electron at the screen is given by 

P{x) = P{l)P(x\l) + P(2)P(x\2). (1.1) 

Astonishingly, performing the experiment we obtain a completely different result for the distribution, 
one which contains the famous interference fringes. There is nothing wrong with probability theory, 
however. If each of the electrons passed through either one slit or the other, then the final distribution 



would have to be given by Eq.(l.l). This is because in this case each electron could be labelled as 
having gone through either one slit or the other, and the two sets of electrons would possess the 
distributions P(x\l) and P(x\2). Then Eq.(|l.l|) follows immediately. If we set up a measuring 
apparatus to detect unambiguously which slit each electron passes through during the experiment, 
we discover, as we must, that this time the resulting distribution is indeed given by Eq. ([□]). That 
is, the act of measuring the electrons during their passage disturbs them, and the experiment is 
altered. We are left with the conclusion that in the absence of the 'which-path' detection apparatus, 
the electrons do not pass through one slit or the other, but must pass in some sense through both. 

Quantum mechanics successfully models the world by allowing a physical system to be in a 
superposition of any of its many possible states. The coefficients in the superposition are called 
probability amplitudes, and the probability that the system is found to be in any given state is 
the square modulus of the respective probability amplitude. While in probability theory it is the 
conditional probability densities that obey the Chapman-Komolgorov equation, in quantum theory 
it is instead the conditional probability amplitudes which obey this equation. To describe the two- 
slit experiment, we assume that when the electron passes the slits, it exists in a superposition of 
being at slit 1 and slit 2, with the probability amplitudes and ip(2) respectively. Denoting the 
conditional probability amplitude density for the electron at the screen, after passing through slit 
1, as iji(x\l), and after passing through slit 2 as ip(x\2), we obtain the final probability amplitude 
density as 

ip(x) = i)(l)ip(x\l) + ?/>(2)V>(x|2). (1.2) 
The final probability density, being the square of the probability amplitude density is then 

P(x) = |V(1)| 2 Wz|l)| 2 + |<K2)| 2 |^|2)| 2 + 2Re[^(l)^(2Mz|lMz|2)] 

= P(l)P(x\l) + P(2)P(x\2) + 2Re[il>(l)i>(2)t(j(x\l)i>(x\2)}. (1.3) 



Comparing this with Eq.(l.l) we see the appearance of an extra term, and it is this which provides 
the interference fringes. 
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To describe the world, quantum mechanics requires that physical systems exist in superpositions 
of their possible states. These superpositions are completely outside of our experience, as systems 
always appear to us to be in one state or another. It is this which lies at the heart of the quantum 
mechanical measurement problem. Somehow all the things which we sense in the world around us 
must 'collapse' to one state or another, because this is how we experience them. But so far we have 
found no physical mechanism to bring this collapse about. In all quantum mechanical measurement 
theory we simply postulate that at some stage during the measurement process a physical system 
under observation collapses to one of its many possible states, with respective probabilities dictated 
by quantum mechanics. It is for this reason that quantum mechanical measurement theory is 
considered philosophical unsound even though is is extremely successful as a practical theory. So 
far it has predicted without fail every experiment undertaken, even though current experiments are 
now being performed at the level of a single atom Q). The answer to the measurement problem 
appears to hold some great secret as to the nature of the universe and our relation to it, but as yet 
no-one knows this answer. 

The measurement problem is synonymous with Schrodinger's famous thought experiment in- 
volving a cat ||. The crux of the matter is that while a microscopic particle appears to exist in 
superpositions, a macroscopic object such as a cat never does, and certainly we never do (have you 
ever experienced being in a superposition?). The problem is that there is no boundary separating 
a cat from a microscopic particle, there is only a difference of scale. So where do quantum super- 
positions end and macroscopic objects like cats begin? This is another statement of the quantum 
mechanical measurement problem. We will show below how quantum mechanical entanglement may 
be used to solve this problem in practice, although it does not provide a complete resolution. 

1.2.2 Entanglement and Decoherence 

When two or more quantum mechanical systems interact, the final state of one of the systems may 
well depend upon the final state of the others, and this situation is referred to as entanglement. To 
see how entanglement comes about we will need to first examine how to treat a system which is 
composed of subsystems. 

Let us consider two two-state systems, denoted by a and b, with basis states |<zi) and |a,2), and 
and |6 2 ) respectively. The basis states of the composite system, which consists of both a and 6, 
are therefore given by the four combinations where system a is in either of its basis states and b is 
in either of its basis states. We may denote these four basis states by \ai)\bj) = \a,i, bj) where i and 
j may be 1 or 2. If we assume some arbitrary interaction between the two subsystems, then this is 
equivalent to an arbitrary Hamiltonian in the four dimensional space of the composite system which 
we will denote by S = a<8> b. The resulting evolution will explore arbitrary states of S. In particular 
consider the state 

\B) s = ^\ai)\b 1 ) + ^\a 2 )\b 2 ). (1.4) 

In this case it is clear that if a measurement is performed on a, then the state of b will depend upon 
the result. For the purposes of rigour let us pause to obtain the result formally. A measurement on 
a alone corresponds to a projection onto a single state in the space of a, while projecting onto the 
entire space of b. The projection operator corresponding to a measurement resulting in the state 
\ai) is therefore given by 

Pa(ai)=J&®|ai)(ai|, (1.5) 
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where lb is the identity operator in the space of system b. Because the action of If, does not alter 
system b, it is standard practice to drop it, a convention which simplifies the notation. We will follow 
this practice here, and write the projector simply as V a (a-i) = \a\) (di\- Performing the measurement 
by operating on \B)s with V a {<Xx), we obtain 

V a {a 1 )\B) s = ^\a 1 )\b 1 ). (1.6) 

The probability that this particular result is obtained is then the norm of this state (being 1/2). 
We see that if a is found in state |eti) then b will be projected into state and conversely. When 
the composite system is in a pure state, we can say that the two systems are entangled when a 
measurement on one system changes our state of knowledge about the other. When the composite 
system is in a mixed state, however, the definition of entanglement is more complex, and we refer 
the interested reader to ||. 

We now wish to consider the relationship between the state of the combined system, S, and the 
states of the subsystems, a and b. If we assume that S is in the pure state 

\^)s = X! a ij\ a ibj) = \a\){otu\bx) + a\2\bi)) + |a 2 )(a 2 i|i>i) + CX22I&2)), (1.7) 

ij 

ibenpij = \oiij\ 2 is the probability that S will be found in the state \<n,bj). Therefore, a measurement 
of a (in the basis {|<ii), (02)}) will result in the state \a%) with probability pn +P12, and in the state 
\a2) with probability p 2 i +P22- If the result is to project a into the state |ai), then the state of b is 

' (aii|&i} + a 12 |6a», (1.8) 



\/pTx +P12 



which is the state of b 'multiplying' \a\) in Eq.(1.7), appropriately renormalised following the pro- 



jection. Similarly the state of b following a measurement of a resulting in the state |a 2 ) is the 
term which 'multiplies' \a,2) with appropriate renormalisation. Writing this in the density matrix 
formalism we have 



Ps = 



( Pii c\{ c{\ cf{ \ 

11 21 72 

C\2 P12 C{| 

11 12 22 

c} 2 { 4i P21 4i 

C22 c 22 c 22 P22 



(1.9) 



where c*j- = a^a* kl is the 'coherence' between the states |oj, bj) and \ctk,bi). The states of b resulting 
from a measurement of a which produces the results |ai) and I02), are, respectively, 

1 ( Pii c\\\ If P21 41 \ 

Pi = 1 11 ) , P2 = ; 21 ' (1-10) 

Pll+PU \ C 12 Pl2 J P21+P22 V C 22 P22 J 

That is, the top left hand corner of ps, and the bottom right hand corner, respectively, both 
appropriately renormalised. If the measurement is made on a by another observer who does not tell 
us the result, then we can only predict the final state with certain probabilities. We must in this 
case describe the state of b as a mixture of the two possible states, each weighted by the probability 
that the measurement of the other observer produces them. The density matrix is therefore 

P=( Pl A Cll ) + ( P2 2 1 i 4l )=Tr a [ Ps] . (1.11) 

V C 12 P12 J V C 22 P22 J 
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Note that the subsystems may have become separated by some arbitrary distance in between the 
interaction in which they became entangled, and the measurement on a. As a consequence, if the 
measurement by the (remote) observer on a changes the state of b from our point of view (that is, 
from the point of view of an observer at b), then this would constitute an instantaneous transfer 
of information from the observer at a to the observer at 6, and would be at odds with the theory 
of relativity [Q. For quantum theory to be consistent with relativity it is therefore necessary for 
the observer at b to be able to treat the state of b as p = Tr a [ps] regardless of whether or not a 
measurement has been performed on a. We will illustrate why this is the case using the example 



of the maximally entangled state given in Eq.(1.4). First we note that if the state of system S was 



\ip)s = l/\/2(|6i) + |&2))|oi), then the density matrix would be 

„_ i ( 1 1 



2 



1 1 



(1.12) 



If, however, the combined state of a and b was \B)g (given in Eq. (|l.4| )), then using p = Tr a [ps] we 
would obtain the state of b as 

,{10 
1 



P=l( n I" (1-13) 



While the density matrix in Eq.( 1.12| ) describes a superposition of \b±) and I&2), the density matrix 



in Eq.(1.13) describes a mixture in which b is either in state \bi) or in |6 2 ), each with a one-half 
probability. To see that this is the correct description as far as an observer at b is concerned, even in 
the absence of a measurement upon a so that S is still in a superposition, consider a measurement 



performed on b in the basis {|+)&, |— )&}, where |±)& = (J fox ) ± |&2))/\/2- In the first case (Eq.( 1.12 )), 
the result is always |+}& because b is simply in that state. Writing the projector as Vb(+) = |+)(+| 
we obtain the probability for this result formally as 

(Mi+W+Ms = i(a 1 |((6 1 | + <6 2 |) l+X+l (|6i) + |6 2 »|ai) = <a 1 K+||+)|a 1 > = l. 

However, for the second case, we have 

(B\vl(+)V b (+)\B) s = i((a 1 |<6 1 | + <a 2 l<6 2 |) |+)(+| (|oi>|6i) + |<z 2 >|& 2 » 



2 

|(ax|(6i| l+X+l |a 1 )|6 1 ) + i(a 2 |(6 2 | |+)<+| |a 2 )|6 2 ) 
I + I = I 

4 ' 4 2' 



(1.14) 



We see that in the second line the orthogonality of \a\) and |a 2 ) has cut the expression into two 
parts, each one corresponding to a different subspace of a. Consequently the projector onto the 
state \+)b in b acts separately in each subspace, there being no interference between subspaces. The 
probability of obtaining the state \+)b is one-half in each subspace of a, and as the probability that 
the particle is in each of these subspaces is one-half, we obtain the total probability of one-half. This 



is, of course, precisely the result obtained with the density matrix given in Eq.( 1.13 ), which assumes 
that there is no interference between the states |6i) and |& 2 ). 

We have seen now that entangling a quantum system, b, with another quantum system to which 
the observer has no access, destroys the interference, or coherence between the states of b as far as 
the observer is concerned. This may be used for a partial solution of the quantum measurement 
problem, introduced in the previous section, in the sense that it explains why macroscopic objects 
do not appear to observers to exist in superpositions; they are entangled with many other systems 
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constituting the environment in which they exist. This entanglement mechanism for a loss of coher- 
ence is referred to as environmentally induced decoherence (EID). However, EID does not actually 
solve the measurement problem. It does not explain how it is that the observer experiences one 
outcome or the other when performing a measurement on b. For a single outcome to result, the 
state of the composite system, S 1 , must collapse, and as a result the measurement problem remains 
unsolved. 



1.2.3 Entanglement and Generalised Measurements 

In the previous section we examined the result of entangling two quantum systems from the point 
of view of an observer which has access to only one of the systems. An alternate question that 
we may ask is, what is the change in our state of knowledge about one of the systems given that 
we know the result of a measurement performed on the other. This is the theory of generalised 
measurements [^|, |^, |l0| . This formalism is general enough to describe every possible operation to 
which a system may be subjected, and is essential because real physical measurements are very 
rarely described accurately by direct projection measurements. 

We will refer to the system which undergoes the direct projection measurement as the probe, and 
the other simply as the system. We assume that the system is initially in some unknown state, pi, 
the probe is prepared in some known pure state, p p — \ip){i/>\p, and their interaction up until the 
measurement is described by the evolution operator U. The state of the combined system just prior 
to the measurement on the probe is then 

Ps = U p l <g> p p u\ (1.15) 

The measurement consists of a projection onto a state in a chosen basis of the probe. Denoting this 
basis by {|ra) p }, where the index n may be discrete or continuous, the unnormalised state of the 
system following the measurement is 

p n = (n\U pi <g> p p t/ |n) = {n\U\ip) p pi (^IpC 7 * \n) = n n pi Q n , (1.16) 

where we have defined the measurement operators fi„ = (n\U\ip} p , and the tilde signifies that the 
density matrix is not normalised. The probability for obtaining this result is 

Pr[n] = Tr[p n ] = Tr[tt„ Pi rf n ] = Tr[f>^„ p t ]. (1.17) 
It follows immediately from this last equation that 

XX fi ™ = 7 ' (1-18) 

n 

as the sum of the probabilities for each of the possible results must be unity. In fact, it can be shown 
that this is the only restriction which the operators f2„ must satisfy. The set of operators is referred 
to as a Positive Operator- Valued Measure, or POVM, because it associates a positive operator, 
H, 1 n fl n , with every measurement outcome, rather than just a probability [ pT| . It is possible to show 
that given any set of operators tt n which satisfy Eq.( |l.lq ), it is possible to find a corresponding 
indirect measurement which is described by that set fi [10]. The action of the operators in the 



last line of Eq.(1.16) transforms density matrices (bounded positive operators) to density matrices, 
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and may therefore be called an operation. The positive operator W n = & n Q n , which appears in the 
expression for the probability of detecting the result n, is referred to as the effect for that result. This 
general formulation of quantum measurements was known originally as the 'theory of operations and 
effects' |l(J. However, it is common to refer to a generalised measurement simply as a 'POVM'. Note 
that if all the fl n are chosen to be orthogonal projection operators, then this formalism reduces to 
that describing simple direct projection measurements. 

The theory we have introduced above is not completely general however, because it assumes that 
we have complete information regarding the initial state of the probe, and complete information 
regarding the result of the measurement. A measurement for which this is the case is referred to by 



Wiseman J12J as an efficient measurement. Because in this case p n (Eq.( |1.16|) ) can be written as the 
outer product of state vectors, an efficient measurement on a system initially in a pure state results 
also in a pure state. To describe the more general case, in which we have only partial information 
about the initial probe state, and/or partial information about the result of the measurement, we 
need specify measurement operators subtended by two indices, f2 n ,m- The first index specifies our 
possible states of knowledge after the measurement, and the second the various possible actual 
measurement results consistent with each particular state of knowledge. The state of the system 
after the measurement is then 

II I 

and the set of effects for this kind of measurement is then 

W n = J2^n < Am. (1.20) 
n i 

To see that this is also sufficient to describe the situation in which the initial state of the probe is 
mixed, we need only substitute a mixed state into Eq.( 1.16| ). If we were to make the measurement 



but ignore the result completely, then the final state of the system, pf, would be a mixture of all 
the possible resulting states, weighted by their respective probabilities: 



Pf 



n \.r n \ „ „ 



This is referred to as the non-selective evolution generated by the measurement. We note finally 
that the non-selective evolution, that is, pi and p/, is not enough to determine the measurement 
operators. In particular, any two sets of measurement operators which are related by a unitary 
transformation will generate the same non-selective evolution. That is, given 

n I 

we have 

m,n,/ n,l n 



E 



1.2.4 Entanglement and Non-Locality 

The existence of superpositions leads to another incredible feature of quantum mechanics: that it is 
non-local. By this we mean that the results it predicts cannot be explained by a theory that is real, 
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causal, and and also local. The term real in this context is used to mean that probability distributions 
are non-negative, and by causal it is meant that later events cannot influence earlier ones. By local 
we mean that events at spatially separated points cannot affect each other instantaneously, in that 
any effect from one to the other will be delayed by at least the time that light takes to travel 
the intervening distance. However, the non-locality of quantum mechanics coexists peacefully with 
special relativity because it does not imply superluminal signalling. That is, it is not possible to use 
the non-locality inherent within quantum mechanics to communicate faster than the speed of light. 

The reason we know that the predictions of quantum mechanics cannot be explained by a real, 
causal and local theory is due to the work of Bell, who showed that the correlations between quantities 
measured at remote points must obey certain inequalities if they are to be explained by such a 
theory |13 |. Because it is the correlations between remote measurements which are non-local in 
quantum mechanics, the separated observers must communicate their results to each other in order 
to see any non-local effect, and it is for this reason that this kind of non-locality cannot be used for 
superluminal communication. 

Bell's argument runs as follows. Assume that two particles (or equivalently two physical sys- 
tems) interact. The state of the combined systems is represented by some quantity A, and at the 
end of the interaction we allow this state to be uncertain so that it is given by some probability 
distribution -P(A). The particles are then separated, and a measurement is made on each. We allow 
the measurement at each location to depend on local variables 9 and <f>. The assumption of locality 
implies that the result of the measurement on the first particle may depend upon the local experi- 
mental arrangement at the location of the first particle, given by 9, but not upon the arrangement 
at the location of the second particle, which is given by (f>. Bell was able to show that under these 
conditions certain inequalities govern the correlation functions of the measurement results. 

We examine now briefly a quantum optical example which violates Bell's inequalities and requires 
just a single photon [[l4| The configuration consists of a single photon incident upon a beam 
splitter, in which the other input port is the vacuum. The two output modes, denoted by a and 
b respectively, are subjected to homodyne detection. The combined state of the output modes 
following the arrival of the photon is given by 

M = ^(|0,l>+i|l,0», (1.24) 

where the ket \n, m), denotes n photons in mode a and m photons in mode b. The photon can have 
either been reflected or transmitted, and the actual state is a superposition of these two possibilities. 



From subsection 1.2.2 we see that this is an entangled state of the two modes. Let us now set the 
homodyne detector at A to measure the quadrature xg = \[ae %e + a) e~ t6 ) and the detector at B 
to measure = ^(be 1 ^ + We^ 1 ^). In this case 9 and (j> are the local variables upon which the 
measurement result depends. Calculating the joint probability distribution for the measurement 
results x and y, given a choice of 9 and <p, we have 

P(x, 9; y,4>) = \ {tp\x, 9, y, 0) | 2 = -{x 2 + y 2 + 2xy sin(0 - ^e^+A (1.25) 

where \x, 9, y, </>) is a joint eigenstate of stg and y^. The third term in the expression for the proba- 
bility distribution shows that x and y are correlated, and this is a direct result of the interference 
between the two states in the superposition which makes up \ip). We have, therefore, a correlation 
between measurements which are performed at spatially separated points. However, the crucial 
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point is whether these correlations can be described by a local theory, and for this we turn to Bell's 
inequalities. 

Each homodyne detector consists of a beam splitter with which the field to be measured is mixed 
with a local coherent field. The intensity of each output port is measured and the difference between 
them is proportional to the measured quadrature. Which quadrature is measured is determined 
by the phase of the local coherent field. By considering the four intensities which are measured in 
the two homodyne detectors as the quantities on which locality is imposed in the manner described 
above, it is possible to show that the correlation function jl6| 

m<l,) = wfi' (1 - 26) 

where a is the amplitude of the local coherent field, must satisfy the inequality 

- 2 < E(6, 0) - E(6, <j>') + E(6', 0) + E{6', 0') < 2. (1.27) 
Evaluating the correlation function {xei)4>) usm g P(x,0;y,(p), we obtain 

E{9,(p)= , 1.28) 

\a\* + 1 

and it is readily verified that for small a (in particular for \a\ 2 < y/2 — 1) the inequality is violated. 

The investigation of non-locality does not end with Bell's inequalities however. By considering 
both direct photo-detection along with homodyne detection, Hardy has shown that non-locality may 
be seen with this system in a much more direct way without the need to use inequalities ]lE{] , In 
addition, Greenberger, Horne and Zollinger have shown that when entangled states of three systems 
are used, it is possible to demonstrate non-locality in a single measurement, rather than with the 
many runs which are required to determine the correlation functions used in Bell's inequalities Jl7[ . 

Recently it has been realised that remote entangled systems may be used to transfer quantum 
states from one place to another, referred to as 'quantum teleportation' JlS| ), and for secure commu- 
nication, or 'quantum cryptography' It has also been discovered that many weakly entangled 
remote systems may be used to obtain much fewer highly entangled systems by local operations 
and classical communication, a process which is referred to as 'state purification', or 'entanglement 
distillation' p0[ . These recent discoveries have shown us that our current understanding of entangle- 
ment is far from complete |(| , and has lead to the exciting and rapidly expanding field of quantum 
information theory. 



1.3 Open Quantum Systems: Master Equations 

A quantum system is referred to as open when it is coupled to a large system possessing many degrees 
of freedom. In the limit in which the large system contains a continuum of natural frequencies, 
we find that the evolution of the reduced density matrix of the system is irreversible. We can 
understand this by considering a classical pendulum, S, coupled to a set of classical pendulums with 
various frequencies of oscillation. If pendulum S is set in motion, and is coupled to just one other 
pendulum, then the energy initially possessed by S oscillates back and forward between S and the 
other pendulum. If S is coupled to a number of pendulums with different frequencies, then it takes 
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longer for all the energy initially in S to return. At first the energy flows to the other pendulums, and 
will return back to S from each at different times due to the different oscillation frequencies. Thus 
for all the initial energy to return to S we must wait until the various oscillations have come back 
in synchronisation. Clearly as the number of oscillators increases and their respective frequencies 
become closer together the time that it takes for the energy to return to S increases. In the limit of a 
continuum of frequencies the initial energy never returns to S and we have an irreversible dissipative 
evolution. In addition, we saw in subsection 1.2.2| that the entanglement of the small system with 



the large system, which is the result of the interaction, will cause the reduced density matrix of 
S to decohere. Hence open systems undergo both irreversible dissipation and decoherence. In the 
theory of open systems the large system is referred to as a (heat) bath or reservoir, and models the 
environment of the small system. Naturally all systems interact with their environment, and the 
study of open systems is important in all cases in which this interaction cannot be ignored. 

When the correlation time (memory) of the reservoir is much shorter than the time scale of 
the dynamics of 5*, a Markovian equation may be derived for the reduced density matrix of S by 
tracing over the reservoir, and this is referred to as a master equation. The exact form of the master 
equation will depend upon the nature of the system-reservoir interaction. In 1976 Lindblad showed 
that in order to preserve the trace and positivity of the density matrix, a master equation must be 
able to be written in the form 

p=-^[H,p}+Y / V[c n ]p, (1.29) 

n 

where T>[c n ] is the superoperator defined as 

T>[cn] = c n pc n - \c n c n p - p\c n c ni (1.30) 

and the c n are any bounded operators plj]. This is now referred to as the Lindblad form. Each 
superopertor V[c n ] represents a source of decoherence, and may also represent a source of loss. The 
master equation describing a single mode of an optical cavity in which one of the end mirrors is 
partially transmitting, and therefore possesses one source of loss, is 

p=-jr[H C) p]+lD[a], (1.31) 

where H c — hu) c a'a is the free Hamiltonian of the cavity mode, a is the annihilation operator 
describing the cavity mode, to c is the cavity mode frequency, and 7 is the rate at which photons 
leave the cavity via the partially transmitting mirror. To see that this master equation does describe 
the irreversible loss of photons from the cavity, we may use it to obtain the equation of motion for 
the average number of photons in the cavity, which is 

— {a a) = Tz\pa a] = —7(0*0), (1-32) 
at 

and this clearly identifies 7 as the photon decay rate. We have now said all we wish to say about 
master equations at this point. We will derive the master equation for a lossy optical cavity in 



section 2.3 in the next chapter, via the use of quantum Langevin equations. Standard derivations 



using the Schrodinger picture may be found in references 
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1.4 Open Quantum Systems: Quantum Trajectories 
1.4.1 Direct Detection 

Master equations, which we introduced in the previous section, are adequate to describe the evolution 
of an open system if the observer has no access to the environment. In this section we describe 
how we may determine the evolution of the state of the system when we have information about 
the environment. For the case of a damped cavity this means performing measurements on the 
light that leaks out. In this description the state of the system will therefore be conditioned on 
the results of our observation of this output light. We describe in what follows the treatment of 
Carmichael ]2^, p5|| , and continue to use the damped cavity as an example. 

Consider an open system consisting of a cavity which is damped through just one end-mirror. 
The light travelling away from the system may be written as the sum of a contribution emitted 
by the cavity through the end-mirror, and the reflection of the incident field. In particular, if we 
examine the field just outside the lossy mirror, we may write 

bout(t) =b sys (t)-b in (t), (1.33) 

where 6 Q ut (i) is the negative frequency part of the electric field operator for the field travelling away 
from the system, b in (t) is the equivalent operator for the field travelling toward the system (the input 
field), and & sys is the contribution radiated by the system. In particular we have b sys (t) = ^/ja(t). 
For a detailed derivation of this relation, the reader is referred to Appendix A. Glauber's theory of 
photo-detection gives the average photo-counting rate as proportional to the expectation value of 
the product of the positive and negative frequency components of the electric field operator, and 
the input and output field operators in the above equation have been scaled such that the detection 
rate for a perfectly efficient detector is given by (b on t(t)b^ out (t)} . 

We want to consider the detection of the light radiated by the system, so we will choose the 
external field to be in the vacuum state, and as a consequence the photon detection rate is given 
by the system contribution alone. Let us consider now the exclusive probability densities for the 
detection of photons. By this we mean the probability density that in a time interval [t, t + T] 
we detect a photon at the times t\,t2, . . .t m and no times in between, and let us write these as 
Pm(ti, tii ■ ■ ■ tm', [t, t + T}). Carmichael was able to show, using photo-detection theory J2(| ^7|, that 
the exclusive photo-detection probability densities may be written as 

p m ( tl , t 2 , ...t m ;[t,t + T}) = Trle^-Wr-^J . . . Je^~^-^Je^^ p(t)}, (1.34) 

where C is the superoperator which gives the time evolution of the reduced density matrix of the small 



system, (for the cavity this is the master equation given by Eq.(1.31)), and J is the superoperator 



defined by J p = b sys (0) pb sys (0) , which, for the damped cavity is given by japa'. Here, and in what 
follows, operators without a time argument denote Schrodinger picture operators. This expression 
leads directly to a treatment of the evolution of the system in the following way. Let us assume that 
we have the density operator for the state of the system at time t, given that we have detected photons 
at times ti,t-2, ■ ■ - tm in the interval [0,t), and let us call this density operator p c . The probability 
density for the detection of a photon at time t should, by the formula presented above, be given by 
Pm(t\ti, . . . ,t m ) — TrjjTpc]. We can now find p c by calculating this conditional probability density 
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by using the ratio of two probability densities: 

Pm{t\h,...,t m ) = — — =Tt[J Pc \. (1.35) 

Pm(tl,t2, ■■ - t m , [0,t)) 

From this we obtain the conditioned density matrix at time t, which is p c = (/5 c )/(Tr[/5 c ]), where the 
unnormalised density matrix is given by 

~ c = e (c-j)(t-t m)j Je (c-j)(t 2 -t 1 ) Je (c-j)t lp{0) _ (1 36) 

We see that we may interpret this equation for the conditioned density matrix in the following way. 
In between photo-detections the density matrix evolves according to the evolution super-operator 
exp[(£ — J)t], where r is the time between detections, and upon a detection the density operator 
'jumps' by the action of the super-operator J. In this case we have perfect detection, and the 
two types of evolution preserve the purity of an initially pure state. As a consequence, choosing 
an initially pure state allows us to write a trajectory evolution for the system state vector. It is 
readily shown that the evolution by exp[(£ — J)t] is equivalent to evolution by the non-Hcrmitian 
Hamiltonian , H e g = H — ihbl ys b sys /2 and the quantum 'jump' is given by the action of b sys . At 
time t the probability per unit time for a jump to occur is given by Tr[p c (t)] = (b sys (t)bl ys (t)} . This 
gives us a method with which to simulate the evolution of the conditioned density operator. We may 
chop time into small intervals, and during each calculate the probability that a detection occurs. 
If it does we apply the jump operator and renormalise the state vector. If it docs not we evolve 
using the non-Hcrmitian Hamiltonian and also renormalise appropriately. In fact, we may use this 
simulation process to cast the description of the evolution of the state vector into quite a different 
form by deriving a stochastic Schrodinger equation as follows. We define a stochastic increment, dN 
which is usually zero, but takes the value unity at a random set of discrete points. This is described 
as a point process. The probability per unit time for dN to be unity is given by (b sys b Bys ). If dN is 
zero then in an infinitesimal time interval dt the change in the state vector is 

m + at)) = w)) + ^ ^<'^<'» - \n - \mh (i.37) 

where H is the Hamiltonian of the system, and to derive this expression we have simply written 
down the re-normalised state after evolution through dt using the non-Hcrmitian Hamiltonian and 
expanded to first order in dt. In the case that dN = 1 the new state vector is 

\^{t + dt)) - . (i.3 8 ) 

(blys(t)bsys(t)) 

The stochastic Schrodinger equation may now be written as the combination of these two possibilities: 



d\m) 



' (bl ys (t)b sys (t)) bj ys b S ys \ + m ^ ( b sys \j;(t)} 



^(blys(t)b S y S (t)) 



mm- 

(1.39) 

Note that this equation is nonlinear because (b sys (t)b sys (t)} is a non-linear function of \ip(t)). This 
stochastic equation gives the evolution of the state vector conditioned on the measurement record, 
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and the equivalent master equation gives the unconditional evolution of the density matrix, being 
the average over all possible measurement records. Consequently, the master equation is readily 
recovered from the stochastic Schrodinger equation by averaging over the possible values of dN for 
each infinitesimal time interval dt. 

Before we complete this section, we wish to note that the equation for the unnormalised con- 
ditional density matrix, Eq.( 1.36| ), provides an alternative formulation of quantum trajectories in 



terms of a linear stochastic Schrodinger equation. We see from this equation that the unnormalised 
conditional density matrix is given by the action of linear operators on the initial density matrix, 
and we can therefore write a linear stochastic Schrodinger equation for the equivalent unnormalised 
state vector using the same procedure as for the derivation of the non-linear version above, but this 
time not bothering to normalise at each step. The resulting stochastic Schrodinger equation is 



d\i>{t)) = \dt(-iH/ti - b sys b sys /2) + dN(t) (b sys - l)j (1.40) 

where we use a tilde to denote the unnormalised state vector. While this equation now looks linear, 
a non-linearity remains in that the probability per unit time for the stochastic increment, dN(t), 
still depends in a non-linear manner upon the state vector. However, it turns out that we may also 
eliminate this non-linearity by simply choosing a fixed probability per unit time for the stochastic 
increment. This means that the stochastic equation will no longer choose the trajectories with the 
correct probabilities. Nevertheless, this equation provides a complete description of the quantum 
trajectories because the correct probabilities can be obtained from the norm of the final unnormalised 



state vector, as indicated by Eq.(1.34). The linear formulation is useful as it is suggests methods of 
solution which are obscured in the non-linear formulation. 

In this section we have shown how a consideration of direct photo-detection leads to a formulation 
of the evolution of a system in terms of a quantum trajectory, which is said to unravel the master 
equation. In this case we could write a stochastic Schrodinger equation for the state of the system 
in terms of a point process which describes the photo-detections. In the following section we show 
how a different kind of measurement process leads to a different unravelling of the master equation 
in terms of a Gaussian stochastic process. 



1.4.2 Homodyne Detection 

Homodyne detection consists of mixing the field to be measured with a local coherent field at a 
beam splitter, and then performing photo-detection on the resulting field. The photo-current then 
contains information about the quadrature operators rather than just the field intensity, as is the 
case with direct photo-detection. The beam splitter is arranged so that the field to be measured is 
transmitted into the output which is detected, and the local oscillator is reflected into this output. 
The detected field is then given by 

e(t) = Sy(Jrja(t) - iy/T^a), (1.41) 

where 77 is the transmission coefficient for the beam splitter, and we have been able to replace the 
cavity mode operator for the local oscillator by the coherent amplitude a because a coherent state is 
a eigenstate of this operator. We now assume the limit in which both the transmittivity of the beam 
splitter and the coherent amplitude of the local oscillator are large so that we may approximate 
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rj 1 and set (3 = yl — 77a. The jump operator for the quantum state is now C = y/j(a(t) + i(3), 
and the average rate of photo-detection is 



where we have chosen the phase of the local oscillator so that (3 is real. With this choice it is the phase 
quadrature Y = —i(a — a ) that is measured. Selecting an arbitrary phase for the local oscillator 
allows the measurement of an arbitrary quadrature. The effective Hamiltonian which generates the 
smooth evolution between jumps is 



The photo-detection rate due to the local oscillator [0) is taken to be much larger than the 
rate for the measured field, so that the term (a)(t)a(t)) may be ignored in the photo-current signal. 
In this case, apart from the constant rate of detection due to the local oscillator (which may be 
subtracted) , the photo-current signal is proportional to the phase quadrature of the measured field. 
Naturally the measurement process could be simulated as in the previous section, with the jump 
operator applied at each photo-count. However, as we have chosen the local oscillator to have a much 
larger intensity than the measured field, most of the photo-counts are due to the local oscillator, 
and cause only a small change in the state of the system. On the time scale upon which the system 
changes, the photo-counts look therefore more like a continuously varying signal. It would also be 
more advantageous to represent the photo-current as a continuous, albeit stochastic, signal, rather 
than as a series of photo-detections. The method we now use to achieve this was due originally to 
Carmichael, and was re- worked more rigorously by Wiseman and Milburn [ ^8| . 

Consider the evolution of the state vector over a time increment At, in which the number of 
photo-detections, m, is much greater than unity. In particular we have (m) ~ r yP 2 At 3> 1. After 
this time step the unnormalised state of the system is given by the action of m jump operators, and 
m + 2 smooth evolution operators for a time totalling At. As we wish to take the limit in which 
At — > 0, (but in which we scale (3 with At so that (m) tends to infinity), we may swap the order of 
the smooth evolution and jump operators because they all commute to first order in St. This allows 
us to write 



Careful inspection of the first part of this expression (the part in the square brackets) shows that 
it does not go to unity as (m) goes to infinity. However, as it is just a number which affects the 
normalisation and overall phase we may discard it. The rest of the expression operating on the initial 
state tends to unity as At — > so long as we scale 7 At and (3 to ensure that while (m) w r yf3 2 At — > 00, 
we have j/3At — ► 0. Now we need to consider the statistics of the number of counts, m. This is 
naturally governed by a Poisson distribution which may be well approximated by a Gaussian when 
the mean is large. In particular we may write 



(C\t)C{t)) = 1 (3 2 + 1 (3(Y{t)) + 1 (a\t)a(t)) 



(1.42) 



H eS = H - ih-y(<Ja - 2i(3a + (3 2 )/2. 



(1.43) 




(1.44) 



m = 7/5' 



' 2 At + 1 (3{Y{t))At + y/^(3AW, 



(1.45) 



where AW is Gaussian distributed with zero mean and variance (AW 2 ) = At. Expanding the 
expression for \ip(At)) to first order in ^f3At, and replacing the increments with infinitesimals, we 
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obtain a stochastic Schrodinger equation for the unnormalised state vector. This is 



\4'{t + dt)) = 



I - ( j-H + lata + ia-yiY) 



where the integrated photo-current is 



N(t) 



1(3 I (/?■ 
lo 



(Y(t'))) dt' 



dt — ia^/jdW 



V7/3 / dW. 
lo 



\m), 



(1.46) 



(1.47) 



In these equations dW is the Wiener increment, and requires the use of stochastic (Ito) calculus 
which will be introduced in the next chapter. Note that this equation is for the unnormalised 



state vector. However, an equation for the normalised state vector, the equivalent of Eq.(1.39) for 
homodyne detection, is easily derived by including in the equation a renormalisation at each time 
increment. Calculating the norm of \ip{At)) we have 



(i/j(At)\ip(At)) = 1 + j(Y) 2 At + ^/(Y)AW. 



(1.48) 



Dividing \ip(At)) by the square root of the norm, and expanding to first order in At we obtain the 
stochastic Schrodinger equation for the normalised state vector as 



\ip(t + dt)) = 



1 - 



-H 



7 



(a^a + (Y(t)) 2 /4) + ia 1 {Y{t)) I dt - (ia + (Y(t))/2)y/jdW 



As with the stochastic Schrodinger equation based on the Poisson process (Eq.(1.39)) the density 
matrix at time t is given by the average over the states generated by all the possible trajectories. 
We may readily recover the master equation by averaging \ip(t + dt))(ip(t + dt)\ over the possible 
values of the stochastic increment dW . 



1.4.3 Other Approaches 

In presenting quantum trajectories we have used the treatment of Carmichael to make it clear that 
they may be obtained from physical considerations such as the measurement of the field radiated by 
the system. However, the treatment of systems in terms of quantum trajectories did not develop in 
this way. We complete our discussion of this topic by giving a brief overview of the various approaches 
which have been used to arrive at quantum trajectories and stochastic Schrodinger equations. 

The beginning of the quantum jump formulation can be traced back to the Photo-detection 
theory of Srinivas and Davies [^9| in 1981, to the calculation of waiting-time distributions for the 
emission of photons in resonance fluorescence by Cohen- Tannoudji and Dalibard pQ] in 1986, and the 



more general calculation of exclusive probability densities for this case by Zoller et al. 31 in 1987. 
The approach of Srinivas and Davies was to develop a description of photo-detection from the theory 
of generalised measurements. The quantum jump description was essentially completely captured 
by their theory, although it was not clear how it was connected to a physical measurement process, 
and they did not employ it as a procedure to simulate the master equation. The approach of Zoller 
et al. was quite different. In 1986, Cohen- Tannoudji and Dalibard [ pO| had suggested using waiting- 
times rather than second order correlation functions to examine resonance fluorescence. Zoller et 
al. built on this by using a method for treating resonance fluorescence developed by Mollow J3^] to 
derive expressions for the exclusive emission probability densities. They found that the expressions 
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they derived had the mathematical form of the Srinivas and Davis theory. In 1992 Hegerfeldt and 
Wilser |}3| and Dalibard et al. [Q showed that it was possible to simulate a master equation 
using a stochastic wavefunction evolution in which jumps occurred at the emission times (ie. a 
wavefunction equivalent of a Monte-Carlo method), and the quantum jump formulation proper was 
born. Very quickly this became widely used because it meant that a pure state could be employed 
in the simulation rather than a density matrix, reducing the required memory. 

The formulation in terms of wavefunction stochastic equations driven by the Wiener process was 
first formulated by Gisin |3^] in 1984, using formal measurement theory arguments, and a little 
later by Diosi |}6), Belavkin J37j and Barchielli |38|]. The latter applied the theory to homodyne 
and heterodyne detection, and so in this sense his treatment is closest to Carmichael's formulation. 
Rather than just interpreting quantum trajectories as the evolution of a system during a measure- 
ment process, Percival |39| has also suggested that stochastic equations might be used to solve the 
measurement problem. If they were viewed as a fundamental theory of quantum dynamics, they 
would provide an intrinsic source of collapse. 

In 1992 Carmichael showed that the detection theory of Srinivas and Davies, and hence the 
quantum jump formulation, could be derived from the standard theory of photo-detection |2^, p7|, 
and in a sense this completed the theory. Carmichael's treatment has been filled out and made more 
rigorous by Wiseman and Milburn |28[ , and in the last few years the theory has been applied by 
many authors p0[ . 

1.4.4 A Continuous Measurement of Position 

In Chapter |3| we will be concerned with deriving evolution operators for stochastic Schrodinger equa- 
tions, and in particular with those corresponding to continuous measurements of physical observ- 
ables. In this section we would like to give an example of a real physical scheme for the measurement 
of position which is described by such an equation. This will also serve to show that the experimental 
arrangement considered in Chapter || corresponds to a continuous measurement of position. 

The measurement scheme consists of bouncing a light ray off a mirror and measuring the phase of 
the reflected light to determine the position of the mirror. To treat the problem using the language 
of open quantum systems discussed above, we will choose the mirror in question to be one side of 
an optical cavity. The other mirror of the cavity is fixed, and is not perfectly reflecting so that 
the cavity may be driven by the light beam through this fixed mirror. In the bad cavity limit we 
expect that the phase of the output light will provide a continuous measurement of the position of 
the movable mirror. Allowing the position co-ordinate of the movable mirror to be confined by some 
arbitrary potential given by i? m , the Hamiltonian of the cavity-mirror system is pl[ 

H = H m + fiLUoa^a - hga^aQ - ME{ae iu}ot ~ aV^ *), (1.49) 

where loq is the frequency of the cavity mode, g is the coupling constant between the cavity mode 
and the mirror co-ordinate (the expression for which need not concern us here), Q is the position 
operator for the mirror and the term proportional to E results from the coherent driving (E 2 /j gives 
the rate of photons impinging on the cavity) . Dropping the free Hamiltonian of the cavity mode by 
moving into the interaction picture, the master equation for the cavity-mirror system becomes 



W = -~[H m ,W\ - [E{a~a)),W] + i[g<JaQ ) W] +-fD(a), 



(1.50) 
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where 7 is the decay rate of the cavity. We are interested in the dynamics of the mirror when the 
light is subjected to phase sensitive detection. We are also interested in the limit of large cavity 
damping, and so it is convenient that in this regime we may adiabatically eliminate the cavity mode 
and derive a master equation for the mirror alone. To proceed we note that the steady state of 
the driven cavity in the absence of the interaction with the mirror, is the coherent state \ao) where 
ao = 2E/ r ). Allowing the strength of the interaction to be much less than the damping rate, we 
transform the state of the system using 

W' = D{-a Q )WD ] {-a Q ), (1.51) 

where D is the displacement operator pi[ , so that the steady state of the system is very close to 
the vacuum. Transforming the master equation using this 'displacement' picture, we obtain 

W' = -~[H m , W] + ig[(ala + \a a \ 2 )Q, W] + tga [(a - o+JQ, W] + >yV(a)W'. (1.52) 
n 

Writing the density operator explicitly in the number basis we may use an approximate solution of 



the form |28 



W = p |0)(0| + ( Pl |l)<0| + H.c.) + P2 |l)(l| + (p' 2 |2)(0| + H.c). (1.53) 

Writing out the equations for po, pi, P2 and p' 2 we find that the off-diagonal elements p\ and p' 2 may 
be adiabatically eliminated (slaved to the on-diagonal elements). The resulting master equation for 
the density matrix of the mirror, given by p = pa + p 2 , becomes 

P = -~[H m - hg\a \ 2 Q, p] - k[Q, [Q,p]], (1.54) 

where k = 2g 2 \ao\ 2 /j. The second is just a linear potential that the mirror experiences due to the 
pressure from the light. The last term, however, is the classic form for a continuous measurement, 
the effect of which is to diagonalise the density matrix in the position basis. To complete our analysis 
we need to see how the measurement proceeds during a quantum trajectory. For this we require 
the jump operator corresponding to a photo-detection. As it is the phase of the output light that 
provides a readout of the position of the mirror we will use homodyne detection. The jump super- 
operator when written in terms of the full density matrix is therefore JW = 7(0 + i/3)W(a< — i/3). 
In keeping with the adiabatic approximation we wish to refer only to the density matrix of the 
mirror, and we therefore want to find the equivalent jump super-operator which acts on the mirror 
density matrix alone. To do this we slave p 2 to po, and this allows us to write W in terms of p by 
substituting into Eq.( 1.53| ). The requisite jump super-operator is then given by 



JmP = Tr c [JW] = Tr c [ 7 (a + a Q + i(3)W'(a? + a - i/3)] = CpC\ (1.55) 

where 



C= (V7/3-*V7«o + V2fcQ). (1.56) 

We can see that if we choose direct photo-detection (j3 — 0) the average count rate is (C'C) — 
7|ao| 2 (l + (4g 2 /7 2 )(Q 2 )), which contains no information about the average value of the mirror posi- 
tion (only the second moment). Using homodyne detection, however, and defining (3 appropriately, 
we have 

C=0 + V2kQ), (1.57) 
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so that the average rate of photo-detection is 

(C^C) =P 2 + 2V2fc/3(Q) + 2k(Q 2 ), 



(1.58) 



and provides a continuous measurement of the position of the mirror. Taking the homodyne detection 



limit as in section 1.4.2, we achieve our goal, which was to obtain a stochastic Schrodinger equation 
for the mirror co-ordinate undergoing a continuous position measurement: 



\i>(t + dt)) = 



+ kQ 2 dt + [Ak{Q)dt + V2kdW}Q 



\m)- 



(1.59) 



This is an equation for the unnormalised state vector, although it may be readily converted to an 
equation for the normalised state vector in the manner of section 1.4.2. 
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Chapter 2 



Quantum Trajectories and 
Quantum noise 



2.1 Stochastic Calculus 
2.1.1 The Wiener Process 

Stochastic calculus is the calculus required to handle differential equations which are driven by 
random processes. We have seen in the previous chapter that quantum trajectories are described 
by differential equations of this nature. We now examine these equations in some detail, what they 
mean, and how to handle them. Our treatment is not intended to be mathematically rigorous, but 
rather is aimed at providing the unfamiliar reader with an understanding of this often rather obscure 
topic. 

A general stochastic equation for a variable x includes both a deterministic increment propor- 
tional to dt, and a random increment |j42f . For the purposes of introducing these equations it is 
best to examine a small time step At and take the infinitesimal limit later. We therefore write the 
general stochastic differential equation as 



where f(x, t) and g(x,t) are arbitrary functions of x and t, and AW is the random increment which 
will take a different value at each time step. We will choose AW to be Gaussian distributed with a 
zero mean (any non-zero mean may be absorbed into the deterministic part in any case), so that 



where a 2 (At) is the variance, being a function of At. Now, how should this variance depend on At? 
Consider the stochastic increment for a time step of At/n, so that the variance is a 2 (At/n). We now 
want to make the stochastic increment over a time step At consistent in the sense that it is given 
by the sum of n of these stochastic increments. Probability theory tells us that when we add an 
arbitrary number of random variables, the variance of the sum is the sum of the individual variances. 



Ax = f(x, t)At + g(x, t)AW, 



(2.1) 




(2.2) 
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The variance for a time step At should therefore satisfy the relation a 2 (At) = na 2 (At/n). This 
linearity condition is satisfied if a 2 (At) is proportional to At. Because scaling the variance may be 
achieved by simply scaling the increment AW , we may set a 2 (At) = At without loss of generality. 
The result is the Wiener increment, with probability distribution given by 



P(AW) = - 1 exp 
V ; V^At 



(AWf 



2At 



(2.3) 



In fact, it is easy to see that any other dependence of the variance on the time step makes no sense 
in the continuum limit. Let us assume that we take a 2 (At) = At a for some a, and calculate the 
variance of the finite increment over a time t in the limit as At tends to zero. The finite increment 
is given by 

W(t)= / dW(t')= lim V AW n , (2.4) 

J{ > 71=1 



where AW n is the increment over the n th interval of duration t/N . The variance of W(t) is then 

N ( A t\ a 

(W(t) 2 ) = lim V(AHU = lim AM — = lim N l ~ a At a . (2.5) 

N^oo — ' N^oo V N ) N^oo 

We see now that it is essential to choose a = 1 to avoid the result being zero or tending to infinity, 
both of which are absurd. 

To be able to manipulate the stochastic equation we must know how the stochastic increment 
acts when raised to some power. For example, if we want to calculate the equation governing x 2 in 



the case that the equation for x is given by Eq.(2.1), then we may write 



A(x 2 ) = (x + Ax)(x + Ax) - x 2 

= 2x Ax + (Ax) 2 = 2xfAt + 2xgAW + f 2 (At) 2 + 2fgAtAW + g 2 (AW) 2 . (2.6) 

Usually the effect of second order terms vanishes in the infinitesimal limit, and the result is that 
given by the chain rule of ordinary calculus: dx 2 = 2xdx. However, note that the expectation value 
of (AH 7 ) 2 is At, which suggests that terms second order in AW may contribute as terms first order 
in At. In fact, it turns out that in the infinitesimal limit we have precisely (dW) 2 — dt. This appears 
especially strange because dW is a stochastic increment while dt is deterministic. We now examine 
why this is so. 

First note that the random variable AZ = (AW) 2 is non-negative, unlike AW which is positive 
and negative with equal probability. In particular, the probability density for AZ is 

P(AZ) = 6 . (2.7) 
V2nAtAZ 

To find out how AZ behaves in the continuum limit, we must examine what happens when we sum 
many of them together. Let us define a variable, Z^(t), which is just such a sum over a time from 
to t, 

N -AZ n /(2t/N) 

*„(.) = X>z. , P(AZ„) = 7 =__, (2.8) 
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and let us denote this variable in the continuum limit (as N tends to infinity) as Z(t). Clearly the 
average value of Z^{t) is constant for any value of N, and is given by (Z(t)) = t. To obtain the 
result dZ = dt, however, we must show that in the infinitesimal limit, when N — > oo, the fluctuations 
of Zjv(t) about the mean value t become smaller and smaller so that Z(t) — t. We could do this by 
calculating the variance of Z^(t), but it is possible to obtain the full distribution for Z(t) which we 
now do using the characteristic function. 

The characteristic function for a random variable is defined as a Fourier transform of the prob- 
ability density. The characteristic function is useful for our purposes here because when adding a 
number of random variables, the characteristic function for the resulting random variable is given 
by multiplying together the characteristic functions for the random variables in the sum jf3|. To 
calculate the probability density for Zpf(t) we therefore first obtain the characteristic function for 
AZ n by taking the Fourier transform of the probability density given in Eq.(2 



AZ„[is-N/(2t)] i 
Xn(s) = I J= , =— d(AZ n ) = m , _= . (2.9) 



y/2n(t/N)AZ n y/(l-ist/(2N)) 
The characteristic function for Z(t) is therefore 



X(s) = lim 

N— »oo 



7(1 - ist/(2N)) 



TV 



N/2 

N*L{ 1+ist {w2 ]] ^ = '"' (2 ' 101 



where we have used the binomial approximation in the second step because N is large. Taking the 
inverse transform we obtain the probability distribution for Z(t) as 

1 r°° 

P{Z{t)) = — e-' ls{z - t Us = 5{z-t). (2.11) 

The distribution for Z(t) is a delta function centred at Z = t, and therefore Z(t) — t. Over any 
time step the integral over the squares of the Wiener increments is equal to the time step, and 
this is true for any time step, in particular an infinitesimal one. As a consequence we may use the 
shorthand rule dW 2 — dt, which is known as the Ito calculus relation. It may be shown that all 
higher powers of the stochastic increment contribute nothing in the continuum limit, and this is also 
true for dtdW @. 

Now that we understand the rules required to handle the Wiener process, we will take a look 
at the solutions to these equations. Before we do however, let us justify our choice of a Gaussian 
random variable for the stochastic increment. Note that we are interested in the continuum limit, so 
that in each arbitrarily small time interval the system gets a kick which is uncorrelated with the kick 
received in the previous time interval. This is a good approximation if the noise driving the system 
is uncorrelated on the time scale of the system dynamics. If this is the case then the central limit 
theorem tells us that if the statistics of the kicks come from any distribution with finite moments, 
when we sum enough of them together the result will be kicks with a Gaussian distribution. We 
can then say that the Wiener process provides a good model for any reasonably behaved source of 
continuous physical noise in which the correlation time is much smaller than the dynamical time 
scale of the system. Physically this means that the noise which is driving the system of interest 
contains very high frequency fluctuations compared to the time scale of the system dynamics. 

In the remainder of this section we examine the solutions of two examples of Eq. ( |2.]| ) which will 
be relevant in later chapters. For our first example we take f(x, t) — kx, where k is some complex 
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number, and take g to be constant, which gives us the Ornstein- Uhlenbeck process j42| 

dx = kx dt + g dW. (2.12) 

This is simply an ordinary linear differential equation with a driving term which is the Wiener 
process. In this case the singular nature of the Wiener increment is not important. To see this we 
may find the solution by changing variables using y = xe~ kt so that 

dy = d(xe~ kt ) = dx e~ kt + x d(e~ kt ) + dx d(e" fc *) = ge~ kt dW. (2.13) 

In this process of changing variables there are no contributions from terms involving dW 2 . Conse- 
quently the solution we obtain is correct for any continuous driving. The solution for y is now 

V(t) = 2/(0) + 5 / e- kt ' dW{t') = 2/(0)4-3 rim V e - fe (»-i)At AW n . (2.14) 
Jo N ^°°^i 

Using x — ye kt we obtain the solution for the Ornstein-Uhlenbeck process, which is 

x(t) = e kt y{0) + ge kt f e' kt ' dW(t'). (2.15) 



(i 



While this solution is no different than that for any ordinary driving term, in this case the integral 
over the driving is a random variable. This random variable is simply a sum over Gaussian variables 
and is therefore itself Gaussian. The variance is easily evaluated to obtain 



lim 1 



N 



e -*(»-i)At AWn 



V e -a*(»-i)A* (AW^) = / e- 2kt ' 

n=l Jo 



dt', (2.16) 



where we have used the independence of the Wiener increments which gives (AW n AW m ) — 8 nm At. 

The Wiener process, due to being of order \fd~i, is not differentiable (it is not of bounded variation) 
even though it is continuous. Nevertheless, it is still possible to write the Ornstein-Uhlenbeck 
equation in the form 

x = kx + ge{t), (2.17) 

where we have used the derivative rather than the differential. It is possible to show that solutions 
to this equation are equivalent to solutions to the differential form if the noise term, e(t), is delta 
correlated. That is, if (e(t)e(t')) — S(t — t') [Q. In particular the moments of J Q e(t')dt' are equal 
to those of J Q dW(t') as required. When written in this form the equation is usually referred to as 
a Langevin equation. The delta correlated (Langevin) noise source is naturally highly singular (its 
variance at any time is infinite), a fact which reflects the non-differentiability of the Weiner process. 
Because spectrum of the noise is the Fourier transform of the correlation function, the spectrum of 
e(t) is flat. For this reason the Wiener process is said to describe white noise. 

While we have concentrated above on a differential equation involving a single variable, the 
Ornstein-Uhlenbeck process is easily generalised to a multi-variable linear equation, which may have 
multiple independent noise sources. We will use equations of this form for our analysis in Chapter^. 

For our second example we consider the time independent linear stochastic equation given by 



dx = [f dt + g dW]x, 



(2.18) 
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in which / and g are constant. This differs from the Ornstcin-Uhlcnbeck process in that the Wiener 
increment now multiplies x, and for this reason it is referred to as linear multiplicative white noise. 
This equation is also readily generalised to multiple variables. In this case / and g become matrices, 
and x the vector of the variables. It is this generalisation which describes linear quantum trajectories, 
and with which we will be concerned in Chapter |^. For the case of a single variable, however, it is 
easily solved. The standard method of solution is to change variables to y = ln(x) , so that 

dy = \dx - ^(dxf = [/ - (g 2 /2)]dt + g dW. (2.19) 

We are now able to integrate directly to obtain 

y(t) = y(0) + [f - (g 2 /2)]t + gW(t), (2.20) 

and the full solution is then x(t) = e y ^\ It is possible to obtain the solution using a different method 
however. This is useful because, as we show in Chapter |^, it may be generalised to equations 
containing multiple variables, and hence to quantum trajectories. In this approach we write the 
equation in the form 

x{t + At) = [l + fAt + g AW]x{t) = e {f - {g2/2))At+aAw x(t), (2.21) 

where the last equality is easily shown by expanding the exponential to first order in At and using 
AW 2 = At. The exponential now acts as an operator which propagates x(t) forward in time by At. 
The solution may therefore be obtained using 

/ N \ 



x(t) = lim J] e</-(» 2 /2))At+ 9 AW„ x(0 ) 



N N 



exp < lim 



(f-(g 2 /2))J2^t + 9j2 AW - 



Af-(9 2 /z)]t+gW{t) 



x(0), (2.22) 

which is indeed the solution obtained previously. This concludes our discussion of stochastic equa- 
tions driven by the Wiener process. For a comprehensive treatment of stochastic equations and 
stochastic calculus the reader is referred to reference f42f . 



2.1.2 The Poisson Process 

The Poisson process consists of a series of instantaneous events occurring at random times, and 
processes of this nature are referred as point processes. For the Poisson process there is a constant 
probability per unit time for an event to occur. We saw in the previous chapter that in the right 
limit the Poisson process may be well approximated by the Wiener process. However, we are not 
always interested in this limit. A general stochastic equation driven by the Poisson process is 

dx = f(x,t)dt + g(x,t)dN, (2.23) 

where this time dN is the Poisson increment. The absence of an event in the time interval dt is 
signified by dN = 0, while an event is signified by dN = 1. The probability for more than one event 
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to occur in any time interval is higher than first order in the time interval, so we need only consider 
these two possibilities. Denoting the probability per unity time for an event by A, the probabilities 
for time interval dt are 

P(dN = 0) = 1 - Xdt , P(dN = 1) = Xdt. (2.24) 

The Poisson increment satisfies dN 2 = dN due to the fact that it takes only the values zero or one. 
The total number of events that have occurred in the interval [0, t] is the integral of the Poisson 
increments over that time. More generally, the integral of a function of time over the Poisson 
increment is 

/ f(t')dN(t') = £ f(t k ) (2.25) 
Ja k=i 

where N(t) is the total number of events up until time t, and tk is the time of the k th event. The 
evolution of x is clearly smooth except when interrupted by discontinuous jumps at the times at 
which dN =1. 



2.2 Quantum Trajectories 

2.2.1 Quantum Trajectories and Generalised Measurements 

In Chapter |l| we introduced master equations by considering a real open quantum system (an optical 
cavity) and motivated quantum trajectories by considering a real continuous measurement (photo- 
detection) performed on this system. However, it is possible to start with the theory of generalised 
measurements, which we also introduced in Chapter ^ and, using this to formulate a continuous 
measurement, obtain both master equations and quantum trajectories in a simple manner. To 
this end consider a generalised measurement in an infinitesimal time interval dt. The non-selective 
evolution caused by this measurement process is 

p(t + dt)=Y^ tt n {dt)p(t)nl(dt). (2.26) 

a 

We now choose a measurement with two possible results, where the measurement operators are 

n Q = l-(j-H + \c c)dt, (2.27) 
Oi = Vdic. (2.28) 

The resulting non-selective evolution is 

p(t + dt) = p{t) - | [H, p(t)]dt + V{c)p{t)dt. (2.29) 

Here T>(c) is the Lindblad super-operator. Recall from Chapter [l] that any master equation may 
be written as a sum over an arbitrary number of Lindblad super-operators, T>(c n ), with different 
measurement operators c„. From the above result we see that each Lindblad operator corresponds 
to a continuous measurement process, alternatively referred to as a decay channel. It is clear that 
this continuous measurement process is a point process, as the probability for the measurement to 
return the result 1 in an interval dt is Ti[c^ cp{t)]dt , being proportional to dt. The result of the 
measurement is therefore usually zero, and is unity only at a discrete series of random times. 
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If we choose to look at the selective rather than the non-selective evolution, then we obtain 
a quantum trajectory rather than a master equation which is the average over all the possible 
trajectories. Clearly the point process we have considered above corresponds to photo-detection 
of the light output from a lossy optical cavity, where a result of zero signifies that there was no 
detection, the result of unity corresponding to a detection. This measurement process may be 
modelled explicitly by an interaction with a two state probe system for time dt, where the probe is 
measured at the end of the interaction p4| . Denoting the two states of the probe as |0) and |1), and 
defining a probe transition operator, a — |0)(1|, we choose the interaction Hamiltonian to be 

H int = ng(ca + ac), (2.30) 

This choice of the interaction Hamiltonian means that the system will cause a transition in the 
probe such that the rate of transition depends upon c. A subsequent measurement of the probe tells 
us whether a transition has occurred, and this provides information about the system. Taking the 
initial state of the probe to be |0) (in fact the interaction is symmetric in the probe state so either 
initial state is acceptable), and evolving the system and meter for a time At, we obtain, to second 
order in At 

Ptot (t + At) = p tot - j-[H totlPtot }At- ^[H tot ,{H tot ,p tot }]At 2 + ... 

= -|[if,p(f)] ® |0) (0| At + {cp® |1)(0| + P c ® |0)(l|) 5 At 

+|(2cpc t ® |1)(1| -Jcp® |0)(0| -pJc® |0)(0|)( 3 At) 2 + ... (2.31) 

where we have written the total Hamiltonian for the system plus probe as -f/tot = H + H m t. At 
the end of the interval we perform a projection measurement onto the probe states |0) and |1). The 
state of the system upon obtaining the result |0) is 

p(t + At) = p(t) - p{t)]At - (l/2)(ct C p + P Jc)(gAt) 2 , (2.32) 

and upon obtaining the result |1) is 

p(t + At) = cp{t)c\gAt) 2 . (2.33) 
Taking the limit as At — > dt, and allowing g to scale as 1/^/At, this is just the result obtained by 



applying the measurement operators introduced in Eqs.( 2.27 ) and ( 2.28 ) 



The set of measurement operators given by Eqs. fl2.27j ) and ( 2.2§| ) is not the only set which will 



generate the correct non-selective evolution. In Chapter pi we noted that applying a unitary transfor- 
mation to the measurement operators leaves the non-selective evolution unaffected. A transformation 
of particular interest is that given by 



\ / 1- h\P\ 2 dt -fj*Vdi \ ( n 



n[ ) v /Vdt i- W\ 2 dt ) v fti 



(2.34) 



This is equivalent to the transformation 

c^c + f3 H -»• H -±iK0*c-Pct). (2.35) 

In the limit in which |/3| 2 3> (crc) this corresponds to homodyne detection, in which the field to be 
measured is first mixed with a coherent field at a beam splitter before being subjected to photo- 
detection (see section 1.4. 2| ). It is also possible to generate the same non-selective evolution using 
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a completely different set of measurement operators, although we are not concerned with these 
here. For an example see reference p5[ . In the following section we will see how we may formulate 
linear and non-linear stochastic equations to describe the selective evolution, or quantum trajectories, 
generated by generalised measurement operators. 



2.2.2 Linear and Non-linear Formulations 



We have seen in the previous section how the non-selective evolution produced by a Lindblad super- 
operator may be unravelled by the continuum limit of a generalised measurement process. In each 
time interval, At, a generalised measurement is performed, and the manner in which the state of 
the system is altered depends on the result of the measurement. This is repeated in the next time 
interval and so on. We want to examine now how we might simulate this process. In particular we 
are interested in (i) the probabilities of the possible measurement sequences, and (ii) the state of the 
system during and after the measurement process. In the following discussion we use the approach 
of Wiseman [|6| . 

Let us assume that the initial state is pure, and denote it by In this case the state of the 

system given the result r is 

Q r (Ai)|^(*)) 



m + At)) 



where P r (t) is the probability for obtaining the result r, being 

p r (t) = (ip(t)p r (At)n r (At)\i>(t)). 



(2.36) 



(2.37) 



To simulate this procedure, we could calculate the probabilities for the possible results using using 



Eq.( 2.37 ), and choose one at random using these probabilities. We then calculate the new state using 
Eq.(2.36), and repeat the procedure for this new state. The probability of obtaining a particular 
final state at the end of a sequence of measurements is then given by the probability that we actually 
obtain that state at the end of our simulation. Call this method A. We may formulate this simulation 
procedure as a stochastic differential equation which turns out to be non-linear. Applying this to the 
measurement process described by the measurement operators Q' and we have two possibilities 
at each infinitesimal time step, and we may use the Poisson process introduced in the previous 
section. By considering the change in the state vector for each of the two possibilities the stochastic 
equation is readily found to be 



\ip(t + dt)) = 




-l\dN 



((c f +p*)( c + /3)) 
(etc) etc (c^ + /3*c) 



/3*c- 



dt 



(2.38) 



It is possible to use another method to simulate the measurement procedure. Instead of choosing 
the outcomes at random using the correct probabilities, we choose them with fixed probabilities 
independent of the state. Following reference ]46]] we will refer to these as the ostensible probabilities. 
In this case the actual probability for obtaining a given final state is not the probability that we 
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obtain that state at the end of the simulation. However, we can preserve this information by omitting 
to normalise the state vector at each step. In this case the state after a measurement is given by 

\4>{t + At)) = fl r (At)\tp(t)), (2.39) 



where we use the tilde to denote an unnormalised state. Re- writing Eq.( 2.37 ) we see that the 
probability for obtaining the final state is simply the norm of that final state: 

P r (t) = $(t + At)\ip(t + At)). (2.40) 

Following this process through a series of measurements we see that the probability of obtaining 
a particular final state is just the norm of that final state. This simulation process is now linear 
because the final (unnormalised) state is just the initial state transformed by the action of a series of 
linear operators. Call this method B. This method is not quite sufficient for our purposes, however. 
Note that in method A we divided the final state by the square root of the probability of obtaining 
that state, and this resulted in normalisation. In order to take the infinitesimal limit to obtain 
a continuous measurement (and in our case, a Poisson process), we will need to maintain this 
procedure of dividing the final state by the probability with which we choose it, although this time 
the probability will be independent of the initial state, as in method B. Thus we choose among the 
possible final states at random, using probabilities chosen to be independent of the initial state, this 
time taking the final state to be 

, ~, . v . n r (At)\ib(t)) 

m + At))= 2.41 

where A r (t) is the ostensible probability with which we (artificially) choose the outcome r. Call this 
method C. Clearly the true probability of obtaining the result r is still contained in the norm of the 
final state. In this case it is 

P r (t) = (ip(t + At)\i>(t + At))A r (t), (2.42) 

and the true probability for a given final state at the end of a series of measurements is 

P(t + r) = m + r)|# + r))P(i), (2.43) 

where P(t) is the ostensible probability that we obtain the given final state at the end of our 
simulation procedure. This method will allow us to write a linear stochastic equation for the quantum 
state. 

To derive the linear stochastic equation describing photo-detection performed on the light emitted 
from a damped optical cavity, the measurement operators are fio and f2i with c = y/ja. Choosing 
the ostensible probabilities Ai(t) = jdt, and Ao(t) = 1 — Ai(t), to give a Poisson process with a 
constant average rate of events 7, the infinitesimal increment of the state at time t becomes 

\Mt + dt)) = = am , ( 2 .44) 



\Mt + dt)) = "oWlffl = [1 - Igig- 2) dt ] $(t)). (2.45) 
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The equivalent stochastic equation is therefore 

\Mt + eft)) = [1 - (|# + 2 a t a _ 2) dt + ( a _ 1)^] |&(t)). (2.46) 

We will examine how solutions to this equation may be found in the next chapter. 

To derive a linear stochastic equation corresponding to homodync detection we use the measure- 
ment operators tt' ft[ and choose the ostensible probabilities Ao(i) = 1 — |/3| 2 (ft, and Ai(t) = |/3| 2 cft. 
The infinitesimal increment of the state at time t becomes 

n[(dt)\ij(t)) 



\Mt + dt)) = 
and the stochastic equation is 
\Mt + dt)) 



c 

1 + - 



P 



W)) 



l-{iH+yc + 0*c)dt |#)), 



!&(*)>• 



j-H+^c c + (3*cj dt + ^dN 



(2.47) 
(2.48) 

(2.49) 



Note that this is much simpler than its non-linear equivalent Eq.(2.3S). To obtain an equation for 
homodyne detection with a local phase set so as to measure the phase quadrature, we set c = ^/ja 
and = ?a/7/3 where (3 is real. We then take the Wiener process limit of the Poisson process as 
discussed in Chapter [j]. The result is the linear equivalent of the non-linear stochastic equation 
derived for homodyne detection in that chapter, 

\$i(t + dt)) = 1- UH+^a ajdt-iay/jdW \4>i(t)}. (2.50) 

More generally, given a master equation of the form Eq.( 1.2S| ), a non-linear unravelling in terms of 
the Wiener process is given by 



\tp(t + dt)} = 



1 - IHdt 
h 



N 

E(K C « ~ ( C n>C« + i(ct)<C„»<ft + ( C „ - (c n ))dW n (t) 
n=l 



\mh (2-51) 



where the dW n are independent Wiener increments. A linear unravelling for the same master 
equation is given by 



\i>(t + dt)) 



1 



N 

-Hdt - ^(\c\ lCn dt - c n dW n (t)) 
n=l 



(2.52) 



In the linear case, the true probability measure for the system to have evolved to a particular state 
at time t is given by (47) 

(iP(t)\ip(t)) w dP w , (2.53) 

where dP w is the Wiener measure (the measure for the ostensible probabilities). That is, it repre- 
sents integration over the joint probability density for all the random variables that appear in the 
expression for \ip(t)) w . It follows therefore that moments of system operators calculated with the 
equivalent master equation at time t are given by the expression 

(o) = [ (m\o\m)^p w . 



(2.54) 

where O is the system operator in question, and dP w represents integration over all possible values 
of the random variables. In Chapter ^ we will show how expressions for final states and their 
corresponding probability densities may be obtained. 
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2.3 Quantum Noise 

So far we have considered the dynamics of an open quantum system using the Schrodinger picture. 
However, it also is possible to use the Heisenberg picture. In this case the system operators obey 
equations of motion driven by a quantum version of the Wiener process. We now examine this 
formulation which we use for our analysis in Chapter ^. 

We will consider here the specific example of a lossy optical cavity damped through one mirror, 
and will therefore require the part of the external field propagating along the axis of the cavity, 
and on one side of the cavity only. That is, a one-dimensional field in half-space. For a detailed 
discussion of the representation of this field in terms of mode operators the reader is referred to 
appendix A, where we examine the input-output relations for the cavity. Here we need merely know 
the Hamiltonian which gives the coupling between the cavity mode and the field mode operators, 
and this is 

/■oo 

H = H ays -ig (b\k)a- a%(k)) dk + H RcM . (2.55) 

J — oo 

Here the Hamiltonian for the cavity mode is denoted by -ff S ys, and that for the external free field 
is denoted by -ffficid- In the expression for the interaction Hamiltonian, g is the coupling constant 
between the mode and the external field, a is the mode annihilation operator and b(k) is the field 
annihilation operator for the field mode with angular wave number k. In order to arrive at this 
interaction, which clearly models the exchange of photons between the cavity mode and the field, 
the rotating wave approximation has been employed, and this is discussed in appendix A. Using the 
above Hamiltonian we may calculate the equations of motion for an arbitrary system operator c, 
and for the field mode operators b(k). These are 

c = -^[c,H sys ]+g[c,a^(J ' b(k,t)dkj - g ^J°tf(k,t)dkj [c,a], (2.56) 

b(k,t) = -ickb(k,t) + ga. (2.57) 
Solving the equation for the mode operators in terms of the system operator a we have 

b(k,t) = b(k,0)e~ tckt +g f e~ lck{t - tr) a(t')dt'. (2.58) 

Jo 

Substituting this back into the equation of motion for c we obtain 

c = -~[c, H sys ] + [c, a] (z a t + V7ftL(t)) - (£a + VtM*)) « + ], (2-59) 
where we have defined 7 = Trg 2 /h 2 , and the operator 

M*) = \ — / b(k,0)e- ickt . (2.60) 
V c J_ oc 

We see that the equations of motion for the system operators are driven by 6i n (i), which may therefore 
be referred to as the field 'input' to the system. For a detailed explanation of the relationship between 
the initial field state and the input and output fields (being those parts of the field travelling to and 
from the system respectively), the reader is referred to appendix A. We have scaled b- m (t) specifically 
so that it satisfies the commutation relations 

[b in (t),bl(t)]=t(t-t'), (2.61) 
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and which is a direct result of the commutation relations of the field mode operators. Note that the 
expectation value of this 'input' operator at all future times is specified completely by the state of 
the field at t — 0. If the initial state of the field is the vacuum, then (b in (t), b in (t)) — 0. Using the 
commutation relation Eq.(2.61), this gives the anti-normally ordered correlation function as 

(b in (t),bl(t))=5(t-t>). (2.62) 

We see that the terms that drive the equations of motion for the system are delta correlated! They 
are therefore a quantum analogue of the classical Langevin noise source that we introduced in 
section 2A . It is useful for us to put this equation in the form of an explicit quantum analogue of 
the Ito stochastic increment. From this we will be able to derive the master equation of section F3 
and in so doing show that this quantum Langevin formalism is equivalent to the master equation. 
Defining the stochastic increment operator as dB{t) = b m (t)dt, we have 



(/ dB(t') / dB\t")) 



(b- m (tf)bUt"))dtfdt" = t, 



(2.63) 



and in the infinitesimal limit the (quantum) Ito relation dBdB^ = dt. Further discussion of these 
noise operators may be found in reference |4S[ ] . 

To derive the Ito equation we move into the interaction picture with respect to -fffioid, and write 
the resulting interaction Hamiltonian in terms of b m {t). We have 

H = H sys + ih^(abl(t) - b in (t)a^). (2.64) 

The evolution of the system operator c is therefore 

C(t + dt) = e - i («sys/ft)dt+\/7(adB t (t)--dB(0a t ) c (i) e -i(ffs.v S /?i)<ii- v ^(adS t (t)-dB(t)a t )_ (2.65) 

Expanding this to first order in t, using the Ito relation and the fact that the stochastic increment 
operators commute with the system operators at time t, we obtain the Ito equation 



c(t + dt) = - — [c, H sys ]dt + y(atca - \a) ac - \ca}a)dt + y/j[c, a dB\t) - dB(t)a f ] 



(2.66) 



From this we may derive the master equation. We do this by calculating the expectation value of 
the derivative of c using the Ito equation, and by equating this to that calculated using the density 
operator: (c) = Tr [cp\. Using (dB) — 0, the expectation value of the derivative of c from the Ito 
equation is 



([c,-Hsys]) +7((a f ca) - \{a) ' ac) - Uco) a)) 



= Tr 



+ r y(a i ca — ^a)ac — ^ca)a) I p 



- = [H sys ,p\+i(apa) 



it i t - 
±a ] ap — 2 pa a. 



= Tr 
= Tr [cp] . 

This is true for any operator c, and so we obtain the master equation 

P = --[H sys ,p] +7(apa f - \a*ap- \pa)a). 



(2.67) 



(2.68) 
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Every master equation has an equivalent quantum Langevin equation and vice versa. 

We may extend the quantum Langevin equation to include other kinds of input. In our analysis 
above we have set the initial state of the field operators b(k) to the vacuum. An almost trivial 
extension is to take the initial state of each mode to be coherent, so that b(k, 0)\(3(k))k — (3(k)\f3(k))k- 
Using a displacement operator to transform each of the initial states to the vacuum jls]. [32 
transformed input operator becomes 

b' in (t) = b in (t) + p(t), (2-69) 

where (3(t) is a c-number, being an appropriately scaled Fourier transform of (3(k). As we are free to 
choose the /3(h), we are free to choose any (3{t). For a coherent state input (3(h) is a delta function 
about k — uj c /c, and consequently f3(t) is constant. In Chapter || we will want to consider the 
input from a laser. While this is basically coherent, it will naturally contain some noise so that the 
amplitude and phase will fluctuate. We can include this classical fluctuation by allowing the modes 
on either side of k = uj c /c to have a non-zero coherent amplitude. Assuming small fluctuations about 
the mean, we may write the input as 

P + 6x(t) +iSy(t) + b iD (t), (2.70) 

where (3 is the coherent amplitude of the laser, 5x(t) are fluctuations in the amplitude quadrature and 
8y(t) are the fluctuations in the phase quadrature. The autocorrelation functions for the fluctuations 
may be left arbitrary, and can be tailored to match the noise on the laser used in any particular 
realisation. 
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Chapter 3 

Evolution Operators for Linear 
Quantum Trajectories 

In this chapter we present a method for obtaining evolution operators for linear quantum trajec- 
tories, and apply this to a number of physical examples of varying mathematical complexity. In 
particular we use as examples quantum trajectories describing the continuous projection measure- 
ment of physical observables. Using this method we calculate the average conditional uncertainty 
for the measured observables, being a central quantity of interest in these measurement processes. 



3.1 Introduction 

We have seen in the previous two chapters that real continuous measurements may be described by 
stochastic equations for the quantum state vector. That these may be written in a linear form has 
been known in the mathematical physics literature for some time pS| , but has only fairly recently 
seen exposure in the physics literature [ ff7[ 50, ^3), where it has been common to use non-linear 



stochastic equations [|4], [h], [52). The advantage of writing master equations as LSE's, rather 
than the more familiar non-linear version, is that in certain cases it has been found that explicit 
evolution operators corresponding to these equations may be obtained in a straightforward manner. 
Naturally the non-linear form does not admit a solution in terms of evolution operators, as these 
operators are by definition linear. However, as far as we are aware, the only method that has been 
used previously to obtain evolution operators for quantum trajectories driven by the Wiener process 
is to choose an initial state which allows the stochastic equation for the state to be written as a 
stochastic equation for an eigenvalue, or which simplifies the action of the evolution operator [^7], [53| . 
In this chapter we present a more general method for obtaining explicit evolution operators for these 
equations which makes no reference to the initial state. We also show that this method may be 
used to derive evolution operators for quantum trajectories driven by the Poisson process, although 
we include this essentially just to show that a unified treatment may be used for both cases; the 
resulting operators do not extend significantly those which have been derived elsewhere |54|]. We 
note also that various authors have found solutions to certain non-linear equations for a single 
particle |p3, fia, bq, BTj . These require however that the initial state takes a Gaussian form, and do 
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not easily generalise to other initial states due to the non-linearity of the equations. For treatments 
of non- linear stochastic state- vector equations see also |37], Q . 

Naturally the resulting evolution operators contain classical random variables. The complexity of 
the stochastic equations which govern these classical random variables depends upon the complexity 
of the commutation relations between the operators appearing in the LSE. If the complexity of 
the commutation relations is sufficiently high then the stochastic equations governing the classical 
random variables become too complex to solve analytically. Nevertheless, even if this is the case, 
the form of the evolution operator provides information regarding the type of states produced by 
the LSE, and the problem is reduced to integrating the classical stochastic equations numerically. 
We also note that the solution to an LSE provides additional information to that contained in 
the solution to the equivalent master equation, because it gives the state of the system for each 
trajectory. For example, the variance of a system operator may be calculated for each final state 
(ie. for each trajectory), and this is referred to as the conditional variance as it is conditional upon 
the results of the measurement. The overall average of these variances may be then be calculated. 
The solution to the master equation allows us to calculate only the variance which is obtained by 
first averaging the final states over all trajectories, which is, in general, quite a different quantity. 

In the following treatment of stochastic equations involving the Wiener process we use as ex- 
amples LSE's corresponding to the continuous measurement of physical observables. A term of the 
form 

(> = ■■■ -k[0,[0,p}}--- (3.1) 

in a quantum master equation for the evolution of a density matrix, p, for a quantum system S, de- 
scribes a continuous projection measurement of an observable O of S. We refer to this measurement 
process as a continuous projection measurement because in the absence of any system evolution, the 
sole effect of this term is to reduce the off-diagonal elements of the density matrix to zero in the 
eigenbasis of that observable. That is, it describes, in the long time limit, a projection on to one of 
the eigenstates of the observable under observation. The rate at which information is gained regard- 
ing the observable is determined by k which is a positive constant. Recall that in Chapter |j] we gave 
a specific example of a continuous position measurement which took this form. Various examples of 
measurements of this kind may be found in references [^8, 45, 5^, ^0|. If, in addition, the observable 



commutes with the Hamiltonian describing the free evolution of the system under observation, then 
the free evolution does not interfere with this process of projection, and the measurement is referred 
to as a continuous Quantum Non-Demolition (QND) measurement |24|, pl|| . 



3.2 The Wiener Process: General Method 

We will explicitly treat here LSE's which contain only one stochastic increment. However it will be 
clear that this treatment may be easily extended for multiple stochastic increments. Let us write a 
general LSE with a single stochastic increment as 

d\rp(t)) = (Adt + B dW(t))\ip(t)). (3.2) 

In this equation A and B are arbitrary operators. We will see that the complexity of the evolution 
operator will depend upon the complexity of the commutation relations between A and B. 
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As a first step in obtaining an evolution operator for the LSE in Eq.(3.2) we rewrite it in the 
form 

\${t + dt)) = eS A -^ 2 ^ dt e Bdw ^m)) = e Adt e Bdw ^\i>{t)), (3.3) 

were we have defined A = A— B 2 /2. It is easily verified that this is correct to first order by expanding 
the exponentials to second order and using the Ito calculus relation dW(t) 2 = dt. To first order the 
state at time t + At is therefore 

\rj>(t + At)) = e AAt e BAW W\4>(t)), (3.4) 
so that the state at time t may be written 

N 

= Jim H(e AAt e BAW ") |<A(0)), (3.5) 

where 

AW n = / dW{t), (3.6) 

and N — * oo as At — > so that NAt = t is always true. We have included the subscript 'w' to 
remind us that the state is dependent on the Wiener process. To complete the derivation of the 



evolution operator we must take the limit in Eq.(3.5). To do this we must combine the arguments 



of the exponentials which appear in the product, so that we may sum the infinitesimals. We will 
choose to do this by first repeatedly swapping the order of the exponentials containing the operator 
A with those containing the operator B. The simplest case occurs when A and B commute so 
that the problem essentially reduces to the single variable case, and we treat this in Sec. |3.2.1 , The 



simplest non-trivial case occurs when the commutator [A, B], while non-zero, commutes with both 
A and B, and we treat this in subsection [3.2.2| . In the final part of this section we examine a more 
complicated example in which the commutator [A, B] does not commute with either A or B. 

3.2.1 A QND Measurement of Photon-Number 

The mathematically trivial case occurs when A and B commute. A non-trivial physical example 
to which this corresponds is a QND measurement of the photon number of a single cavity mode. 
Denoting the annihilation operator describing the mode by a, the cavity field Hamiltonian is given 

by m 

H = ftu(ata+~), (3.7) 

in which u> is the frequency of the cavity mode, and the observable to be measured is O = a^a. With 
this we have 

A = -iui(a!a+ -) - 2fc(a f a) 2 , (3.8) 
B = \/2fca f a, (3.9) 
in which k is the measurement constant introduced in Eq. ( [3.1| ) . As A and B commute the exponen- 



tials in Eq.(3.5) combine trivially and we obtain 

B^AWn 



|^)) w = A lim o e AJVAt exp 



|V(0)> =e At e BW{t) |^(0)). (3.10) 
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As the Wiener process, W(t), is a sum of independent Gaussian distributed random variables, W n , 
it is naturally Gaussian distributed, the mean and variance of W(t) being zero and t respectively. In 



a particular realisation of the stochastic equation Eq.(3.2), the Wiener process will have a particular 
value at each time t, and as we mentioned above, the set of all these values corresponds to the 
trajectory that is taken by that particular realisation. The fact that to obtain the state at time t we 
require only the value of the Wiener process at that time means that we do not require all the tra- 
jectory information, but just a single variable associated with that trajectory. For more complicated 
cases, in which the operators do not commute, we will find that other variables associated with the 
trajectory appear in the evolution operator. 

As the situation we consider here is a QND measurement, the phase uncertainty introduced by 
the measurement of photon number does not feed back to affect the measurement, so that the result 
is simply to decrease continuously the uncertainty in photon number, and the state of the system as 



t tends to infinity tends to a number state. If we denote the evolution operator derived in Eq.(3.1C) 
by V(t), and start the system in an arbitrary mixed state p(0), then at time t the normalised state 
of the system may be written 

_ V{t)p{Q)VHt) 
P[ ) ~ Tr{V(t)p(0)Vi(t)}- ( ' 

As V(t) is diagonal in the photon number basis, we only require the diagonal elements of the initial 
density matrix to calculate moments of the photon number operator. Denoting the diagonal elements 
of the initial density matrix by p n , and the diagonal elements of V(t)V'(t) by V n , the variance of 
the photon number, for a given trajectory, is given by 



_ En" PnVn (En n PnVn) , . 

a n( l )w - ^ v —= 77TT- \ a - LZ ) 



The uncertainty in our knowledge of the number of photons is the square root of this variance. 
Averaging this uncertainty over all trajectories therefore tells us, on average, how accurately we will 
have determined the number of photons at a later time. To calculate the value of the uncertainty 
for each trajectory, averaged over all trajectories we must multiply cr n (t) w by the probability for 
each final state and average over all the final states. The probability measure for the final states, 
pit), is given by the Wiener measure multiplied by the norm of the final state, Tr {V(t)p(G)V^ (£)}■ 
This probability measure is not in general Gaussian in W(t), but a weighted sum of Gaussians, one 
for each n. Performing the multiplication, we obtain the average conditional uncertainty in photon 
number as 

(<J n {t)w) = / . n(m - n)p n p m V n Vm dP w , (3.13) 

J V nm 

in which 



dP w = -^e^/^dW. (3.15) 

We note that (cr n (t) w ) may be written as a function of r = kt, being the time scaled by the 
measurement constant. Hence, as we expect, increasing the measurement time has the same effect 
on (a n (i) w ) as increasing the measurement constant. We evaluate (<J n (T) w ) numerically for an initial 
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Conditional Uncertainty for a QND Measurement of n 



Initial Photon Number Distributions 




C - 
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X 

Figure 3.1: The conditional uncertainty in photon number averaged over all trajectories, (a„(i) w ), 
is plotted here against the dimensionless scaled time, t = kt. The dotted line corresponds to an 
initial coherent state, and the solid line to an initial thermal state. Both initial states were chosen 
to have a 2 = 20, giving the thermal state a mean photon number of (n) — 4, and the coherent state 
a mean photon number of (n) = 20. The photon number distributions for the two initial states are 
displayed in the inset. 



thermal state, and an initial coherent state, and display the results in figure 3.1. We have chosen 
the initial states so that they have the same uncertainty in photon number, with the result that the 
mean number of photons in each of the two states is quite different. The results show the decrease 
in uncertainty with time, which is seen to be only weakly dependent upon the initial state. 

3.2.2 A Measurement of Momentum in a Linear Potential 

The simplest mathematically non-trivial case occurs when the commutator between A and B, while 
non-zero, is such that it commutes with both A and B. A physical situation to which this corresponds 
is a continuous measurement of the momentum of a particle in a linear potential. If we denote the 
position and momentum operators for the particle as Q and P respectively, then the Hamiltonian 
is given by 

H=^-P 2 -FQ 1 (3.16) 
2m 

in which m is the mass of the particle and F is the force on the particle from the linear potential. 
In this case we have 

A = (^-2^ + fg, (3.17) 



B = V2kP, (3.18) 



in which k is again the measurement constant. 

Returning to Eq.(3.5) we see that to obtain a solution we must pass all the exponentials containing 
the operator B to the right through the exponentials containing the operator A. In order to perform 
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this operation wc need a relation of the form 



e B e A = e A e c . 



(3.19) 



For the present case the required relation is simply given by the Baker-Campbell-Hausdorff for- 



mula 62, 23 



e B e A = e A e B e-^ 



(3.20) 

Using this relation to propagate successively all of the exponentials containing B to the right in the 
product in Eq.(3J3) we obtain 



N 



no 



AAl RAW, 



) = exp [AN At] exp 



N 



exp 



N 



-[A, B]At^2(n-T)AW„ 



(3.21) 



All that remains is to calculate the joint probability density for the random variables. The first is 
the Wiener process, and the second is 



Y(t) = 



lira AtYVrc- 1)AW„ = / t'dW(t') 



(3.22) 



Ai 

Clearly these are both Gaussian distributed with zero mean and all that we require is to calculate 
the covariances (Y(t) 2 ) and (W(t)Y(t)) to determine completely the joint density at time t. Using 
(AW n AW n ) — S nm At these quantities are easily obtained: 



(Y(t) 2 ) 



71=1 



t u dt' = t73, 



(W(t)Y(t)) = lim J2{(n - I) At) At = f t' dt' = t 2 /2. 

The state at time t, under the evolution described by the stochastic equation, is therefore 

\^(t)) w = e At e BW ^e-^ B ^\m), 
where the joint probability density for W and Y at time t is given by 



(3.23) 
(3.24) 

(3.25) 



P W (W,Y) = 



27rf 2 J 



exp 



6 



--W 2 - ^Y 2 + ^rWY 



6 

t 2 



Note that to obtain the probability density for the final state, this must be multiplied by the norm 
of the state at time t. Returning to the specific case of a particle in a linear potential, we may 
now obtain results for various quantities of interest. Writing the evolution operator in terms of the 
momentum and position operators we have 



exp 



exp 



2k(PW(t) + FY(t)) |-0(O)> . 



(3.26) 



Using the Zassenhaus formula [ |63| to disentangle the argument of the first exponential we may 
rewrite this in the more convenient form 



|V>(t)) w = exp 



iF 



-Qt 



exp [i](-P 2 t - PFt 2 - F 2 t 3 /3)] exp \V2k(PW(t) + FY(t))] \if)(0)), (3.27) 
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in which rj — {^Em ^ ^k). For those not familiar with the Zassenhaus formula, it is complementary 
to the BCH formula. While the BCH formula shows how to write a product of the exponentials 
of two operators as an exponential of the sum of the operators and their commutator (or in more 
complicated cases repeated commutators of the operators), the Zassenhaus formula shows how to 
write an exponential of the sum of two operators as a product of the exponentials of the operators 
and repeated commutators. 

Let us first take an arbitrary initial state, writing it in the momentum eigenbasis so that we have 

/oo />oo 
9(p)\p)dp , / \^(p)\ 2 dp=l. (3.28) 



Using Eqs.( [3.27 ) and (3.28) in Eq.(2.54) to calculate the moments of p given by first averaging over 



all trajectories (that is, the moments which would be given by the equivalent master equation) we 
readily obtain 

(p(t) n ) = ((p(0)+Ft) n ). (3.29) 

In particular, for any initial state, |^(0)), the average value of the momentum at time t, (p{t)), is 
simply shifted from the initial value by the impulse Ft. The variance of the momentum at time i, 
CTp(t) = (p(t) 2 ) — (p(t)} 2 , remains equal to its original value. That is, the uncertainty introduced into 
the position of the particle by the momentum measurement does not feed back into the momentum, 
even though the momentum does not commute with the Hamiltonian. This is because while the 
momentum determines the position at a later time, the converse is not true. These results for the 
moments are easily checked using the equivalent master equation. 

Now let us consider the conditional variance of the momentum at time t averaged over all 
trajectories. In the previous section we calculated the conditional uncertainty, being the square root 
of the variance, and averaged this over all trajectories. Here however, we will find that the conditional 
variance is independent of the trajectory taken, and depends only on the measurement time. This 
will also be true of the example which we will treat in the next section. In this case clearly it does not 
matter if we first average the conditional variance over the trajectories, and then take the square root, 
or if instead we average the conditional uncertainty, because the averaging procedure is redundant. 
However, in general the two procedures are not equivalent. We will denote the conditional variance 
by ((7p(£) w ). As the uncertainty in position does not feed back into the momentum, we expect that 
this variance should steadily decrease to zero. This is because during a trajectory our knowledge 
of the momentum steadily increases so that the distribution over momentum becomes increasingly 
narrow. To perform this calculation we take the initial state to be the minimum uncertainty wave- 
packet given by the ground state of a harmonic oscillator of frequency u>. The average values of the 
position and momentum of the particle are both zero in this state and the respective variances are 

= ^ S ^>- (P 2 ) = ^ s oj(0), 
and in momentum space the state may be written 

1-0(0)) = (—!—) 4 re- p2/(2mh ^\ P ) dp. (3.30) 
\nmhujj J_ 00 

The moments of momentum for each trajectory are given by 

(p-% = WyW'»w, (3.31) 
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and we calculate the first and second to give cr%(t) w — (p(t) 2 ) w — (p(t))^. We obtain 

This is independent of W and Y and hence independent of the trajectory. It is therefore unnecessary 
to average over the final states. Indeed (<Jp(t) w ) decreases steadily from the initial value to zero 
as t — ► oo as we expect from the discussion above. This means that while the average value of 
momentum is determined by the measurement results, the error in our estimate of the momentum 
at time t is not. 

3.2.3 A Quadrature Measurement with a General Quadratic Hamiltonian 

We now consider an LSE in which the commutator [A, B] does not commute with either A or B. 
As in the previous example, let P and Q be, respectively, the canonical momentum and position 
operators for a single particle so that they obey the canonical commutation relation [Q, P] = ih. 
With this definition we will take A and B to have the following forms: 

A = aP 2 + 7 Q 2 +£QP + rjP + (Q, (3.33) 
B = kQ + kP, (3.34) 

where a, 7, 77, £, k and n are complex numbers. This example applies to an optical mode of the elec- 
tromagnetic field, including classical driving and/or classically driven subharmonic generation |64| 
and for which an arbitrary quadrature is continuously measured. It also applies to the situation 
of a single particle, which may feel a linear and/or harmonic potential, and which is subjected to 
continuous observation of an arbitrary linear combination of its position pl[ and momentum. 

To obtain an evolution operator for the LSE with this choice of the operators A and B, we 



require, as before, a relation of the form given by Eq.( 3.1S ). To derive this relation we proceed in 
the following manner. 

First we may use the Baker-Campbell- Hausdorff expansion j23|, or alternatively solve the equa- 
tions of motion given by dB/de = [A, B], to obtain an expression for e tA Be~ tA . The result is 



e - eA eBe eA = e/i(e)Q + ef 2 (e)P + e/ 3 (e), (3.35) 



in which 



/i(e) = i(-2/s7 + *0S + kC, (3.36) 

f 2 (e) = i(2fea-K€)S + KC, (3.37) 
A 

/ 3 (e) = ^{k^ + 2kaC-Ka-2^r,)[C-l] + j{kn + K0S. (3.38) 



In these expressions C — cosh(i?iAe), S — smh(ifi\e) and A = \/ £, 2 — Aaj. Using the relation 

e~ eA f(eB)e eA = f(e~ eA eBe eA ), (3.39) 



we obtain from Eq.(3.35) 



e -eA £ eB e eA = ^Q+ef^P fM') . ( 3 . 40 ) 
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Multiplying both sides of this equation on the left by e eA we obtain a relation of the form e e 
e eA e eD ^' , as we require. 



We see from the above procedure that the relation in Eq.( 3.19 ) may be obtained so long as 
a closed form can be found for the solution to the operator differential equation dB/de = [A,B]. 
Clearly this is straightforward if this equation is linear in B, which is true in the example we have 
treated here, and is sometimes possible in cases in which the equations are non-linear. In addition, 
for this example we also require the BCH relation in the form 

e A e B = e A+B+ i^ B l (3.41) 

This is so that we can sum up in one exponential the operators that result from swapping e AWnB 
and e nAtA . 

Using the expressions derived above, with the replacements e = nAt and e = AW n , for each n 
from 1 to N, by repeatedly swapping the exponentials containing B with those containing A as in 
the previous example, we obtain 



N 

lim l[(e AAt e BAW ") = e ^ e ^(t)Q+x 2 (t)P e x 3 (t) j%z(t) > (3 42) 

71=1 

in which the classical stochastic variables Xi and Z, are given by 

Xi(t) = I ' fi{t')dW{t'), (3.43) 
Jo 

Z(t) = I ' h{t')X 2 {t')dW{t')- I ' h^Xx^dWit'), 



where the expressions for the fi are given above, and the integrals are Ito integrals. The Xi are 
Gaussian distributed with zero mean, and their covariances are easily calculated as in the previous 
example: 

{Xi{t)X j {t)) = [f i {t')f j {t')dt'. (3.44) 
Jo 

In addition, the two-time correlation functions for these variables are also easily obtained analytically. 
In particular we have 

/•min(t,r) 

(X^Xjir)) = fiV)fjV)M. (3.45) 

Jo 

However, Z(t) is not Gaussian distributed. We are not aware of an analytic expression for this 
variable, so that its probability density may have to be obtained numerically. We note in passing, 
however, that in some cases double stochastic integrals of this kind may be written explicitly in 
terms of products of Gaussian variables [Q . We note also that Z determines only the normalisation 
of the final state, and not the state itself. The normalised state at time t is therefore independent 



of Z, and we discuss the consequences of this in subsection |3.2.4 
We may now write the state at time t as 

„AtX 1 (t)Q+X 2 (t)P x 3 (t)+mz(t) 



|^(t)) w = e At e A lW y+A 2W ^ e A 3W +^ W |^ (0 ^ (3 46) 



Hence even though values for averages over all trajectories may in general have to be calculated 
numerically, the evolution operator provides us with information regarding the type of states that 
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will occur at time t. In particular, if the initial state is Gaussian in position (and therefore also 
Gaussian in momentum) , then as each of the exponential operators in the above equation transform 
Gaussian states to Gaussian states, we see that the state of the system remains Gaussian at all 
times. The mean of the Gaussian in both position and momentum change with time in a random 
way determined by the values of the stochastic variables. 

We will shortly consider a particular example: that of a harmonic oscillator undergoing a contin- 
uous observation of position, and use this evolution operator to calculate the conditional variance for 
the position at time t. We will take the initial state to be a coherent state, which is a Gaussian wave 
packet. This conditional variance does not depend upon the trajectory, but simply upon the initial 
state and the measurement time, as indeed we found to be the case for the momentum measurement 
in section |3.2.2 . 



Let us first show that for an initial coherent state the conditional variance of any linear com- 
bination of position and momentum is independent of the trajectory for all of the cases covered 



by the evolution operator in Eq.(3.46). To do this we must calculate the effect of this evolution 
operator on a coherent state. Clearly the effect of the right-most exponential operator is at most to 
change the normalisation, which effects neither the average values of position and momentum, nor 
the respective variances. The effect of the next exponential, being linear in P and Q, is calculated 
in appendix We find that it changes the mean values of the position and momentum, and al- 
ters the normalisation, but the state remains coherent in that the position variance (and hence the 
momentum variance) is unchanged. Finally, the effect of the exponential quadratic in P and Q is 
calculated in appendix |^. We find that this operator modifies the position variance. However, as 
the operator does not contain any stochastic variables, and as the manner in which it changes the 
position variance is independent of the mean position and momentum, we obtain the result that the 
effect on the position variance, and hence the variance of any linear combination of position and 
momentum, is trajectory independent. 

Let us now consider a harmonic oscillator in which the position is continuously observed. This sit- 
uation has been analysed by Belavkin and Staszewski using the equivalent non- linear equations J32| . 
The operators A and B in this case are given by 



—i \ 2 { —imu) 



,2 



a = u^r ~ 2fc J Q ' (3 - 47) 

B = V2kQ, (3.48) 



in which m is the mass of the particle, uj is the frequency of the harmonic oscillation, and k is 
the measurement constant for the continuous observation of position. Taking the initial state to be 
coherent, and denoting it \a), the initial position wave- function is given by 



2s- 



2\ V4 



e -s 2 x 2 +2sxa-±(\a\ 2 +a 2 ) ^ (3.49) 



where s 2 = mw/(2fi). Using the results in appendix |b| we find that the coefficient of x 2 at a later 
time t is given by 



S ' 2 = s 2 



1-21 
3-21 



1 - 21 

1 + 2 

1 + 2 



(3.50) 



where 

-1/2 
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and we have defined the parameters 



-A/7- 1 ' r=zM- ( 3 - 52 ) 



After some algebra this may be written as 



,., 2 - izta,nh(zujt) — 1 



s' z = s z iz / / , (3.53) 

tanh(zwt) — iz 

in agreement with that derived by Belavkin and Staszewski. The conditional variance for x at time 
t is given by 

<*(*)w = ^Ty (3-54) 



As t tends to infinity, Eq.( 3.53| ) gives a steady state value for the conditional variance, which is 



-l 



^ = 4linM = (^ S Vv / V^TT+l) ■ (3-55) 

The parameter r is a dimensionless quantity which gives essentially the ratio between the frequency 
of the harmonic oscillator, and the rate of the position measurement. We may view the dynamics 
of the position variance as being the result of two competing effects. One is the action of the mea- 
surement which is continuously narrowing the distribution in position, and consequently widening 
the distribution in momentum. The other is the action of the harmonic motion, which rotates the 
state in phase space, so converting the widened momentum distribution into position. Depending on 
the relative strengths of these two processes, determined by the dimensionless constant r, a steady 
state is reached in which they balance. If the rate of the measurement is very fast compared to 
the frequency of the oscillation (corresponding to r <C 1), then the localisation in position is much 
greater than it would be for an unmonitored oscillator, and in that case we succeed effectively in 
tracking the position of the particle. However, if the frequency of oscillation is much greater than 
the rate of localisation due to the measurement, then the steady state position variance remains 
essentially that of the unmonitored oscillator. 

3.2.4 An Open Question 

The results obtained in the previous section are curious in the light of the solution obtained by various 
authors for the non-linear equation describing a continuous position measurement for the case of an 
initial Gaussian state |5^, |55| . It is clear from this solution that the probability density for the final 
states is Gaussian. It is surprising therefore that the evolution operator for the equivalent linear 
equation describing this situation contains random variables which are not Gaussian. These two 
facts are compatible because the non-Gaussian variable in the evolution operator given in 



affects only the normalisation of the final state, and not the state itself. Because of this it is, at 
least in theory, possible to eliminate the non-Gaussian variable from the solution. This may be 
seen as follows. Let us assume that we have an initial state and an evolution operator which 
is a function of the random variables X and Z (which may in general be vector valued). We let 
the random variable Z determine only the normalisation of the final state, so that the evolution 
operator may be written 

V(X,Z,t) = 0(X,t)f(Z,t), (3.56) 
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where O is an operator valued function, and / is simply a complex valued function. The unnormalised 
state at time t is then given by 

\m) w = 0(X,t)f{Z 1 t)\i,). (3.57) 

Clearly once we have normalised that state at time t, it is no-longer dependent upon Z. In particular 
the normalised state is given by 

m)w _ mmk . (3.5 8 ) 

The probability density for the final state is 

P(X, Z, t) = (i,(t)m)) w P w (X, Z, t), (3.59) 

in which P W (X, Z, t) is the probability density given by the Wiener measure for the variables X and 
Z. However, seeing as the normalised state depends only upon X, we require for all calculations 
only the marginal probability density for X. Denoting this marginal density also by P, we have 

P(X,t) = J P(X,Z,t) dZ. (3.60) 

In certain cases the probability measure for the normalised state may therefore be Gaussian, even 
though the measure for the unnormalised state is not. However, as P(X, Z, t) contains a factor of the 
norm of \ip(t)) w , the probability measure for the output process will, in general, only be Gaussian 
if the norm is Gaussian in X. Clearly the norm is Gaussian in X for initial Gaussian states in the 



case we investigate in section 3.2.3 



In the case considered in subsection 3.2.3 it is not clear how to evaluate the marginal probability 
density in which the non-Gaussian variable has been eliminated because we do not have an expression 
for the joint distribution of all the variables including Z. Whether it is possible to obtain an evolution 
operator in which the variable Z has been eliminated remains therefore an open question. It seems 
almost certain, however, that once the non-Gaussian variable has been eliminated the resulting 
distribution for the final state will be Gaussian in the remaining variables. 



3.3 Extension to the Poisson Process 

A linear stochastic equation describing the evolution of a system undergoing quantum jumps takes 
the form 

\i)(t + dt)) n = [1 + Adt + BdN]\$(t)) n (3.61) 

in which dN is the increment of a Poisson process, and A and B are arbitrary operators. The 
Poisson increment, dN, is cither equal to zero or one, with the result that (dN) 2 = dN. If during 
any infinitesimal time interval dN — then no jump occurs, whereas if dN = 1 a jump occurs at 
that instant. The rate of jumps is given by A, so that (dN) — Xdt. The average rate of jumps may 
be chosen so as to render the expressions involved in deriving the solution in their simplest form. 
We now note that 

1 + BdN = e Ms+i)<w (3.62) 
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where we have used (dN) 2 = dN . Using this we may write the linear differential equation as 

\4>(t + dt)) N = e Adt e^ B+1 ^ N \m) N (3.63) 
where we have used dtdN = 0. The state of the system at time t may now be written as 

M 

|#)) N = lim Y\ e AAt e l ^ B +^ AN ^(0)) (3.64) 

m— 1 

where we have cut the time into M slices so that At = t/M. In order to derive a closed form for 
the evolution operator we must swap the exponentials in A and B repeatedly so as to move all 
the exponentials containing AiV to one side. Once that is achieved the infinitesimals must then be 
summed. Then, in the limit as M — > oo these sums form integrals. The complexity of the result will 
depend upon the commutation relations between A and B, and this procedure will not result in a 
closed form solution for the evolution operator if the commutation relations are two complex. To 
illustrate this method we turn now to a concrete example. 



3.3.1 A Damped Cavity with Coherent Driving 

We consider first an undriven optical cavity mode which is allowed to decay through one of the end 
mirrors, and then extend this to include coherent driving. A photo-detector is placed in the beam 
which is output from the cavity. The photons are detected one-by-one as they leak out of the cavity, 
and a quantum jump in the evolution of the system corresponds to the detection of a photon. A 
linear stochastic equation describing this situation is 

\ip(t + dt)) N = [1 - ((mj+ 7/2)0+0 -7/2)dt + (a- l)dN] \4>(t))x, (3.65) 

in which u> is the cavity resonance frequency, 7 is the cavity decay rate, a is the cavity annihilation 
operator, and we have chosen A = 7 to select the simplest expression. The operators A and B are 
therefore 

A = -(^ + 7/2)0+0-7/2, (3.66) 
B = o-l. (3.67) 



To swap the exponentials in Eq.(3.64) we now require a relation of the form 



ae ea ' a = e Ea ' a f(a). (3.68) 



Multiplying both sides on the left by e ea+a , we see that / is given by 

f(a) = e- £afa ae £ata . (3.69) 

This is simply the time evolution of the operator a in the situation in which the Hamiltonian is 
proportional to a^a. This is easily calculated with standard quantum optical techniques, and we 
obtain f(a) = ae 6 . We see that the first requirement for our procedure to produce an evolution 
operator in closed form is that we can find the evolution of the annihilation operator when A is 
taken as the Hamiltonian. The second condition is that operators that result from performing the 
swapping (that is, f(a)) allow the Poisson increments to be summed together. 
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Swapping all the terms in Eq.(3.64), and taking the limit as M — > oo, we obtain the evolution 
operator: 

|-0(£)> - e [-(-+T/2) a t a+T/2]* e -^+n'/2)^( i ) a JV(t)|^ (0)) _ (3 7Q) 

Here N(t) is simply the Poisson process, being the sum of the Poisson increments, and therefore the 
total number of photon detections up until time t. The second random variable which appears in 
the evolution operator, Y(t), is given by 

N(t) 

Y(t) = I t'dN{t') = V t k , (3.71) 



pt "W 

/ t'dN{t') = Y^t k , 
J " k=i 



in which the tk are the times at which the photo-detections have occurred up until time t. The 
probability distribution for N(t) is simply the Poissonian 

P(N,t)=e- xt ^. (3.72) 

The probability distribution for Y(t) may be calculated as follows. First note that the probability 
for each of the N photon-detection times, given that N detections have occurred, is uniformly 
distributed in the interval [0, t]. The probability density for the sum of the emission times, is the 
convolution of the probability densities for each emission time. The probability density for Y(t) 
given that there have been N emissions within time t, which is denoted by P(Y\N,t), is therefore 
the convolution of N probability densities of the form 

P(t') = ( T ° t * < * • (3-73) 
y ' \ otherwise v ; 

The joint probability density for N and Y is then P(N, Y, t) = P(N, t)P(Y\N, t). 

Examining the evolution operator, we see that the first term gives the evolution of the system 
resulting from the free Hamiltonian, including a smooth decay. The second term, containing Y, 
accounts for the alteration to this evolution from the loss of photons in the quantum jumps. The 
final term, being a product of annihilation operators, describes the effect of the jumps. 

The evolution operator may now be used to obtain the final state of the system for any initial 
state, given a particular trajectory. In particular, as the operators which appear in this evolution 
operator act on coherent states in a simple way, the final state of the system may be obtained most 
readily by writing an initial state in the coherent state basis. Averages over all trajectories may 
then be calculated by using the distributions for N and Y , along with the norm of the final state, 
which is also a function of N and Y. The most complex aspect of the result is the distribution for 
Y. However, due to the fact that Y determines only the norm of the final state, we may eliminate 
this variable from the solution. We now perform this elimination. 

We note that once we have normalised the final state, and ignored any overall phase factor, the 
final state will not depend upon Y(t). We do not require, therefore, the full probability density for 
the final states which is a joint probability in N and Y , but only the marginal probability obtained 
by integrating over Y. The probability density for the final states is given by 



p(N,Y,t) = m)m)) N p(N,Y,t). 



(3.74) 
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We may write (V>(i)|^0)>N = h{N,t)e-^ Y ^. The marginal probability may now be calculated as 
follows: 



P(N,t) = h(N,t)P(N,t) I P{Y\N,t)e-~> Y(t) dY 

n N 



h(N,t)P(N,t) 
h(N,t)P(N,t) 



P(t')e 



(1 - e-*) 
7* 



-7* 



A' 



(3.75) 
(3.76) 
(3.77) 



where in the second line we have used the fact that the tk are all independent. Now we have the 
probability density for the final states in terms of the norm calculated excluding Y. We may now 
write the state at time t using the simplified evolution operator which excludes Y, as 



|^(t))„ = e [-(^+7/2)« t «+7/2]* a JV(t)| V)(0))) 

where the probability for the final state may be written as 



p(N,t) = m)$(t)) n p(N,t) 



It 



N 



(3.78) 



(3.79) 



This evolution operator was first obtained by Garraway and Knight |34[ , using a somewhat different 
approach. They applied this to an initial 'Schrodinger Cat' state, which is a superposition of two 
coherent states \a) and |/3). As they did not give the probability density for the final states, we will 
use this as our example. Setting [3 = —a for simplicity so that the magnitude of the initial coherent 
states are the same, the initial state is 



|V(0)> 



V2(l + e- 2 M 2 ) 



(3.80) 



Applying the evolution operator in Eq.(3.78), and subsequently normalising, we obtain the state at 
time t to be 

U e -(-+7/2)A + (_l)JV(i)| _ ae -(iw+7/2)t\ 

Iw))n = — 



^2(1 + (-l^We^M^*) 
The probability density for the final states, given in Eq.(3.79), is found to be 



P(N,t) 



-|a| 2 (l-e-^) 



(M 2 (i 



-~/t\\N 



-|a| 2 (l+ e --')\ (_| a |2 (1 _ e - 7t)) 



f 



o-2\a\ 2 



AT! 



-2|a| 2 



Nl 



(3.81) 



(3.82) 



which is the sum of two Poissonians. We see from Eq.( [3.81 ) that at all times during the evolution 
the cavity mode remains in a 'Cat' state, although the magnitude of the 'Cat' decays smoothly away 
over time. In addition, the state toggles between an even and odd 'Cat' at each photon detection. 

Now that we have elucidated the procedure, let us extend the model to include coherent driving 
of the cavity mode. The operators appearing in the linear stochastic equation are now 



A = -(iw + 7/2)a t a+ {Ea) - E* a) -7/2, 
B = a-1, 



(3.83) 
(3.84) 
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where E is the amplitude of the driving field in units of s . We obtain the evolution operator for 
this equation as 

N(t) 

I#)>n = e At Y[ ((a - a s )e^/ 2t * + « S )|V(0)> , (3.85) 

k=l 

where a s = (2E/j). Note that this time the spontaneous emission times cannot be eliminated 
as in the previous case. Also, in this case the evolution operator is not strictly in a closed form. 
Nevertheless, it is sufficiently simple to be useful. If we take a coherent state, as the initial state, 
then noting that all the parts of the evolution operator which depend on the stochastic variables 
determine only the normalisation, we see that the evolution is completely deterministic, and we 
obtain for the normalised state at time i, 

|V(t)) N = I (a - a s )e-« t/2 + a s ) . (3.86) 

The state of the cavity therefore changes smoothly from its initial value to the steady state amplitude, 
a s . The jumps in this case have no effect upon the evolution. If the initial state was not coherent, 
however, the jumps would play a role in the evolution. 

3.3.2 A Damped Cavity with a Kerr Non-linearity 

For our final example we consider a damped cavity with a Kerr nonlinearity p5| . The operators 
appearing in the stochastic equation are now 

A = -{iw + 7/2)a f a - i|a t2 a 2 - 7 /2, (3.87) 

B = o-l. (3.88) 

In this case the evolution operator becomes 

|^(t)) N = e At e- (M+lxa+a+7/2)r(t) a Arw |V'(0)). (3.89) 

In fact, we note in passing that a double integral of the Poisson process appears in the derivation 
of the evolution operator for this case, but that it may be immediately discarded as it determines 
only an overall phase factor. Let us now examine the evolution of an initial coherent state, \a), in 
this system. To proceed we first note that the action of exp[irj(a^ a) 2 } on a coherent state results in 
superpositions of coherent states [ |66| . The number of coherent states in the superposition depends 
on the value of r\. For r\ — n/n the superposition contains n distinct coherent states. If we choose 
the evolution time so that %i = 7r we will obtain a superposition of two coherent states. Applying 
the evolution operator for this time, we find that the normalised state is 

m))n = (j^j lae-^W-^) + | _ ae -^mt-ixr(t) } _ (3 . 90) 

in which 

J\f(i) = 2Vl + exp(-2H 2 e-r f ) (3.91) 

is the normalisation. We see that the final state is a 'Schrodinger Cat', irrespective of the number of 
jumps. However, the orientation of the 'Cat' in phase space depends on Y(t), the sum of the times 
at which the jumps occurred. If we record the photo-detections, we can therefore keep track of the 
location of the 'Cat' in phase space. However, if we do not detect the emitted photons, the result 
of averaging over the emission times will phase average the final state. Consequently, in the limit of 
many photon emissions the final state will tend towards a mixture of number states. 
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3.4 Possibilities for Future Work 

So far we have examined the use of evolution operators for linear quantum trajectories in calcu- 
lating final states, their distributions, and the averages of conditional variances over the possible 
trajectories. More generally, however, quantum trajectories are useful whenever we are interested in 
properties of the state of the system during a continuous measurement. One example in particular is 
the question of the purity of the state of the system, which, to my knowledge, is only just beginning 
to be explored |67[ |. This is interesting because it tells us how much information the continuous 
measurement is providing us about the state of the system, rather than just certain system vari- 
ables. If the state of the system is initially impure, then we may view this as characterising a lack of 
knowledge about the initial state. As the measurement proceeds we would expect the state to purify, 
reflecting the fact that we are obtaining information. This has applications, among other things, to 
the preparing of pure states from thermal mixtures. Bose et al. [jss|| have presented a scheme for 
preparing a superposition state of a macroscopic object. This requires, however, that the object be 
initially prepared in a pure state. As all macroscopic objects exist in thermal states a process of 
purification would be required, and continuous measurement would provide a possible method for 
achieving this. 

3.5 Conclusion 

We have presented in this chapter a method for deriving explicit evolution operators for quantum 
trajectories when they are cast in a linear formulation. We have shown a number of things re- 
garding these solutions. First, for linear stochastic equations in which the operator multiplying 
the time infinitesimal and the operator multiplying the Wiener infinitesimal commute with their 
mutual commutator, it is possible to obtain an explicit evolution operator in which the classical 
random variables are Gaussian. In general this provides an explicit solution for equations of this 
form. Second, we have examined the situation in which the corresponding master equation contains 
measurement operators which are linear in position and momentum, and for which the Hamiltonian 
is quadratic in these operators. Note that in this case the quantum Langevin equations are linear, 
and can therefore be solved. The corresponding statement we are able to make about the Wiener 
process linear quantum trajectories is that in this case explicit evolution operators may be found. 
However, we have also found the surprising fact that the classical random variable appearing in 
the evolution operators are not all Gaussian. We have indicated that it may well be possible to 
eliminate the non-Gaussian variable, and obtain a completely analytic solution to equations of this 
form, although we have not been able to do this. We have also shown that it is possible to extend 
the technique to provide a unified approach to deriving evolution operators for both the Wiener 
process and the Poisson process. 
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Chapter 4 

Quantum Noise in a Cavity with a 
Moving Mirror 



In this chapter we consider an experimental realisation of the position measurement of a moving 
cavity mirror. We perform this analysis for two reasons. The first is to perform a quantum mechanical 
treatment of phase modulation detection, which is often used in practice, and which we show has a 
different shot noise level than homodyne detection. The second is to include the effects of various 
experimental sources of noise, such as intracavity losses and classical laser noise, so that we may 
examine how the quantum back-action noise appears among these sources. 



4.1 Introduction 

Interferometers provide a very sensitive method for detecting small changes in the position of a 
mirror. This has been analysed extensively in the context of gravitational wave detection Jo9| , [7C| , 
fnl |7^| and atomic force microscopes E], [73| . A key limit to the sensitivity of such position detectors 
comes from the Heisenberg uncertainty principle. The reduction in the uncertainty of the position 
resulting from the measurement is accompanied by an increase in the uncertainty in momentum. 
This uncertainty is then fed back into the position by the dynamics of the object being measured. 
This is called the quantum back action of the measurement, and the limit to sensitivity so imposed 
is referred to as the standard quantum limit. 

In real devices which have been constructed so far, the quantum back action noise in the mea- 
surement record is usually small compared to that arising from classical sources of noise. However, 
as the sensitivity of such devices increases it is expected that we will eventually obtain displacement 
sensors that are quantum limited. The quantum back-action noise has not yet been seen experi- 
mentally for macroscopic devices, and is therefore a topic of current interest. Once the standard 
quantum limit has been achieved, this will not be the end of the story, however. Various authors 
have shown that it is possible to use contractive states Q, or squeezed light fn|, to reduce the 
quantum back action and therefore increase the sensitivity of the measurement even further. 

The interferometer we consider here for measuring position consists essentially of a cavity where 
one of the mirrors is free to move. This system is also of interest from the point of view of cavity 
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QED. Usually cavity QED experiments require optical cavities where the atomic excitations and 
photons in the optical modes become entangled. The dynamics follows from the interplay between 
these quantum variables. However, a challenging realm for cavity QED experiments involves instead 
an empty cavity (that is, a cavity containing no atoms or optical media) where the photons in the 
cavity mode interact with the motion of one of the cavity mirrors. In this scheme, the position 
of at least one mirror in the optical resonator is a dynamic variable. The coupling between the 
photons and the mirror position is simply the radiation pressure that stems from the momentum 
transfer of 2hk per one reflected photon with the wavcnumbcr k. It has been been shown that this 
system may be used to generate sub-Poissonian light in the output from the cavity |7(], |77|, [T8|. 
The moving mirror alters the photon statistics by changing the optical path length in a way that 
is proportional to the instantaneous photon number inside the cavity. This system may also be 
used to create highly non-classical states of the cavity field, such as Schrodinger cats (6^, (79), and 
might even be used to create cat states of the mirror j68|. In addition, it has been shown that 
such a configuration may be used to perform QND measurements of the light field ]7£| |8(], |8l| , and 
to detect the decoherence of the mirror, a topic of fundamental interest in quantum measurement 
theory ||8^| . Due to recent technological developments in optomechanics, this area is now becoming 
experimentally accessible. Dorsel et al. have realised optical bistability with this system |f33| , and 
other experiments, particularly to probe the quantum noise, are now in progress. 

In order to create displacements that are large enough to be observed, one is tempted to use a 
mirror having a well-defined mechanical resonance with a very high quality factor Q. Thus, even 
when excited with weak white noise driven radiation pressure, the mirror can be displaced by a 
detectable amount at the mechanical resonance frequency v. For such a mirror to behave fully 
quantum mechanically one needs to operate at very low temperatures since the thermal energy kT 
very easily exceeds fiv. For example, a v/2ir — 100 kHz resonance is already significantly excited 
thermally at 5 (iK. However, it is not necessary to reach the fully quantum domain to observe 
the quantum back action. By simultaneously combining high optical quality factor (ie. by using 
a high-finesse cavity) and a specially designed low mass mirror with very high mechanical quality 
factor one can, at typical cryogenic temperatures, create conditions where the radiation pressure 
fluctuations (which are the source of the quantum mechanical back-action referred to earlier) exceed 
the effects caused by thermal noise. In this chapter we discuss considerations for detecting this 
quantum back-action noise. 

There are already a number of publications dealing with quantum noise in optical position 
measurements. Our main purpose here is to extend this literature in two ways which are important 
when considering the detection of the quantum noise. The first is the inclusion of the effects of 
experimental sources of noise, such as the classical laser noise and the noise from intracavity losses. 
The second is to perform a quantum treatment of phase-modulation detection, so that the results 
may be compared with those for homodyne detection. While this method of phase detection is often 
used in practice, as far as we are aware it has not previously been given a quantum mechanical 
treatment. In addition to these main objectives, we also show that the standard Brownian motion 
master equation is not adequate to describe the thermal damping of the mirror, but that the corrected 
Brownian motion master equation derived by Diosi B4| rectifies this problem. 



In section 4.2 we describe the configuration of the system. In section 4.3 we perform a quantum 
mechanical analysis of phase modulation detection. In section 4.4 we solve the linearised equations of 
motion for the cavity/mirror system, using a non-standard Brownian motion master equation which 
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is of the Lindblad form 84 . In section LI we use this solution to obtain the noise power spectral 
density (which we refer to simply as the spectrum) for a measurement of the phase quadrature using 
phase modulation detection. In the first part of this section we discuss each of the contributions 
and their respective forms. Next we compare the spectrum to that which results if the standard 
(non-Lindblad) Brownian motion master equation is used to describe the thermal damping of the 
mirror, and also to that which would have been obtained using homodyne detection rather than 
phase modulation detection. Finally we show how the error in a measurement of the position of 
the mirror may be obtained easily from the spectrum. We evaluate explicitly the contribution to 
this error from various noise sources, and plot these as a function of the laser power. Section 4.6 
concludes. 



4.2 The System 

The system under consideration consists of a coherently driven optical cavity with a moving mirror 
which will be treated as a quantum mechanical harmonic oscillator. The light driving the cavity 
reflects off the moving mirror and therefore fluctuations in the position of the mirror register as 
fluctuations in the light output from the cavity. In the limit in which the cavity damping rate is 
much larger than the rate of the dynamics of the mirror (characterised by the frequency of oscillation 
v and the thermal damping rate 7) the phase fluctuations of the output light are highly correlated 
with the fluctuations of the position of the mirror and constitute a continuous position measurement 
of the mirror |4l[] . 

An experimental realisation will therefore involve a continuous phase-quadrature measurement 
of the light output from the cavity to determine the output spectrum of the phase-quadrature fluc- 
tuations. The nature of the detection scheme used to measure the phase quadrature is of interest to 
us, as we shall see that it will effect the relationship of the shot noise to the other noise sources in 
the measured signal. Quantum theoretical treatments usually assume the use of homodyne detec- 
tion |Il], 76 , f77|, [78| . However this is often not used in practice p3, |6| . Many current experiments 



use instead phase modulation detection ]87| , which was developed by Bjorklund |88| in 1979. Before 
we treat the dynamics of the cavity field/oscillating mirror system, to determine the effect of various 
noise sources, we will spend some time in the next section performing a quantum mechanical treat- 
ment of phase modulation detection. We will focus on this scheme throughout our treatment, and 
compare the results with those for homodyne detection. A diagram of the experimental arrangement 



complete enough for the theoretical analysis is given in Fig. 4.1. We note that in practice a feedback 
scheme is used to lock the laser to the cavity so as to stabilise the laser frequency. For an analysis 
of the method and an expression for the resulting classical phase noise the reader is referred to 
references |8f| and Q. We do not need to treat this feedback explicitly, however. Its effect may be 
taken into account by setting the value of the classical laser phase noise in our analysis to the level 
it provides. 



4.3 Phase Modulation Detection 



The laser which drives the cavity is isolated from the cavity output, and the entirety of this output 
falls upon a photo-detector. In order that the photo-detection signal contain information regarding 
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Figure 4.1: The light output from the laser is modulated at frequency A, and from there drives 
the cavity. The front mirror of the cavity is fixed, while the back mirror is a mechanical harmonic 
oscillator. The diagonal line is a shorthand representation for the arrangement which isolates the 
laser from the cavity output, and the light that is reflected back off the front mirror. All this light 
falls upon a photo-detector, and the photo-detection signal is demodulated (to pick out the phase 
quadrature signal) before going to a spectrum analyser. 

the phase quadrature, the laser field is modulated at a frequency A, which is chosen to be much 
greater than the natural frequency of the harmonic mirror. The sidebands that result from this 
modulation are far enough off-resonance with the cavity mode that they do not enter the cavity 
and are simply reflected from the front mirror. From there they fall upon the photo-detector. The 
result of this is that the output phase quadrature signal appears in the photodetection signal as a 
modulation of the amplitude of a 'carrier' at frequency A. This is then demodulated (by multiplying 
by a sine wave at the modulation frequency and time averaging) to pick out the phase quadrature 
signal, and from there the spectrum may be calculated. 
First let us denote the laser output field by 



In this expression /3 is the average coherent amplitude of the field, which we choose to be real. The 
deviations from this average are given by Sx(t), being the (classical) amplitude noise, and Sy(t), 
being the (classical) phase noise. The quantum noise, which may be interpreted as arising from the 
vacuum quantum field, is captured by the correlation function of the field operator <$a; n (i). Here the 
subscript refers to the field's relation to the cavity, and not the laser. The correlation functions of 
the various noise sources are 



We have left the classical noise sources arbitrary, as this allows them to be tailored to describe the 
output from any real laser source at a later time. The units of (3 are such that, if a photo-detector 
was placed in the beam, 1 would give the average rate of photo-detection. The average values of 
the three noise sources are zero, as are all the cross correlation terms. 
After modulation the field becomes 



P + 5a in {t) + Sx(t) +i5y(t). 



(4.1) 



(5a in (t)Sal(t + r)) = 6(t), (Sal{t)5a in {t + r)> = 0, 
(Sx(t)Sx(t + r)) = G x (r), (6y(t)6y(t + r)) = G y (r). 



09 + Sx(t) + iSy{t)){\ + ee 



-iAt 



ee tAt )+Sa in (t). 



(4.2) 
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Note that modulation affects the coherent amplitude of the input field, and has therefore no effect on 
the quantum contribution. In order to achieve phase detection it is necessary to modulate the phase 
quadrature, which is the reason we choose the modulation to be exp(— iAt) — exp(iAi) = 2zsin(Ai). 
Using now the input-output relations of Collett and Gardiner [^9| , the field output from the cavity 
is 

flout (t) = -(j3 + Sx(t)+t6y(t))(l + se- iAt ~ee iAt )-6a in {t) + ^f(6a(t)+a), (4.3) 

in which a(t) = Sa(t) + a is the operator describing the cavity mode and 7 is the decay constant 
of the cavity. We are interested in the steady state behaviour, and we choose a to be the average 
steady state field strength in the cavity. In addition, in order to solve the equations of motion for the 
cavity we will linearise the system about the steady state, which requires that (5cv (t)Sa(t)) <C |a| 2 . 
The operator describing the photo-current from the photo-detector is 



I(t) = a out (i) t a out (t) 

= a 2 + a(5X out + 25x) + 2e(3 sin(At) ( SY C 



'out 



2^/ady 





+ (2 £ /3) 2 sin 2 (At)(l + ^), (4.4) 

in which 

5X out (t) = ^6X{t)-SX in {t) , SY out (t) = ^6Y(t)-5Y in (t), 

and 

6X(t) = Sa(t) + <Ja + (t), SY(t) = -i(5a(t) - 6at(t)), 
5X in (t) = 6a in (t)+5al(t), SY in (t) = -i(5a in (t) - 5al(t)). (4.5) 

We have also set a = {—(3 + y/ja), and assumed this to be real. Now, we see that the intensity 
signal, I(t), contains the phase quadrature multiplied by sin(Ai). That is, the phase quadrature 
signal appears as a modulation of a carrier at frequency A. We therefore choose A to be much larger 
than the bandwidth of the phase quadrature fluctuations (this bandwidth being determined by the 
dynamics of the cavity-mirror system), and the classical laser noise, so that the phase quadrature 
signal appears alone, centred at the frequency A, with the sole exception of the quantum noise, 
which, being white, also has components at the frequency A. To extract the phase quadrature 
signal, and also, unavoidably, the quantum noise, we perform a demodulation procedure. This 
involves multiplication by a sine wave at frequency A, and subsequent averaging over a time, T. 
The multiplication by the sine wave shifts the signal both up and down by the frequency A. As a 
result of this there is now a component of the phase quadrature signal back at zero frequency. We 
then extract just the signals around zero frequency by applying a low pass filter, being accomplished 
by averaging. The averaging time T must be much larger than 1/A so as not to pass any of the 
extraneous signals now appearing centred at the frequencies ±A and ±2A, but much smaller that 
the time-scale of the system dynamics so to pass all of the signal around zero frequency without 
distortion. The final signal is therefore given by 

R(t) = ^ / sin(Ai)J(t + r) dr. (4.6) 
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We now must evaluate this to obtain R(t) explicitly in terms of the phase quadrature. Writing out 
the integral, and dropping everything which averages to zero (that is, which is not passed by the 
low-pass filtering) we obtain 



R(t) 



where 



<&(*) 
92 (t) 



SY out {t + T)dT 



5y{t + t) dr 



le -^{t+T) sx in {t + T) dr 



Re 



e- l2A(t+T) SY in (t + t) dr 



+ 91 (*) + *(*), (4-7) 

(4.8) 
(4.9) 



Note that we choose T to be much smaller than the time-scale upon which SY and 8y change, so that 
the integration is essentially equivalent to multiplication by T, an effect which is cancelled by the 
division by T. However, we should also note that i5y ou t contains <5Yi n , so that in replacing the first 
term in R(t) by e(3 SY(t) ou t we must remember that this only contains the frequency components of 
5Y[ n in a bandwidth of 1/T around zero frequency. The result of this is that SY(t) ont is uncorrelated 
with qi and q 2l being the quantum noise in the bandwidth 1/T around the frequencies A and 2 A 
respectively. We need to know the correlation functions of these noise sources, and whether or not 
they are correlated with any of the other terms in R(t). It is clear that qi and q2 are not correlated 
over separation times greater than 2T. Using Eq. (|4.9|) to evaluate the correlation function of q 2 , for 
example, we have 







otherwise. 



On the time-scale of the fluctuations of SY we can approximate this as a delta function, so that qi 
and q2 (and also 5Yi n ) are still effectively white noise sources. We may therefore write 



R(t) = e/3SY(t) out + qi (t)+q 2 (t)-2e^a5y(t), 

and the correlation functions of q\ and q 2 are 

(gi(t)gi(f + r)) = (1/2)09 - ^jaf 5{t), 
(q 2 (t)q 2 (t + T)) = (l/2)(e/3) 2 5(t). 



(4.11) 



(4.12) 



The signal therefore contains the phase quadrature of the output field, 5Y out (t), plus three noise 
terms. While the last term, being the input classical phase noise, is correlated with 5Y out (t), q\ and 
q 2 are not. Taking the Fourier transform of the signal to give 



'2tt 



R(b)e 



dt. 



(4.13) 



we may write 

R(u) = e/3 6Y out ((j) + qi {uj) + q 2 (w) - 2 £v ^aSy(uj). (4.14) 
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This is the Fourier transform of the signal in the case of phase modulation detection. If we were 
to use ideal homodyne detection this would be instead pq] 

R h (co) = K0(8Y out (u) - 2Sy(u)), (4.15) 

where $ is the amplitude of the local oscillator and k is the reflectivity of the beam splitter used 
in the homodyne scheme. Thus, in the case of phase modulation detection, there are two white 
noise sources which do not appear in homodyne detection. They stem from the fact that the phase 
quadrature detection method is demodulating to obtain a signal at a carrier frequency. Because 
the quantum noise is broad band (in particular, it is broad compared to the carrier frequency) the 
demodulation picks up the quantum noise at A and 2A. There is also a term from the classical 
phase noise in the sidebands. For real homodyne detection there will also be an extra contribution 
from the noise on the local oscillator. 

Once we have solved the equations of motion for the system operators, we will obtain 5Y(u>) in 
terms of the input noise sources. We can then readily calculate (R(uj)R(uj')}, which appears in the 
form 

(R(lu)R(lu')) = S{lu)5(lo + uf). (4.16) 

The delta function in w and u>' is a result of the stationarity of R(t), and S(u>) is the power spectral 
density, which we will refer to from now on simply as the spectrum. This is useful because, when 
divided by 27r, it gives the average of the square of the signal per unit frequency (The square of 
the signal is universally referred to as the power, hence the name power spectral density). Since 
the noise has zero mean, the square average is the variance, and thus the spectrum provides us 
with information regarding the error in the signal due to the noise. The spectrum is also a Fourier 
transform of the autocorrelation function |l2|. The specific relation, using the definitions we have 
introduced above, is 

/oo 
(R{0)R(T))e- lWT dT, (4.17) 
-oo 

and as the autocorrelation function has units of s~ 2 , the spectrum has units of s _1 . To deter- 
mine the spectrum experimentally R(t) is measured for a time long compared to the width of the 
auto-correlation function, and the Fourier transform is taken of the result. Taking the square mod- 
ulus of this Fourier transform, and dividing by the duration of the measurement obtains a good 
approximation to the theoretical spectrum. We proceed now to calculate this spectrum. 



4.4 Solving the System Dynamics 

Excluding coupling to reservoirs, the Hamiltonian for the combined system of the cavity mode and 
the mirror is J8l| 

p 2 1 

H = Tuujqo) a H 1 — mv 2 q 2 ~ Tiga) aq 

2m 2 

+h {% [E + ^5x{t) + iy/j5y{t)} a f + H.c.} . (4.18) 

In this equation u>o is the frequency of the cavity mode, q and p are the position and momentum 
operators for the mirror respectively, m and v are the mass and frequency of the mirror, g = loq/L 
is the coupling constant between the cavity mode and the mirror (where L is the cavity length), 
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and a is the annihilation operator for the mode. The classical driving of the cavity by the coherent 
input field is given by E which has dimensions of s -1 , and is related to the input laser power P by 
E = y Pj/(hLu ) — The classical laser noise appears as noise on this driving term. 

The moving mirror is a macroscopic object at temperature T, and as such is subject to thermal 
noise. While it is still common to use the standard Brownian motion master equation (SBMME) 
p8| , p2[ to model such noise, as it works well in many situations, it turns out that it is not adequate 
for our purposes. This is because it generates a clearly non-sensical term in the spectrum, which will 



be shown in subsection 4.5.1. As far as we know this is the first time that it has been demonstrated 
to fail in the steady state. Discussions regarding the SBMME, and non-Lindblad master equations 
may found in references (9^, |9l| . We will return to this point once we have calculated the spectrum. 
We use instead the corrected Brownian motion master equation (CBMME) derived by Diosi [pi[ , to 
describe the thermal damping of the mirror, as this corrects the problems of the SBMME. Using this, 
and the standard master equation for the cavity losses (both internal and external), the quantum 
Langevin equations of motion for the system are given by 

a = ~[a,B\- (^^a + ^fSa ta (t)+^/jib in (t), (4.19) 

q = -Uq,H]+h(T/QmkT)^r)(t), (4.20) 
Ti 

V = ~\p,H]-Tp + (2mTkT)iat), (4.21) 

in which the correlation functions for the Brownian noise sources are 

mm) = 5(t-t% (v(t)v(t'))=S(t-t'), _ (4.22) 
mv(t')) = -i(V3/2)S(t-t'), (r,(t)Z(t'))=i(V3/2)6{t-t'). (4.23) 

In these equations all internal cavity losses including absorption, scattering and transmission through 
the movable mirror are included via the input operator bi n (t) and the decay constant fi. The effect of 
mechanical damping and thermal fluctuations of the mirror are given by the classical noise sources 
£(i) and rj (t) and the mechanical damping constant T. 

We note here that if we were to use the standard Brownian motion master equation [B2l |48|, 



Eqs.(4.20) and (4.21) would instead be given by 



q = ~\q,H\, (4.24) 
p = -^[p,H}-Tp+(2mTkT)^C(t). (4.25) 

where (CMC(O) = S(t — t'). These Langevin equations do not preserve the commutation relations of 
the quantum mechanical operators, and as a result it is clear that the description cannot be entirely 
correct. 



Calculating the commutators in Eqs.(4.19) to (4.21) 



a = E- (j^^j a + igaq + ^5a iu {t)+^b iu {t) + ^5x{t)+i^5y{t), (4.26) 

q = — + h<r/6mkT)iT](t), (4.27) 
m 

p = -mv 2 q + hga^a-Tp+ (2mTkT)?£(t), (4.28) 
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in which g is the strength of the coupling between the cavity field and the mirror. Introducing a 
cavity detuning A (that is, setting the cavity resonance frequency in the absence of any cavity field 
to lo c = uo Q + A), and solving these equations for the steady state average values we obtain 



(a) ss = 

(<Z}ss = 
(P)ss = 



2E 



0, 



(4.29) 

(4.30) 
(4.31) 



where we have set the detuning to A = g(q) ss to bring the cavity on resonance with the driving field 
in the steady state. Linearising the quantum Langevin equations about the steady state values, and 
writing the result in terms of the field quadratures, we obtain the following linear equations 
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In this set of equations we have scaled the position and momentum variables using 



SQ 
SP 
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(q - (g) ss ), 



(p - (p) ss ), 



I 

(4.32) 
(4.33) 



and we have defined \ = g{2h/mu) 1 ^ 2 , which has units of s 1 . The quadratures for the input noise 
due to intracavity losses are given by 



5X bt [ 
6Y Wl 



(4.34) 
(4.35) 



Without loss of generality we have chosen the input field amplitude to be real (Im[/3] =0), so that 
the input phase quadrature is given by Y[ n . We now solve these in the frequency domain in order 
to obtain the spectrum directly from the solution. To switch to the frequency domain we Fourier 
transform all operators and noise sources. In particular we have, for example 



1 r 00 

Sa (u) = -j= J Sa(t)t 

Sa\uj) = -= / 5(J(t)e iajt dt = [5a(-u)]t. 

Rearranging the transformed equations, the solution is given by 

(SX(lu),SY(lu),SQ(lu),SP(lo)) t = M(lu)u(lu), 



(4.36) 
(4.37) 

(4.38) 



where n(w) is the vector of transformed noise sources. If we write the matrix elements of M(u>) as 

Hj 



Mij(ui) — niij{uj) / D{lo), then D{lo) = ((7 + ji)/2 — iui) 2 (v 2 — lo 2 — iTto), and the non-zero ma are 
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given by 

mn = m 2 2 = ((7 + A*)/2 - iu>) (v 2 - lu 2 - iTuj) , m 2 i = X 2 a 2 ^, 

= X"(r - i<*>) ((7 + M)/ 2 - iu) , rn 2 A = x av ((7 + m)/ 2 - iu) = m 31 , 

TO33 = (r - iw) ((7 + fj)/2 - iu f , m 34 = v ((7 + /i)/2 - iu) 2 = -TO43, 

TO41 = — zxaw ((7 + yu)/2 — iw) , m4 4 = — iu: ((7 + fi)/2 — iu) 2 . 

We have now solved the equations of motion for the system in frequency space. The spectra of the 
system variables may now be calculated in terms of the input noise sources. Using the input-output 
relations, which give the output-field in terms of the system variables and the input noise sources, 
the spectra of the output field, and hence of the measured signal, may be obtained. Note that 
quantum mechanics plays no role in the solution of the motion of the system. The linear equations 
of motion may as well be equations for classical variables. The only part that quantum mechanics 
plays in determining the spectra of the system variables is that some of the input noise sources are 
quantum mechanical. That is, their correlation functions are determined by quantum mechanics. In 
fact, if all the noise sources had purely classical correlation functions, then the SBMME Langevin 
equations would not lead to any problems, as they are perfectly correct as equations of motion for 
a classical system. 



4.5 The Power Spectral Density 

To calculate the spectrum of the signal, we require the correlation functions of the input noise 
sources. To reiterate, these are 

(6X in (w)6X in (uj')) = (SY in (u;)SY in (w')) = S(lu + <J), 
(SX in (w)6Y in (uj')) = -(SY in (w)6X in (uj'))=iS(w + w'), 

and similarly for SX^i n (u>) and £Yb,i n (w). The correlation functions for the classical laser noise, and 
thermal noise sources are 

(Sy(w)6y(uj')) = G y (u)S(w + J), 

mm) = (v(tHt')) = s(t- t % 

mm) = -mw))=i(V3/2)s(t-n (4.39) 



After some calculation we obtain the spectrum of the signal for phase modulation detection as 
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where 

|£M| 2 = [( 7 + + c 2 ] 2 - u?f + r 2 ^ 2 ] , (4.41) 

and T s is a dimensionless scaled temperature given by T s — (fcs / '{Kv))T '. This phase-fluctuation 
spectrum may be thought of as arising in the following way. The mechanical harmonic oscillator, 
which is the moving mirror, is driven by various noise sources, both quantum mechanical and classical 
in origin, and the resulting position fluctuations of the mirror are seen as fluctuations in the phase 
of the light output from the cavity. 

Let us examine the origin of the various terms in the spectrum in turn. The first two terms, which 
appear in the first set of square brackets, are independent of the frequency, and are the contribution 
from the (quantum mechanical) shot noise of the light. 

The next three terms, which multiply the second set of square brackets, are the back-action of the 
light on the position of the mirror, noise from internal cavity losses, and the classical amplitude noise 
on the laser, respectively. Note that the only distinction between the back-action and the internal 
losses is that the former is proportional to the loss rate due to the front mirror, and the latter is 
proportional to the internal loss rate. It is easily seen that these noise sources should have the same 
effect upon the position of the mirror: the back-action is due to the random way in which photons 
bounce off the mirror, where as the internal losses are due to the similarly random way in which 
photons are absorbed by the mirror, (or anything else in the cavity). The amplitude fluctuations of 
the laser also affect the mirror in the same manner, all though since these fluctuations are not white 
noise (as is the case with the quantum noise which comes from the photon 'collisions'), the response 
function of the mirror is multiplied by the spectrum of the amplitude fluctuations. 

The term which appears in the third set of square brackets, is due to the classical phase fluctu- 
ations of the laser. Clearly this has quite a different form from that due to the quantum noise and 
the classical amplitude fluctuations. In particular, it is not dependent upon the coupling constant, 
<7, because it is derived more or less directly from the input phase noise. Conversely, the noise that 
derives from the amplitude fluctuations has its origin from the fact that the amplitude fluctuations 
first drive the mirror, and it is the resulting position fluctuations which cause the phase fluctuations 
in the output. The classical phase noise term has a contribution from the phase noise in the side- 
bands and from the phase noise on the light which has passed through the cavity. The final two 
terms, which multiply the forth set of square brackets, are due to the thermal fluctuations of the 
mirror. Note that these terms are only valid in the region in which k^T 3> hu. 

Finally we note that we do not see squeezing in the spectrum of phase quadrature fluctuations. 
This is because squeezing is produced when the cavity detuning is chosen so that the steady state 
detuning is non-zero j77j . We have chosen to set the steady state detuning to zero in this treatment 
as we are not concerned here with reducing the quantum noise. 

In what follows we examine various aspects of the spectrum which are of particular interest. 
Before discussing considerations for detecting the back-action noise, we compare the spectrum with 
that which would have been obtained using the SBMME, and for that which would result from the 
use of homodyne detection. We then write the spectrum at resonance as a function of the laser 
power, and plot this for current experimental parameters. So far we have been considering the noise 
power spectrum, and have made no particular reference to the limit this implies for a measurement 
of the position of the mirror. In the final part of this section we show how the spectrum tells us the 
limit to the accuracy of position measurement in the presence of the noise sources. 
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4.5.1 Comparison with Standard Brownian Motion 



As we have mentioned above, to obtain the spectrum we have used the corrected Brownian motion 
master equation p4| . This is essential because the standard Brownian motion master equation 
produces in this case the following term in the spectrum: 



2uj"fTx 2 'a 2 v 



\D(w)\* 



(4.42) 



This is clearly incorrect since it is asymmetric in u>. It follows readily from the stationarity of the 
output field, and the fact that the output field commutes with itself at different times, that the 
spectrum must be symmetric inw. In particular, the stationarity of the output field means that the 
correlation function of the signal only depends on the time difference, so that 



(R(t)R(t + t)) =G(t). 



(4.43) 



As the output field commutes with itself at different times, R commutes with itself at different times, 
and we have 

G{-t) = (R(t)R(t - t)) = (R(t - r)R(t)) = G(t). (4.44) 

The correlation function is therefore symmetric in r. As the spectrum is the Fourier transform of 
the correlation function, it follows from the properties of the Fourier transform that the spectrum 
is must be symmetric in ui. 

Diosi's corrected Brownian motion master equation removes the term asymmetric in to by adding 



a noise source to the position (see Eq.(4.21)) which is correlated with the noise source for the 
momentum. In doing so it produces an additional term in the spectrum proportional to 1/T. For 
temperatures in which his is much smaller than ksT (T s 3> 1) this term is very small and can 
be neglected. However, the question of probing this term experimentally is a very interesting one, 
because it would allow the CBMME to be tested. 



4.5.2 Comparison with Homodyne Detection 

Let us now briefly compare the spectrum derived above for phase modulation detection to that 
which would be obtained with homodyne detection. Firstly, if homodyne detection had been used, 
the overall scaling of the spectrum would be different, as it would be proportional to the strength of 
the local oscillator. Thus the factor of l/(e/3) 2 would be replaced by 1/(/3k) 2 , in which j3 and k are 



as defined in Eq. ( 4.15 ). This overall factor aside, two terms in the spectrum would change. The 
shot noise component would be reduced to unity, and the classical phase noise contribution would 
become 



4G» 
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(4.45) 



There would also be a contribution from the noise on the local oscillator. 



4.5.3 The Error in a Measurement of Position 



So far we have been considering the noise spectrum of the phase quadrature, as this is what is 
actually measured. However, the reason for performing the phase measurement is that it constitutes 
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essentially a measurement of the position of the mirror. In this section we show how the error in a 
measurement of the position of the mirror may be obtained in a simple manner from the spectrum, 



Eq.(4.40), and give an example by calculating it explicitly for some of the terms. 

We can choose to measure the amplitude of position oscillations at any frequency, but for the 
purposes of discussion, a measurement of a constant displacement is the simplest. First we must see 
how the position of the mirror appears in the signal, which is the phase quadrature measurement 
(that is, convert from the units of the signal into units of the position fluctuations). This is easily 
done by calculating the contribution to the spectrum of the position fluctuations due to one of the 
noise sources (for the sake of dcfinitcness we will take the thermal noise), and comparing this to the 
equivalent term in the spectrum of the signal. This gives us the correct scaling. Performing this 
calculation, we find that the spectrum of position fluctuations of the mirror due to thermal noise is 



given by the thermal term in the spectrum (Eq.( [4.40 )), multiplied by the factor 

From this we see that the scaling factor is frequency dependent. This means, that the spectrum of 
the position fluctuations is somewhat different from the spectrum of the resulting phase quadrature 
fluctuations. For the measurement of the phase to correspond to a true measurement of the position 
the two spectra should be the same. This is true to a good approximation when 7 is much larger 
than the range of u> over which the spectrum of position fluctuations is non-zero, and this is why 
the scheme can be said to constitute a measurement of position when 7»i/,r. 

In performing a measurement of a constant displacement of the mirror (achieved by some constant 
external force), the signal (after scaling appropriately so that it corresponds to position rather than 
photocurrent) is integrated over a time T m . The best estimate of the displacement is this integrated 
signal divided by the measurement time. The error, Ax, in the case that the measurement time is 
much greater than the correlation time of the noise, is given by 

/>oo 

Ax 2 -/ (R x (0)R x {T))dr/T m = S x {0)/T m . (4.47) 



In this equation R x and S x are the appropriately scaled signal and spectrum. To calculate the error 
in the measurement of a constant displacement, all we have to do, therefore, is to scale the spectrum 



using the expression Eq.(4.46), evaluate this at zero frequency, and divide by the measurement time. 
In general, the spectrum evaluated at a given frequency gives the error in a measurement of the 
amplitude of oscillations at that frequency. We calculate now the contribution to the error in a 
measurement at zero frequency and at the mirror resonance frequency, from the shot noise, thermal, 
and quantum back-action noise. In the following we write the expressions in terms of the parameters 
usually used by experimentalists. That is, the cavity finesse, T = nc/(2L"f) (assuming intracavity 
losses are small compared to 7), and the quality factor for the mirror oscillator, Q = v/(2T). The 
contribution due to the shot noise is the same at all frequencies and is 

A 2 7T 2 T3 11 (hc 2 \ 1 
The contribution from the quantum back action for a measurement of a constant displacement is 
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Error in a position measurement vs. Laser Power 
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Figure 4.2: The error in a measurement of the position of the mirror, for a measurement time of 
T m — Is. The dashed curve corresponds to a measurement of a constant displacement, and the dot- 
dash curve to a measurement at the mirror resonance frequency. The contribution of the quantum 
back-action noise to both curves is the sloping section to the right. The range of laser powers used 
in the experiment lies between the solid lines. Clearly the experimental parameters are in the regime 
in which the quantum back-action noise may be observed. 



and for an oscillation at frequency v is Ax| A (^) = 4Q 2 Ax| A (0). The contribution from the thermal 
noise is 

fc R T 



for a constant displacement, and is Ax^ h (V) = 4Q 2 Aa;^ h (0) for an oscillation at the mirror frequency. 
The contribution from the other noise sources may also be readily evaluated from the terms in the 
spectrum Eq.(4.4C). 

Let us examine the total error in a position measurement resulting from these three contribu- 
tions for state-of-the-art experimental parameters. Reasonable values for such parameters are as 
follows j87|. The cavity resonance frequency is ujq — 2tt x 2.82 x 10 15 Hz (assuming a Nd:YAG laser 
with a wavelength of 1064 nm), the cavity length is L = 1 cm, the mass of the oscillating mirror is 
m = 10 -5 kg, and the resonant frequency of the mirror is v = 2ir x 2 x 10 4 Hz. The quality factor 
of the mirror is two million, which gives r = 5 x 10~ 2 s" 1 . With these parameters for the cavity we 
have x — 1-49 x 10 -17 s _1 . The cavity damping rate through the front mirror is 7 = 10 6 s _1 , and we 
will take the damping rate through the oscillating mirror to be negligible in comparison. This gives 
a 2 — 2.13 x 10 9 , and we take e = 0.2. The cavity may be cooled to a temperature of T — 4.2 K, so 
that T s = k-QT/ifiv) = 4.37 x 10 6 , which is certainly in the high temperature regime (T s 1). 

In figure 4.2 we plot the position measurement error as a function of the laser power, both for the 
measurement of a constant displacement, and for a displacement at the mirror resonance frequency. 
We have plotted this for a measurement time of one second. Conversion to another measurement 
time merely involves dividing by the square root of the measurement time given in seconds. 

The uncertainty due to the shot noise falls off with laser power, while that due to thermal noise 
is independent of laser power, and that due to the quantum back-action increases with laser power. 
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These results are already well known. The thermal and shot noise contributions are much greater 
at the resonance frequency of the mirror, due to the high mechanical Q factor. Clearly the optimal 
regime for obtaining a precision measurement of position is off-resonance, where the back-action and 
thermal noise are suppressed. However, the optimal regime for detecting the quantum back-action 
noise is at resonance, as the absolute magnitude of this noise is largest in this case. Reasonable 
experimental values for laser power lie between the solid lines, where the increase in noise due to 
the back-action is visible. However, our analysis of the spectrum shows us that the full situation is 
more complicated. We have shown that the internal cavity losses and the classical laser noise have 
the same dependence on frequency as the quantum back-action, and the size of these contributions 
must therefore be determined accurately for a given experimental configuration so that they may be 
distinguished from the back-action. 

4.6 Conclusion 

We have examined the optomechanical system consisting of an optical cavity containing a moving 
mirror to see how the quantum mechanical back-action appears among the various sources of classical 
noise. We have shown a number of things regarding this question. First of all, the relationship of 
the shot noise to the noise resulting from the oscillating mirror, and hence the limit on a position 
measurement due to the shot noise, is dependent on the phase measurement scheme. In particular, 
the result for phase modulation detection, which is commonly used in experiments of this kind, is 
not the same as that for homodyne detection. We have found that while the signature of the classical 
phase noise is quite different for that of the quantum-back action, the noise due to intracavity losses 
and classical amplitude noise has a very similar signature to the back-action. However, as far as the 
parameters of the cavity and oscillating mirror are concerned, realisable experiments are beginning 
to fall in the region where the quantum back-action may be observed. 

In our treatment of the system we have shown that the standard quantum Brownian motion 
master equation produces a clearly spurious term in the steady state noise spectrum for the phase 
quadrature measurement. We have shown that the corrected Brownian motion master equation, 
derived by Diosi, corrects this error. However, it also produces a new term in the spectrum which 
is small in the high temperature limit. Testing for the existence of this term now poses a very 
interesting experimental question. 

In the introduction to this chapter we pointed out that many authors have shown that there are 
a large number of possible uses for optomechanical cavities, including various sensor applications, 
microscopy, gravity wave detection, the generation of non-classical light both internally and in the 
output beam, basic cavity QED experiments, QND measurements, and for probing the decoherence 
of macroscopic objects. Testing the quantum theory of thermal noise, by probing for the existence of 
the Diosi term in the thermal spectrum, provides yet another reason for the importance of developing 
these techniques. The detection of the quantum back-action noise in such cavities will represent 
another milestone in this development. 
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Chapter 5 



Determining the State of a Single 
Cavity Mode 

In the previous chapters we have been concerned with quantum systems subjected to the continuous 
measurement of a physical observable. In this chapter we consider an alternative problem in quantum 
measurement theory: the determination of the entire quantum state of a system, rather than simply 
a particular observable. We consider a single cavity mode and show that the quantum state of such 
a system my be obtained from the photon statistics measured both before and after the interaction 
with one or more two-level atoms. 

5.1 Introduction 

Determining the quantum state of a system is essentially just the process of obtaining the probability 
distribution of a random variable by taking many samples of that variable. The application of 
this process to quantum states is, however, complicated by the fact that the state determines a 
set of probability amplitudes for the values of any given physical observable, rather than just a 
probability distribution ]9^ |. It is then further complicated if one wants to consider mixed states 
as well as pure states. To determine the quantum state of a single degree of freedom, given that 
it is pure, we must obtain the probability distribution for at least two non-commuting physical 
observables. To see why this is so we need merely note that when expressed in the basis of any 
physical observable, the wavefunction consists of a complex number for each value of the observable. 
The probability distribution for that observable gives only the moduli of this set of complex numbers, 
and consequently we must obtain the distributions for at least two physical observables to obtain 
the quantum state. 

If we wish to determine a mixed state of a single degree of freedom then we must obtain the 
distributions for many observables. To gain an insight into why this is so we need to consider how 
many independent real numbers it takes to specify mixed states. Clearly, to be able to invert the 
measured data to obtain the quantum state, there must be at least as many independent real numbers 
in the data as are required to specify the state. To specify a mixed state requires a density matrix 
or Wigner function. While the wavefunction describing a pure state is merely a function of one 
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variable, these are functions of two variables. As the probability distribution for a single observable 
is one-dimensional, it only gives us information equivalent to a single slice of the two-dimensional 
Wigner function. As a consequence we must measure the distributions for many different observables 
to be able to reconstruct the Wigner function, and hence a mixed state, with any accuracy. 

Schemes for determining the quantum state of a system do not necessarily correspond directly to 
measuring the probability distributions of a set of observables. More generally, we can speak of a set 
of parameters for which data values are taken. Clearly, to obtain enough independent data values 
to reconstruct the Winger function, we must obtain data which is parametrised by two continuous 
co-ordinates. To reconstruct a pure state, on the other hand, only one of the parameters need be 
continuous (or countably infinite). The other need merely take a small number of discrete values. 

There are now a variety of methods which have been suggested for the reconstruction of quan- 
tum states. For example the quantum state of a single optical field mode may be determined by 
tomographic reconstruction from the distributions of all the rotated quadratures obtained using ho- 
modyne detection [ 53 . A number of methods for the measurement of a field mode which involve 
the use of beam splitters have also been suggested J95 |, and tomographic reconstruction has been 
suggested and realised for other quantum systems p6[~97|, [|. In addition, schemes for measuring a 
single cavity mode by using two-level atoms as probes have been suggested by Wilkens and Meystre 
(nonlinear atomic homodyne detection) ps| , and by Bardroff et al. (quantum state endoscopy) [ p9| . 
Both these methods involve passing a resonant two-level atom through the cavity containing the 
state to be measured, and measuring the final atomic excitation probability as a function of the 
transit time of the atom through the cavity. In the former scheme the atom interacts with a refer- 
ence field in addition to the cavity mode, and the final atomic excitation probability is also measured 
as a function of the phase of this reference field. In the latter scheme, while there is no reference 
field, the measurements are performed for four different initial atomic states. We note also that 
recently a scheme for measuring the state of a single cavity mode which involves the use of atomic 
magnetic sublevels has been suggested by Walser et al. [ 100 1 . In this scheme an atom is sent through 
the cavity in order to map the cavity state onto the magnetic sublevels, and then the final atomic 
state is determined from repeated measurements using a tomographic technique. 

Both Homodyne tomography and nonlinear atomic homodyne detection may be used to measure 
mixtures as in these schemes measurements are taken for two continuous parameters. On the other 
hand, in quantum state endoscopy there is one continuous parameter, and one finite discrete param- 
eter (the discrete parameter being the initial atomic state, of which just four are used), so that only 
pure states may be reconstructed with this scheme. However, quantum state endoscopy requires a 
somewhat simpler inversion procedure than the other two methods. 

We present here various schemes for measuring the state of a cavity mode in which either prac- 
tically no inversion, or very simple inversion of the measurement data is required to reconstruct 
the state; the expansion coefficients of the state in the Fock basis are in an essentially one-to-one 
correspondence with the measurement data. We consider only pure states, however, and while these 
schemes are applicable to most pure quantum states currently of interest (number states, coherent 
states, squeezed states, Schrodinger cats) the schemes we present in detail here are not applicable 
to all pure cavity states. The first two, which involve allowing the cavity mode to interact with 
one two-level atom, enable the measurement of all pure cavity states which do not contain 'zeros' 
in the photon probability distribution. By a 'zero', we mean that one of the photon probabilities, 
say P n , is zero, while there are both probabilities P m >n, and P m <n which are non-zero (ie. the 
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photon distribution is split into two or more parts). An example of a state in this category is the 
squeezed vacuum. The two-atom scheme, which involves allowing the cavity to interact with two 
two-level atoms, allows the reconstruction of all pure cavity states which have no contiguous (adja- 
cent) zeros in the photon probability distribution. That is, the photon distribution may have any 
number of parts, so long as each of the non-zero photon probabilities is separated by no more than 
one vanishing photon probability. In general if n two-level atoms are sent through the cavity, then 
pure quantum states with n — 1 contiguous zeros in their photon probability distribution may be 
determined. However, as the number of atoms increases, the data inversion that is required increases 
in complexity. 

The measurement procedure for all these schemes consists of three steps. First the photon 
statistics of the cavity state are determined by measurements of photon number on a set of identically 
prepared cavities. This might be achieved by photon-counting the number of photons which leak 
out of the cavity with a photo-detector (applicable in the optical regime) , or by a variety of other 



methods [101]. Next, a second set of cavity modes, all prepared in the same initial state as before, 
are allowed to interact for a specified time with either one or two two-level atoms. The photon 
statistics of this set is then determined. In the case of a single atom, we may choose whether or 
not we measure the final atomic states for each cavity. Finally, the second step is repeated, but the 
atom(s) are prepared in a different initial state, giving a third set of photon statistics. The three 
sets of photon statistics may then be used to calculate the phase factors of the coefficients of the 
initial state in the number basis. Along with the initial photon statistics, which give the amplitudes 
of the number state expansion coefficients, this determines the complete quantum state. Clearly 
this procedure is similar to that required for both atomic homodyne detection and quantum state 
endoscopy. The main difference is that while in the latter two schemes the final atomic excitation 
probabilities are measured for many different interaction times, in this scheme just one interaction 
time is chosen, and the final photon probabilities are measured. Thus in this case the cavity state 
is determined from a measurement of discrete variables only. 

In the following section we introduce the general principle on which our measurement schemes are 
based, and define some notation. In Section p. 3 we describe the measurement schemes which involve 



an interaction with one atom, and in Section p.4| we describe the schemes involving an interaction 
with two atoms. Section 5.5 concludes. 



5.2 Preliminary Comments 

In the Fock (number) basis, the state \ip) of a single cavity mode of the electromagnetic field may 
be written as 

oo 

|V)=$> n e if, "|n), (5.1) 

where \n) is the quantum state in which the cavity contains n photons, so that P n — r\ is the 
probability that the single mode will be measured to contain n photons. To obtain the photon 
statistics (the photon probabilities P n ), we may measure the number of photons in the cavity for a 
large number of cavities prepared in the state \ip). If we were now to form a new state which was 
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the superposition of \ip) with a known state, for example that given by 

oo 

71=0 

then the photon statistics of the resulting superposition state contains information about the phase 
factors 9 n . If we create the superposition state 

\Xi) = -j=m + \4>)), (5-3) 
then the photon probabilities for this state, 

Qi.r, = Kx|n>| 2 , (5.4) 



are given by 

so that the cosine of each phase is 



Qi,„ = P„(l - cos#„), (5.5) 



cos#„ = <9i,„/P n - 1. (5.6) 
Alternatively, if we form the superposition 

|X2> = ^=(|V>+#», (5.7) 
then the photon probabilities for this state, denoted Q,2,n, are given by 

Q 2 ,n =P„(l + sin0 n ), (5.8) 
so that the sine of each phase factor is given by 

Bin0 n = l-Q 2 , n /P n . (5.9) 

As knowledge of the sine and cosine of the angle completely determines the angle itself, the three 
probability distributions, P n ,Qi, n an d Qi,m which may be measured by photo-counting, allow the 
complete quantum state to be calculated in a simple way. 

Unfortunately, it is not possible to create the superposition of an unknown quantum state with 
a known quantum state; both quantum states must be known. We can see this from the following 
argument. It is clear that given two states \<j>) and \x) we can find a unitary evolution operator U 
with the property 

m = /o^op m m (I0) + (5 - 10) 

y/2 + 2Re[{<j>\x)\ 

To construct such a unitary operator we need merely to choose two bases \4>i) and \ipi) such that 
\4>) = |0i) and (|0) + \il>)) /y/2 + 2Rep|x>] = IV'i)- The desired operator is then U = ^ 
However, the unitary operator which is required to form the superposition of \<f>) and |%) is dependent 
in general upon the initial state \<j)). This follows immediately from the linearity of U. If we assume 
that the same unitary operator will perform the operation 

m = /,Jf/,| M (IV,) + (5 ' n) 
v /2 + 2Re[(V'|x)] 
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where ^ |</>), then using the linearity of U we find that in general 

U{a\<j>) + /%)) ? l(a|0) + PW) + \ X )), (5.12) 

where N is the normalisation constant. It follows that in general for different initial states a different 
evolution operator is required to create the superposition of the initial state with a known state |x)- 
Before we can find the correct evolution operator to create the superposition of two states, we must 
know both of those states. We conclude, therefore, that we cannot create the superposition of two 
unknown quantum states, or even the superposition of a known quantum state with an unknown 



one, just as we cannot clone one unknown quantum state [102]. 

Nevertheless, with a simple atom-field cavity QED interaction we may perform measurements 
along the lines described above. Allowing a cavity mode to interact with a two-level atom in the 



Jaynes-Cummings manner [ 103 forms superpositions of adjacent number states. In this way we may 
use a procedure similar to that described above to measure the phase differences between adjacent 
number states in the superposition which makes up the quantum state to be measured. However, 
this only allows us to measure the complete quantum state if there are no zeros in the photon 
distribution. If a certain number state |n) does not contribute to the state to be measured (ie. 
f n = 0) when there are contributing number states on either side of it, then as we cannot determine 
the phase difference between those number states and the number state that does appear in the 
superposition, we cannot obtain the phase difference between the two states on either side. To do 
this we require to superpose number states which are separated by more than one photon. This 
is why sending two atoms through the cavity allows us to measure quantum states with zeros in 
the photon distribution, so long as none of the zeros are adjacent. Coherent states [ p4| and number 
states are among those which may be measured with a one-atom interaction, and squeezed states 
and 'Schrodinger cats' (being superpositions of two coherent states) are among those which may be 
measured with a two-atom interaction. 

For the calculations in the following sections, it will be convenient to define compact notation 
for certain functions. In the expressions that follow n represents a non-negative integer, and t and 
n„ represent real numbers. We will use the following definitions: 

T f = cos(ft„i) 

T+ = cos(O n+ it), l£ = cos(tt n+2 t) 

Tl = cos(n n _i*) , d = cos(ft„_ 2 i). (5.13) 

We will use the same notation for sine functions, where T will be replaced with T. When the 
cosine or sine functions are squared, the T and T will have tildes added, thus 

f* = cos 2 (fW) , T* = sm 2 (Q n t). (5.14) 

When it is convenient we will use a number subscript for T and T, and in this case the definitions 

r^ = cos(fi n t), r t n = sm(n n t), (5.15) 



will hold. This notation will greatly reduce the space required for the expressions derived in the 
following sections, without introducing any ambiguity. 
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5.3 Measuring quantum states with a one- atom interaction 



As we stated in section p.2| , the state of the single cavity mode to be measured may be written in 
terms of the number states as 

oo 

IV) = $> n e ie "|n). (5.16) 

Once the photon distribution, P„, has been determined, we create a new state in which the photon 
distribution depends on the phases 9 n . This may be achieved by passing a two-level atom through the 
cavity so that it interacts with the cavity mode via the familiar Jaynes-Cummings Hamiltonian |103|, 
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H = oj<7 z + Ma* a + k{<7 + a + a a'). (5-17) 

Denoting the ground and excited atomic states by \g) and |e) respectively, the operators in this 
expression are 

o z = |(|e)<e| - \g)(g\) ,with a+ = \e)(g\ , a~ = (<r+)t, (5.18) 

the latter two being the respective raising and lowering operators for the atom. The annihilation 
operator for the cavity mode is denoted by a, oj is the frequency of the cavity mode (which is resonant 
with the atomic transition), k is the strength of the (dipole) coupling between the atom and the 
field, and we have set Ti — 1 for convenience. As the atom passes through the cavity, it interacts 
with the cavity mode, and depending on its initial state it will either 'begin' to emit a photon into 
the mode, or 'begin' to absorb a photon from the mode. If the atom is sent into the cavity in a 
superposition of ground and excited states, then the state of the mode will be transformed essentially 
into a superposition of the initial mode plus a 'partial' extra photon, or the initial mode minus a 
'partial' photon. (The first case if the atom exits in the ground state, and the second case if the 
atom exits in the excited state.) Hence, allowing the mode to interact with a resonant two-level 
atom will create a superposition of the mode state with a shifted version of itself. We show now how 
the phases 8 n may be determined from photon counts performed on the mode after the interaction. 
We will see that forming a superposition state by using two copies of the original state is not ideal 
for the purposes of quantum state measurement, in that it only allows a certain set of quantum 
states to be measured. In the following analysis we will work in the interaction picture, where the 
interaction Hamiltonian is given by Hj = a + a~ a)). 

Before we give results for the general case in which the cavity may contain a superposition of 
arbitrarily many number states, we examine the simpler situation in which the intra-cavity state is 
in a superposition of zero and one photons. This example contains the basic idea behind our scheme 
which we simply generalise afterwards. Let the initial cavity state be written as 

|x(0))=r |0)+r ie <ei |l)- (5.19) 

We take the initial atomic state to be l/^/2(\g) + e l ^\e)). After the atom has interacted with the 
cavity for a time t the joint 'cavity-atom' state is (after omitting the 1/V2 prefactor) 

r |0, 5 )+ [roe^-zne^T*] |0, e) + [ne" 1 ^ - ir e^T<] \l,g) 
+ ne^ 1+ ^ri|l,e) -ir x e i{ - e ^T\\2,g). (5.20) 
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If we denote the probability that the cavity will be found to contain zero photons after the interaction 
by Qo/2, then from Eq.(5.20) we obtain the following equality 



Qo = P (l + f )+P l r o + 2^/P^T t o r o sm(6 1 -<i>). (5.21) 

Clearly this contains all the necessary information regarding the phase 6\, since we can obtain Pq, Pi 
and Qo by photon counting techniques. This now completely determines the cavity state. 

We turn now to the general case. Taking the initial state of the cavity mode to be |x(0)) = 
and allowing the atom to enter the cavity in the state l/y/2(\g) + e l<6 |e)), the state of the atom plus 
cavity system after an interaction time t is 



CO 

V2\ X (t)) = r o e ie °\0,g) + £ [r„e^"+«lt - ir^e^T^ 



\n,e) 



n=0 



+ ^[ r n+^ ien+ ^ t n-irne i ^ + ^T t n ]\n+l,g). (5.22) 

n=0 

In this expression the T's and T's are functions of f2„ as defined in the previous section, and in 
addition fl n — KyrT+T. 

Denoting the joint probability that after the interaction both the cavity is measured to contain 
n photons, and the atom is found to be in its excited state as Q„/2, and the corresponding joint 
probability for the atom to be found in its ground state by Q^/2, we have, for n > 1 



Qn = P n f * + P„+if * + 2V-P„Pn+ir t T* sin(A0 n+ i - (f>), (5.23) 

Qi = P n fl + P n -irl - 2 v /p i P„-iT t r t sin(A0„ - 0), (5.24) 

where A9 n = 9 n+ \ — 9 n . For n = we have 

Qo = Pofo + Pif + 2^F 1 r a r o sm(A9 1 -q>) 1 (5.25) 

Q 9 o = Po- (5.26) 



We denote the total probability for the cavity to contain n photons after the interaction by Q n /2, 
which is clearly given by Q n /2 = Ql/2 + Q e n /2. 

As the overall phase factor for the state is unimportant, we may set the phase 9q equal to zero, 
and the set of phase differences A9 n gives us complete information about the phase factors. Solving 
the above equations for the sine of the phase differences we see that these are easily calculated from 
the photon probabilities P n , Q e n and Q 9 n by 

P f * + P i T* — O 9 
sin(A#„-0) = " n - x ^ — (5.27) 

-(A^-0) = VPA _ iM ■ ( 5 - 28 ) 

for all n. Therefore, choosing the initial atomic state so that = we obtain the sine of every phase 
difference, and choosing the initial atomic state so that 4> = tt/2 wc obtain the respective cosines. 
Naturally once we have the sines and cosines the phase differences A9 n are well determined. In each 
case we obtain a measurement of all the phase differences from both Q e n and Q 9 n . 
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We may also obtain the phase differences without measuring the final atomic state. In this case 
performing photon counting on the final cavity state gives us the overall photon probabilities Q n . 
The expressions for these are 



Qo = Po(l + r? ) ) + P 1 T o + 2 v AfV^r*T*sin(A0„-0), (5.29) 

for n = 0, and 

Qn = P„(f ' + f *) + P n +if 1 + P^rXf + 2^/P n P n+1 T t T t sin(A0„ +1 - 0) 

-2 v /P„P„_iX. t r_ t sin(A#„ - dp), (5.30) 

for n > 1. We may calculate A6*i from Qo using 

Bin(A0!-0 = Q °~ P " (1+r ° ) ~ Pin . (5.31) 
The other phase differences may then be obtained by recursion, using 







/pi + :' 





sm(A0 n+1 -</.) = 2 v /P tI P„ +1 r*T ~ W 11 ' 1 S " 

It is clear from the above formulae that only phase differences between adjacent states may be 
measured in this way. If there is a hole in the photon probability distribution, so that P n = for 
some n, then we cannot measure the phase difference between the coefficient for that number state 
and the two on either side. As a consequence the above schemes will allow us to measure only pure 
states which do not contain zeros in the photon distribution. 



5.4 Measuring quantum states with a two-atom interaction 

To measure all the phase factors for a state in which there are zeros in the photon number distribu- 
tion, we may pass two atoms through the cavity consecutively. This creates a state which contains 
superpositions of number states separated by both 1 and 2 photons. In this case phase differences 
between number states separated by 2 photons may be measured by photon counting, so that zeros 
in the photon distribution are not an obstacle to the measurement so long as no two are adjacent. 

To analyse this measurement scheme, we may write the initial state of the system, which consists 
of three uncorrelated subsystems, the cavity and two atoms, as 

oo 

l*> = Ey eie "l")®(l-9) + e# l e ))®(lff) + l e ))- ( 5 - 32 ) 

ra=0 

We allow the first atom to interact with the cavity for a time t, and then, at some time after 
that, we allow the second atom to interact with the cavity for a time s. In the interaction picture, 
the expression for both the state and the joint probabilities, Q n 9 ,Q n e ,Q n 9 ,Q n e , following these 
interactions are fairly complicated, and are given in Appendix From these probabilities the 
phase differences may be calculated. If we are measuring the phase factors for a state for which all 
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the probabilities for odd numbers of photons are zero (for example a squeezed vacuum), then the 
expressions for the phase differences simplify considerably. Setting s = t, we have in this case 

P f*f * + P oT*T* — O" 
cos(A0 n -0) = " - - + n - 2 - , (5.33) 

2^P^T 2 WTl 

cos(A0 n -0) = t ^"~ 1 "' 2 ' " , (5.34) 
2VP„Pn- 2 r t T*T* 

cos(A#„-0) = , ~— ' - - , (5.35) 

2 v /P„P„- 2 riX*T* 

p p£pi i p 'Tt'Vt n ee 

cos(A6>„ - 0) = ; , (5.36) 

2 v /PnPn-2CT*Ti 

where A#„ = 0„ — 9 n - 2 , and rt is even. If we do not measure the final atomic state, then we may 
use the recursion relationship 

cos(A0 2 -0) = Ml + 2f ° + ^ + f ° f I ' C ° SW) + P2nf{ , (5.37) 
cos ( A^ +2 -0) = + + m + P n+2 (f^) - Q n 

— , cos(A#„ — <p . 



Similar expressions hold when all the probabilities of even numbers of photons are zero, and these 
are also easily derived by setting the P n for n even to zero in the expressions for the Q„ x given in 
Appendix ^|. 

For the general case, in which there are a number of non-adjacent zeros in the photon probability 
distribution, we suggest two procedures which may be used to calculate the phase factors from the 
joint probabilities. Inspection of the expressions for Q 99 , Qf^Li, Qn-i an d Q^1 2 shows that setting 
4> to zero, each is a linear combination of sin(A#„), sin(A0 n _i) and cos(A#„). Any three of the 
expressions for the Q's may therefore be linearly inverted to obtain these three trigonometric quan- 
tities. Setting <f) = 7r/2, the expressions for the four Q's are then a linear combination of sin(A6>„), 
sin(A# n _i), cos(A(9„), cos(A6>„-i) and sin(A#„). The first two of these have been determined by 
the previous step, so that there are again three unknowns to determine by linear inversion of the 
expressions for three of the Q's. Those two steps determine the sines and cosines of the phase 
differences, and consequently the phase differences themselves. If, for example, P m is zero, then 
the calculation of 9 m +i — @m-l is simpler since the terms containing A9 n and A9 n+ i vanish in the 
expressions for the Q's. 

Alternatively we may use a recursive procedure. Inspection of the expressions for Qq E and Qq 9 
shows that these may be used to determine the sine and cosine of A9i. Once we know that, there 
are only two unknowns rather than three to determine from Q 9 2 s , Qf e , Ql 9 and Qq 6 , making the 
inversion simpler. From these four Q's one of the unknowns we obtain is A9 2 , so that there are 
again only two unknowns that we need to determine from the next set of Q's, which is Qf 9 , Q 2 e , 
Q 2 9 and Q\ e , and so on up the recursion chain. 
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5.5 Conclusion 

We have shown that measuring the photon statistics of a single mode cavity field both before and 
after an interaction with a single two-level atom may be used to reconstruct in a very simple manner 
pure cavity states which contain no zeros in the photon probability distribution. This uses the fact 
that when a two- level atom initially prepared in a superposition of ground and excited states interacts 
with a cavity it creates superpositions of adjacent number states, so that the photon distribution 
after the interaction contains information about the phase factors of the coefficients in the expansion 
of the initial state in the Fock basis. We have also shown that this may be extended to measure 
pure states which contain no adjacent zeros in the photon distribution when two two- level atoms 
are allowed to interact with the cavity consecutively. In general it is clear from our analysis that 
sending n consecutive atoms through the cavity will form superpositions of Fock states separated 
by n photons, and therefor should allow the reconstruction of pure quantum states which contain n 
contiguous zeros, although the linear inversion that is required becomes increasingly complicated. 



Chapter 6 

Conclusion 



In this thesis we have presented work on a number of topics in the field of quantum measurement and 
quantum noise. In Chapter ^ we were concerned with the linear formulation of quantum trajectories, 
which describe the evolution of the state of a system subjected to continuous measurement processes. 
In contrast to the Schrodinger equation, which is sufficient to describe the evolution of a closed sys- 
tem, and for which evolution operators may be obtained simply by exponentiating the Hamiltonian, 
obtaining evolution operators for the linear stochastic equations describing quantum trajectories is 
not straightforward. We presented a method for doing this, and showed that it would render explicit 
evolution operators for a variety of stochastic equations describing various measurement processes. 

The quantum trajectory formulation, because it describes the evolution of a quantum system 
for every realisation of the measurement process, may be used to calculate aspects of the dynamics 
of the measured system not possible with the use of a master equation. One such aspect, which is 
yet to be investigated in detail, is the purification of an initially mixed state during a continuous 
measurement. This purification reflects the fact that information about the state of the system is 
being obtained continually during the measurement, and will undoubtedly be the subject of future 
work. 

Our ability to manipulate and observe individual quantum systems is continually increasing, 
and continuous quantum measurement theory is therefore likely to become increasingly relevant 
experimentally. Experiments are now being developed to probe the quantum noise in the position 
measurement of macroscopic objects, such as the configuration we examined in Chapter ^, and 
position measurements of single atoms in a cavity are now beginning to be realised [105]. The 
continuous tracking and control of the position of a single quantum degree of freedom is unlikely to 
be far away. 

In Chapter |^ we considered the question of the determination of the state of a quantum system 
from a set of measurements. This subject has been a topic of intense activity in the past few years, 
and there are now many measurement schemes which have been suggested for this purpose. These 
are now of practical importance in determining the states of quantum systems under our control. 

It is not by chance that the systems we have been concerned with here are quantum optical. This 
field has proved a rich ground for both testing the efficacy of, and the development of new ideas in, 
quantum measurement theory. Probably the most famous example of the use of quantum optics for 
the former is the testing of Bell's inequalities E3] . Regarding the latter, both the quantum trajectory 
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and quantum noise formulations of open quantum systems were developed in the field of quantum 
optics (although the former was also developed in a fairly parallel fashion by pure measurement 
theorists [p7| ). In addition, the extension of continuous measurement theory to include continuous 
feedback was also developed in this field @. Quantum optics has also proved a fertile ground 
for the development of many kinds of measurement schemes, such as QND measurements, joint 
measurements on incompatible observables, and the reconstruction of quantum states, among others. 
Note that in using the term 'quantum optics' we mean to include also the study of non-relativistic 
atom-laser and atom-atom interactions as well, such as the cooling, trapping and manipulation of 
atoms and ions. These are now providing the most powerful means yet of manipulating quantum 
states f97f . 

Recently a new focus has arisen in quantum measurement theory, due to the emergence of 
quantum information theory. The laws of quantum information theory stem from the dynamics of 
quantum systems, and the rules of quantum measurement theory. As a consequence, the formalism 
of quantum measurement theory, and in particular the formulation in terms of Positive Operator- 
Valued Measures, forms the basis of quantum information theory. However, in this case the focus is 
on the information contained in quantum systems, and ways in which this may be manipulated. 

Clearly a two-state quantum system may contain one classical bit of information. However, 
whether this is all that such a system may contain is not such a trivial question. A two-state 
quantum system is actually described by two real parameters. Specifying these parameters with 
arbitrary accuracy requires an arbitrarily large amount of information, and from this point of view 
the amount of information which may be stored in a two-state quantum system is arbitrarily large. As 
a consequence it is the amount of information that may be extracted during a measurement process 
which gives the maximum storage capacity of a quantum system. Similar arguments apply to the 
amount of classical information which may be transmitted by using entangled quantum systems, 
and quantum measurement theory is now being employed to answer questions of this nature. 

We may also consider a new kind of information, quantum information, as distinct from the 
standard 'classical' kind defined in Shannon's information theory |10£]. That this is a useful concept 
is illustrated by the following two examples. The first of these is the discovery that computations 
may be performed on information contained in quantum two-state systems, or 'q-bits', that cannot 
be perform on the usual classical variety [ 109 . As a consequence, the transmission of quantum states, 
and therefore quantum information, takes on an importance of its own. A second example is the 
use of entanglement between two quantum systems for secret communication, a procedure referred 
to as quantum cryptography Jl8[] . When two classical systems are correlated, they are referred to 
as containing 'mutual information'. However, it is certainly not possible to use two classically 
correlated systems for secret communication. The 'mutual information' between two entangled 
quantum systems is therefore distinct from that between two classically correlated systems, and its 
properties must be studied accordingly. Now that entanglement is seen as a useful resource, quantum 
measurement theory is being used to investigate its manipulation. 

We see that quantum measurement is becoming increasingly relevant, not only to the experimen- 
tal investigation of individual quantum systems, but also to the exploration of new ways in which 
the nature of the quantum world can be put to practical use for the manipulation of information. As 
time goes by no doubt we will begin to see the new theoretical developments becoming a practical 
reality. 
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An Introduction to the 
Input-Output Relations 

Here we present an introduction to the Input-Output formulation of open quantum systems. A 
basic knowledge of Quantum Electrodynamics is assumed. The derivation of the one-dimensional 
Lagrangian, the modal decomposition of the reservoir light field, the nature of the input and output 
fields and their relationship to initial and final conditions are discussed. 

A.l Introduction 

Input-Output theory, or Input-Output formalism, is a particular formulation of the quantum treat- 
ment of open systems. Open systems are systems coupled in some way to an environment, or 
reservoir, which has many degrees of freedom. The result of the coupling is that the system may 
exchange energy with the reservoir, and the reservoir introduces noise into the system. 

Open systems see extensive use in Quantum Optics; an atom coupled to the quantum electro- 
magnetic field, which as a result undergoes spontaneous emission, is an example of an open system, 
and we have found that we may treat the decay of optical cavities in essentially the same manner. The 
Input-Output formulation is a one dimensional theory, and is particularly useful for treating optical 
cavities; the modes of the free quantum electric field outside the cavity, which form the reservoir, 
may be decomposed into travelling fields which are input to, and output from, the cavity. The Input- 
Output treatment of open optical cavities is equivalent to the Master equation treatment p3| , but 
where as the Master equation treats only the evolution of the system, the Input-Output formulation 
allows the calculation of the output, as it shows that this may be written in terms of the input and 
the system operators. In fact, if it is desired, the system evolution may be solved using the Master 
equation (rather than the equivalent quantum Langevin equations of the Input-Output formulation), 
and then the output may be calculated using the relations of Input-Output theory. 

Input-Output theory was first formulated in a paper by Collett and Gardiner |89tl , and their 
treatment appears in the texts by Gardiner [Q and by Walls and Milburn|Q. In treating a cavity 
that experiences damping through one end mirror, we are only interested in the free electromagnetic 
field in half of one dimensional space; from the lossy cavity mirror out to infinity. The approach 
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taken by Collet and Gardiner is to treat the field in only this half-space. This results in a modal 
decomposition of the free field in which each mode operator corresponds to a mode travelling in both 
directions. Other treatments, which instead use a unidirectional field may be found in the book by 
Carmichael |22|], and in Wiseman's PhD thesis 0|. Our treatment here is essentially an expansion 
upon that given by Gardiner. We will, however, follow the notation in Cohen- Tannoudji et al. |106[. 



In Section A.2 we introduce the Lagrangian and Hamiltonian that are used for the treatment 
of Input-Output theory. We show what steps are necessary to arrive at these from the standard 
Lagrangian of non-relativistic quantum electrodynamics. In Section |A.3| we show how the external 
light field may be decomposed into a field travelling towards the system, the 'input' field, and one 
travelling away from it, the 'output' field. By specifying the external light field at some initial time, 
the input field is determined for all later times. The 'output' field may be written as the sum of 
the reflected input field and a field radiated by the system. In the final section we examine the 
input-output relations for the specific example of a damped optical cavity. 



A.2 The Lagrangian and Hamiltonian 
A.2.1 The Standard Lagrangian 

The starting point for our discussion is the Lagrangian of Quantum Electrodynamics in the Coulomb 
gauge. We assume the reader is familiar with this description of the interaction of light and matter. 
For readers unfamiliar with the quantisation of the electro-magnetic field the excellent text by Cohen- 



Tannoudji, Dupont-Roc and Grynberg | 106 | is highly recommended. This text is also recommended if 
the reader is unfamiliar with the treatment of Maxwell's equations in Fourier space. The Lagrangian 
in question is 

L = \Y^ m ^l +^Coui+ [j-Ad 3 r+ fd 3 r£ (A.l) 

a J J 

C = £2 [A 2 + c 2 (V x A) 2 ] . (A.2) 

In these equations C is the Lagrangian density for the fields, and J j • Ad 3 r describes the interaction 
between the particles and the fields. The index a labels the particles, where r Q = (r°,ry,r°) are 
the particle positions and m a the particle masses. The vector potential is denoted by A, and j is 
the current. For a discrete number of particles this may be written as j(r) = J2 a lai^i^r — r a ), 
where q Q are the charges of the particles. The vector potential is transverse, so that in real space 
it satisfies the condition V • A = 0. This condition is much easier to understand in Fourier space, 
which will be introduced below. The Lagrangian determines the equations of motion of the positions 
of the particles and the values of the fields at each point in space via Lagrange's equations []: 

d (d_L_\ dL_ ^ dC 

dt^Aj - QM .^%(a jAi ) ' A - ,i 



d ( dL\ _ dL 
dt \dr a J dr a 



(A.4) 



1 We use the functional derivative here. If the reader is not familiar with the functional derivative, then it is only 
necessary to note that in this case it means differentiate what ever appears under the integral sign in the expression 
for the Lagrangian with respect to the variable in question. 
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The positions of the particles, and the values of the fields at each point in space are the variables 
governed by the Lagrangian. 

If we take the spatial Fourier transform of the fields we obtain new field variables which are 
linear combinations of the old variables. The spatial Fourier transform of a field A(r), along with 
its inverse, is given by 

A(k) = U=\ 3 J A(r)e- lkr d 3 r , A(r) = U=\ " J A(k) e i;kr d 3 fc. (A.5) 



We will refer to the fields which are functions of position as 'real space' fields, and the Fourier 
transformed fields, which are functions of the 'Fourier space' variable k, as k-space fields. We will 
also use an abuse of notation and drop the bar over the k-space field, writing it simply as A(k). Note 
that the inverse transform relation gives the real field as a sum over plane wave 'modes' e . If we 
were to evolve the field in time in the absence of particles, then the plane waves would propagate. 
The direction of propagation of these waves is given by the direction of k. 

The dynamic equations are simplified by the process of taking the Fourier transform. The 
equations which would include spatial derivatives in real space contain only time derivatives in 
Fourier space; the equations are now ordinary differential equations. The result is that the field 
at each point in Fourier space is a dynamical variable uncoupled to the field at any other point 
(actually in the presence of particles the field variables will be coupled to each other indirectly via 
the motion of the particles). 

We also note that transverse fields are easily characterised in k-space. A general vector field has 
three degrees of freedom at each point in k-space. A transverse field on the other hand satisfies 
k • A(k) = 0. That is, the field at each point k in Fourier space is transverse to k. It has therefore 
only two degrees of freedom at each point. 

To obtain the correct Lagrangian for the k-space fields, all we have to do is substitute in the 



Lagrangian Eq.(A.l) for the real space fields in terms of the k-space fields. Note that this is not 
the same as simply taking the Fourier transform of the Lagrangian. If the k-space fields were real, 
then Lagrange's equations (with the real-space fields replaced by the k-space fields) would generate 
the correct equations for these new k-space variables 0. Actually the situation is complicated by the 
fact that the new fields are complex, but a thorough derivation shows that the correct equations 



are generated by treating the field and its complex conjugate as independent variables [ 106 1 . The 
correct Lagrange equations for the fields are now 
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E a* 



d{3 3 A*) 



(A.7) 



2 This can be seen as follows from the fact that Lagrange's equations determine the correct fields to produce a 
stationary point of the action integral. If the action is at a stationary point for a particular choice of the real fields, 
then when we substitute for the real fields in terms of the corresponding k-space fields it will remain stationary for that 
particular choice of the k-space fields. The correct equations for the k-space variables are therefore those equations 
which produce a stationary point of the action integral with the Lagrangian written in terms of these fields. These 
are of course simply Lagrange's equations written using the k-space fields, which is what we wanted to show. 
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A. 2. 2 The Input-Output Lagrangian 
Dropping two Dimensions 

For the purposes of Input-Output theory, we are only interested in propagation in two directions, 
one given direction, and its reverse. We therefore wish to consider only the fields along a line in 
Fourier space which passes through the origin. (In fact, later we will reduce the description further 
to just half the real line.) To do this we choose the x-axis as our direction of propagation, and take 
a finite normalisation area in the x and y directions. Taking this normalisation area to be L 2 , where 
L is the normalisation length for both y and z, the Fourier decomposition is now partially discrete 
and partially continuous: 

A(r) = 1 V] f A nm (k x ) e «*r.+i(2im/£)r„+i(W/£)r. dkx (A.8) 

with 

i f+L/2 f+L/2 , 

A nm (k x ) = -== \ \ / A(r) e -ik x r x -i(2nn/L)r y -i(2* m /L)r„ £ f (A. 9) 

V27rL 2 J-L/2 J 

We may ignore the degrees of freedom which do not correspond to propagation along the x-direction 
so long as the equations of motion do not couple these to those which do correspond to propagation 
in the x-direction. We have noted above that in the absence of particles the field is not coupled 
between different points in Fourier space, so that whether or not we may ignore the field variables 
that are not on the line n = m = will depend upon the coupling of the field to the particles. 

Let us first consider the part of the Lagrangian which remains in the absence of particles: the 
Lagrangian for the free field. Written in terms of the Fourier variables this is now 



L = e Idk 



A nm (k x ) ■ A nm {k x ) c 2 k x A nm (k x ) ■ A nm (k x ) (A.IO) 



As the Lagrangian determines the equations of motion for each of the Fourier variables independently, 
we may ignore the equations of motion for those for which norm are not equal to zero by simply 
dropping them from the Lagrangian. The one-dimensional Lagrangian that results is 

L = e Jdk x [A*(k x )-A(k x )~c 2 k 2 x A*{k x )-A(k x )~] , (A.ll) 

where A(k x ) = Aoo(k x ). In fact the integral need only extend over half the real line because 
the field for negative k x is determined entirely by the field for positive k x . In particular we have 
A(k x ) = A*(—k x ). This is due to the fact that the spatial vector potential is purely real. Also, in 
general there are two polarisations, so A(k x ) = (0, A y {k x ), A z (k x )). However, we will consider only 



one polarisation, the y-component, writing this simply as A{k), where k = k x . From Eq.(A.ll) the 
Lagrangian determining the equation of motion for A{k) is 

L = £ o y"dfc[|i(fc)| 2 -c 2 fc 2 |A(fc)| 2 ] . (A.12) 

Writing x for r x , a notation we will keep from now on, the part of the real-space field that interests 
us is given by A(x)e y , where e y is a unit vector in the y direction and 

A(x) = -jL= J A(k) e ik * dk. (A.13) 
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Cavity 



x=0 



x -> 



Figure A.l: A diagram of the small system (in this case a cavity) which interacts with an external 
field extending from x = out to infinity. Mirror 1 is totally reflecting so that neither the cavity 
modes or the external field to the right of the cavity interact with the external field to the left of 
the cavity. Mirror 2 is partially reflecting, so that there may be an exchange of energy between the 
cavity and the external field extending to the right. 

Dropping Half the Real Line 

Before adding an interaction Lagrangian we make one further alteration, and that is to treat the 
real field on only half of the real space only. This makes sense in the case of optical cavities with 
only one output mirror. Referring to Fig jA.l| , since mirror 1 is chosen to be completely reflecting, 
the field to the left of the cavity is coupled neither to the field to the right of the cavity, nor to the 
field inside the cavity. 

We return briefly to mention the standard treatment of the quantum light field. We have indi- 
cated above that taking the Fourier transform of the field simplifies greatly the equations of motion. 
We may refer to taking the Fourier transform as decomposing the field into travelling wave modes. 
The Fourier fields are not independent for all k however, and the standard approach is to take linear 
combinations of the Fourier fields E and B (or equivalently A and II) to obtain a new operator 
which is independent for all k. This new operator has the advantage that it decouples (diagonalises) 
the equations of motion between A and II, and so obeys a very simple equation of motion in the 
case of a free field. The operator in question is a(k), which satisfies 



and is referred to as the mode operator for the plane wave mode k. The reader is assumed to be 
familiar with this standard treatment of non-relativistic Quantum Electrodynamics. 

Now that we wish to consider only half the real line the Fourier modes we have been using above 
are not the most convenient. Since a cosine expansion is sufficient to represent the even extension 
of A(x,t) on the real line, we may use instead the transform 



[a(k),at(k')]=5(k-k'), 



(A.14) 




(A.15) 



which we will refer to as the cosine transform. The transformed field is real because A(x, t) is real, 
and is only defined for non-negative k. Using the standard Fourier transform relationship, the inverse 
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transform may be calculated as follows: defining A(k) and A(k) by 

**•<>-{ £Zi til ■^-{ AlkJ l til 

we may write 

poo 

A ^ = v^J dkA{k > t) ( elkx + e ~ lkx) (A - 17) 

/OO rOQ 
dkA(k,t)e lkx + -jJ=^ dk A(-k,t)e lkx (A.18) 
-oo J — OO 

OO 



= 7^ I MA(k,ty k *. (A.19) 

This is just the standard Fourier transform. Defining A(x 7 t) to be the even extension of A(x J t) 1 we 
may use the inverse Fourier transform relation to obtain 

/oo poo 

dx A(x, t)e~ lkx = / dxA(x, t) cos(kx), (A.20) 

-oo J 

giving the inverse cosine transform as 

A(k,t) = J^ dxA(x,t)co$(kx). (A.21) 
v w Jo 

We shall denote the cosine transform of the canonical momentum field U(x) by 7r(fc). The Lagrangian 



in terms of this cosine variable is the same as given in Eq.(A.12), but the integration over k is now 
only over non-negative k. 

In choosing to represent the field by a cosine expansion, rather than, for example, a sine expan- 
sion, means that we have chosen the boundary condition at x = such that an incoming wave is 
reflected without a phase change. This can be seen by noting that in this case the field for positive 
x is effectively half of a symmetric field extending over all x. Therefore, a wave travelling towards 
x = from positive x has a twin wave travelling in from negative x. As the first wave passes the 
origin the twin wave also passes the origin, and appears to an observer in the positive x region as a 
reflection of the original wave with zero phase change. If we had used a sine expansion instead this 
would have resulted in a 7r phase change upon reflection. However, we may use any phase change 
we wish, as the correct phase change in any experimental realisation can be correctly accounted 
for simply by multiplying the output field by the correct phase factor. Because of this, different 
conventions are used by different authors. Here we follow Gardiner |48) in setting the phase change 
to zero for the purposes of the derivation. When we consider the input-output relation for an optical 
cavity in the final section we will return to this point to indicate the various conventions. 



The System-Field Interaction 

We now consider the interaction Lagrangian. We have simplified the initial free Lagrangian, in 
a manner that makes sense for the system we are treating. This is somewhat phenomenological, 
since we are assuming that devices such as perfect reflecting mirrors exist, without treating their 
interaction with the light field in a detailed manner. In choosing an interaction Lagrangian we 
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make a further phenomenological departure from the standard electro-magnetic Lagrangian given 



in Eq.(A.l). We drop the matter- field coupling given in that equation and introduce an interaction 
which will again make sense for optical cavities. We assume that there is a small system (perhaps 
a field mode inside an optical cavity) which is coupled solely to the degrees of freedom of the field 
that are propagating along our chosen direction by the interaction Lagrangian 



O a dk K(k)A(k), 
Jo 



(A.22) 



where O a is some operator of the small system. This interaction was presumably motivated by 
the coupling of an atom with the electro-magnetic field in the electric dipole approximation, which 
leads to spontaneous emission. The important feature of this interaction is that it couples a system 
operator to many modes of the electric field so that this field becomes a reservoir for the system. 



When we choose a particular operator for O s to treat the case of an optical cavity in Section [A.4| , 
it will become clear our choice models the exchange of photons between the cavity and the external 
field. 

The system operator O s is coupled to the electric field over a frequency range given the coupling 
constant n{k). The frequency range (band width) of the modes to which the system is coupled is 
inversely proportional to the length over which the system is coupled to the light field in real space. 
We will find it expedient eventually to take the ideal case in which n{k) is flat, as this will simplify 
the treatment. 

The Lagrangian for the system plus reservoir may now be written 



L — i S ys(Z 



£» 
2 



dk 



A(k) 2 + c 2 k 2 A(k) 2 + O s / dfc K (k)A(k) 



(A.23) 



where L sys (Z) is the Lagrangian for the system alone. 



A. 2. 3 The Hamiltonian and Quantisation 

To quantise the field we calculate the Hamiltonian which corresponds to the Lagrangian and define 
the field canonical momentum by 



7r(fc) = 



dL 



dA{k) 

The resulting Hamiltonian is given by 
H = H sys (Z) + ^^ 



= e A(k) + O s K(k). 



-(7T{k) ~ O s n(k)) 2 + c 2 k 2 A 2 (k) 



dk, 



(A.24) 



(A.25) 



Note that if we were using the standard Lagrangian (Eq.(A.l)) the canonical momentum would be 
simply the time derivative of the vector potential, where as due to the interaction we have chosen, the 
canonical momentum contains the extra term O s K,(k). Note also that for our treatment to be valid, 
we require the interaction we have chosen to be a good approximation to the true Lagrangian, ie. 
for the system we are treating, the true, and exceedingly complex, Lagrangian should approximately 
reduce to the one we are using. If this is the case, then the conjugate momentum should also be 
a good approximation to the truth, ie. it should be approximately the correct thing to impose the 
commutation relations upon. 
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The equations of motion for A(k) may now be derived either by using Lagrange's equations, or 
by using Hamilton's equations, given by 

. dH . &H 
A = ^ ' * = (A ' 26) 

Quantisation of the field is accomplished by imposing commutation relations between the field and 
its canonically conjugate momentum such that the equations of motion may be written 

A=^[A,H] , ir=^[7r,H}. (A.27) 

That is, the commutators must be chosen to satisfy the relations 

[A,F(A,n)]=in-F(A,Tr) , [tt, F(A, tt)] = -ih— F(A, tt), (A.28) 

where F(A, n) is an arbitrary function of A and ir. The commutator which satisfies the above 
criteria is 

[A(k, t), w(k', t)\ = ihS(k - k'). (A.29) 

Equations of Motion 

There are now three ways to calculate the equations of motion: using Lagrange's equations, using 
Hamilton's equations, or calculating commutators. It is the third (calculating commutators) which 
is most often used in quantum mechanics. Performing this calculation, the equations of motion for 
the fields are 

A(k,t) = —7T(k,t)-—O s (t)K{k) (A.30) 

7T(M) = e uj 2 A(k,t) (A.31) 

In order to obtain mode operators for the mode functions cos(kx) we need to take a linear com- 
bination of the fields A(k,t) and n(k,t) that both decouples this set of differential equations, and 
satisfies delta commutation relations. The correct linear combination is 

H^t) = J^A(k,t)+i^=L=ir{k,t) , u^k/c, (A.32) 

where the dependence on e Q uj is required to decouple the differential equations, and the scaling by 
\/2h is required to obtain the correct commutation relations. The differential equation governing 
the mode operators is 

b(k, t) = -iujb(k, t) - j-^—O s {t)n{k). (A.33) 

Since this is a first order differential equation, given a choice for b(k) at a particular time to, the 
solution is determined for all t both prior to, and later than to. In particular, the solution for the 
above equation is 

KM) - KMo)e-" (t - fo) - \f^< k )ff lu{t ~ T) °s( T ) dr - ( A - 34 ) 
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In the absence of the system (ie. O s = 0) the field mode operators evolve in the simple manner 

b(k,t) = b(k,t )e-^ {t - to \ (A.35) 

which gives the evolution of the free field. Writing the vector potential in terms of the mode operators 
we have 



A(x,t) 



o V ttsquL 2 
I ft/4 



o V ttEoUjL 2 



[b(k, t) cos(kx) + b^ (fc, t) cos(fca;)] dfc 



[b(k,t)(e ikx + e~ lkx )+h.c.] dfc. 



(A.36) 
(A.37) 



In the absence of the system the solution for the vector potential may be written in terms of the 
mode operators at some arbitrary time to as 



A(x,t) 



I ft/4 
o U ne ujL 2 



b(k, t )( e *( te -"( f -*o)) + e -i(k*+oi(t-t ))) + h. 



dfc. 



(A.38) 



Since e l ( fcx±a, (* is a wave travelling in the positive (negative) x direction, we see that each mode 
operator b(k) is associated with two waves of frequency to travelling in opposite directions. 



A. 3 Input and Output Fields 



Let us first consider the free field which is given by Eq.( A.3£ ). At x — this field is 

b{kM)e~ iu{t ~ to) +h.c. dfc. 



^free(0,t) - , y 

J V TTS a UjL 2 

By a simple rearrangement of A(x,t) it is seen that 

A licc (x, t) = A(0, t - f) + A(0, t + f). 



(A.39) 



(A.40) 



The total field can therefore be written as the sum of an inward travelling field and an outward 
travelling field. These inward and outward propagating fields may be written, as we have done, in 
terms of the field at x = 0. That is, if we know the field at x = for all time, then we may construct 
the entire field for all time. This may be understood as follows. If we know the field at x = 0, then 
this is the result of a field which has propagated along the x-line from positive infinity to x = 0. 
Also, now that the field has this value at x = 0, this will propagate in the positive x direction back 
out to infinity. The field at a particular position x and time t is therefore the sum of a field which, 
at a time x/c later will reach x = 0, and is therefore A(0,t + f ), and another field which will have 
come from x = and taken a time x/c to reach position x, and so is therefore A(0, t — § ). 

Let us denote the free field at x — by Ao(t). With no system interaction, this operator would 
be equal to the actual field, A(0,t), for all time. However, once we include the system interaction 
the evolution of the field is no longer given by multiplying the mode operators by an exponential, 
but is instead given by Eq.( |A.34D . This means that A (t) is only equal to the field at x — when 
t = to, and will in general be different from it for all other times. 
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Writing A(x, t) in terms of the mode operators at t = to, but using the full evolution of the field 
in the presence of the system, it is possible to write, after some careful manipulation |48[ 

x+c(t— to) 

A{x,t) = A Q {t+i)+A {t -*)-!/ {2c) j K {x')Olt-\^\) dx', (A.41) 

x — c(t— to) 

where K(x) is the inverse cosine transform of n(k). We will now examine closely the last term in 
this expression, which is the contribution of the system to the field. Recall first of all that we take 
the interaction of the system with the field to be non-zero only over some small range close to x = 0. 
In particular let us choose k(x) to be zero outside the interval [0,Ax]. Now let us examine the 
third term for values of x outside the interaction region. In this case x' is always less than x in the 
integral, so that we may write the system contribution as 

x+c(t—t ) 

A sys (x, t) = -l/(2c) J K{x')O s (t - a=2l) dx 1 . (A.42) 

x — c{t — to) 

Taking t to be far enough in the past so that the interval [x— c(t — t ), x + c(t — 1 )] covers completely 
the interval [0, Ax], we may simplify the above expression further to obtain 

A sys (x, t) = ~ [ K(x')O s (t - £==i) dx'. (A.43) 
Zc Jo 

The system contribution to the field at position x and time t is clearly an integral over the range 
over which the system operator O s interacts with the system. The contribution from each point 
of the interaction being the value that the system operator had at a time which is earlier than the 
current time t by the amount required for the field to propagate from that point to x. Outside the 
interaction region the field generated by the system is propagating in the positive x direction, which 



is clear since we may write A sys as a function of t — § . Returning to Eq.(A.42), it is clear that the 
lower bound on the integration, namely x — c(t — to), is there because the contribution from points 
which are further away from x than c(t — to) should not be included, as the field will not have had 
time to propagate to x in the time interval t — to. 



What Eq.( A.41 ) tells us is that if we choose the state of the field at t = to, then for times later than 
to, the field is what it would have been from free evolution (given by calculating Aq — A(0, t) using 
free evolution and writing A(x,t) = Ao(t + f ) + Ao(t — f )), plus a term generated by the system, 
which is propagating outwards. This is a fairly intuitive result, given that Maxwell's equations 
generate propagating solutions. 

What we have said so far is, in a sense, all we need say about the input and output fields. The 
in-field may be specified independently of the system, by giving an initial condition for the full 
external field. The out-field is then given by the in-field plus a system contribution. 

Note that as the in- and out-fields are propagating, we may sit at any distance away from the 
system to observe them. Changing the position of observation merely effectively shifts the origin 
of the time axis for the observations. As this is the case, we may evaluate the input and output 
fields at any given position and call this the time varying input to and output from the system. For 
convenience we will choose the position to be x = 0. Naturally any real observation would be at 
some finite distance away from the interaction region. 
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To obtain the central result of this section we need merely modify Eq.(A.41) by taking the limit 
as t — > — oo and setting k(x) — kS(x). After identifying A (t) as the input to the system, (the 
inward travelling field evaluated at x = 0) and the other two terms (also evaluated at x = Q) as the 
output from the system, this gives us 



K 



A out {t) = A iD (t) - —O s {t). 

2c 



(A.44) 



This allows us to calculate the properties of the output field in terms of the input field and the 
system evolution. 

At this point we have essentially completed our task in that we have derived the input-output 
relation. However, in order to gain an understanding of the way in which the standard treatment of 
Input-Output theory links together the output field with a final condition for the external field (as 
opposed to an initial condition) we will finish with a discussion of this point. 

To do this we will now examine the field prior to to. Of course the solution we have written 
down above is valid for t earlier that to as well as t later than to, but we will shortly write it in a 
more intuitive form for the former case. Note first that the field at t = to contains nothing of the 
field emitted by the system at times later than to, which is reflected in the fact that the solution for 
later times contains the outward propagating system contribution in addition to the freely evolved 
terms. However, if we wish to look at the solution for the field at an earlier time t, we should note 
that the field at t does contain the contribution radiated by the system for the time interval t — t. 

Recall that when we use only free evolution, the out-going field is simply a reflection of the 
in-going field. That is, the in-going and out-going fields are the same. Since, as we have mentioned 
above, Ao(t) contains the field emitted by the system (up until to), the in-field that we calculate 
using free evolution also contains this field up until to, but travelling inwards. In the full solution, 
however, the in-field will obviously not have this extra inward travelling system contribution. The 
full solution will therefore be the sum of the terms due to free evolution and an inward travelling 
term which exactly cancels the extra inward travelling term contained in the free solution. Indeed, 



with a little re-arrangement the full solution, given by Eq.(A.41), may be written 



c — c(t— to) 



A(x,t) = A (t + ±) + A (t - 



,) + l/(2c) J K{x')O s {t+\^\) Ax'. (A.45) 



x+c(t — t Q ) 



Note that the last term of this equation is inward propagating. Eq.(A.41) is the intuitive expression 



to use for the field at times later than to, and Eq.( A.45 ) is the intuitive expression to use for the 
field at times earlier than to- 

To sum up, if we specify a state for the external field at some time (ti, say) after the time at 
which we wish to consider the system, then we should note that the out-field given by free evolution 
contains the system contribution up until time t\ (and is therefore the correct output field up until 
that time). Alternatively, if we specify the state of the field at a time (to, say) earlier than the time 
at which we wish to consider the system, we should note that the in-field is the correct input field 
for all times later than to- This prompts the following definitions: 



A in (t) = 
A out (t) = 



H£ UjL 2 



TTS UjL 2 



b{k,t )e~ luj{t ~ to) 
b(k, tije 



-iuj{t—t\ ) 



h.C. 



h.C. 



dfc = Ao(t), 



Ah, 



(A.46) 
(A.47) 
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where to is a time in the past, and t\ is a time in the future. Note that A out is simply A lri with 
to replaced with t-y. From the above discussion we can conclude that A out (t — x/c) is the correct 
output field (up until time ty), and A ln (t + x/c) is the correct input field (for times after to). 

By equating Eqs ( A.41 ) and ( A.45| ) (with A (t) replaced by A[ n (t) in Eq.( A.41 ), and replaced by 
A ou t(t) in Eq.(A.45)) we find, as we expect, that for x outside the interaction region 



1 f°° 

A ont (t + x/c) = A in (t + x/c)- — / dK(x')O s (t+\^L\) dx', 



(A, 



where we have taken the limit to — > — oo and t\ — > oo. Since this is true for any t, even though x 
must be outside the interaction region, we can see that 



A out (t) = A m (t) - 1 / <l/v(.r'!Oj/ - -|i ,1,-'. 



(A.49) 



Now Ai n (t) is the input field at x = 0, and A out (t) is the output field at x = 0. Setting k(x) = kS(x) 
we have 

h ' ■' (A.50) 



Aout(t) =A ta (*)- — O s (t) 
2c 



In addition, using the previous result and Eq.(A.4I) (with A (t) replaced by A in (t)), we find, again 
as we expect, that 

A(x, t) = A in (t + x/c) + A out {t - x/c). (A.51) 

That is, the total field is the sum of the input field propagating in, and the output field (which 
includes the system contribution) propagating out. 



A. 4 Input-Output Relations for a Damped Cavity 

In the previous section we examined the input-output relations for a general open system. To obtain 
the input-output relation for a lossy optical cavity in their final and most useful form, it is necessary 
to make the rotating wave approximation. This requires that the frequency of the cavity mode is 
much larger than the rate at which the system changes due to any other dynamics, such as the 
interaction with the external field. 

To obtain the correct equations for a damped cavity, we choose the coupling operator to be a 
quadrature of the cavity mode, which gives dynamics linear in the mode and the field operators. 
Any quadrature operator will do, as rotating this quadrature simply changes the phase relationship 
between the input field and the cavity mode. We choose O s = p = —i(a — a^) as this gives the 
cavity mode the same phase as the input field. This would be true, for example, if the input 
mirror was a simple dielectric boundary. At this point we make the rotating wave approximation, 
which involves moving into the interaction picture and dropping terms from the Hamiltonian (and 
therefore from the equations of motion) which are oscillating (rotating) very rapidly compared to 
the remaining dynamics (due to the interaction) . These may be dropped because their contribution 
averages to zero in the time-scale at which the dynamical variables are changing. Thus the terms 
ab(k) and a^b^k), which oscillate in the interaction picture at frequencies greater than u>o, are 
removed from the interaction Hamiltonian. As the remaining terms, ab^{k) and a'b(k) will only 
contribute significantly to the dynamics when the difference between lu = ck and u>o is much smaller 
than ll>o, again due to this averaging, we perform the further approximation of setting u = loo 
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in the integral over the field mode operators, so that it may be removed from the integral. This 
may be though of as a narrow bandwidth approximation, and assumes that we only need consider 
frequencies in a bandwidth about Uq which is much smaller than loq. Before obtaining the final 
form of the Hamiltonian we make one further approximation which is to extend the integral over the 
mode wave-number fc to negative infinity. While this is physically absurd, the extra terms are once 
again rapidly oscillating, and contribute negligibly to the dynamics. The resulting Hamiltonian is 



H = hcjoa^a - iJ -^^K [ (& t (fc)a - tfb(k)) dfc + - 



^ 2 {k)+c 2 k 2 A 2 {k) 



dfc. (A.52) 



In fact there is another term which comes from the interaction Hamiltonian, but it consists only 
of system operators, and results merely in a shift of the frequency of the cavity mode. This may 
always be taken into account by redefining uiq. The equations of motion for the system and reservoir 
operators are now derived easily by taking commutators with this Hamiltonian. We have 



/>oo 

-i^a + * I - / b(k) dfc (A.53) 



2e hL 2 ~ 



(A.54) 



where we have written the solution for the reservoir operators in terms of some initial conditions and 
the system operator. Substituting this solution in the equation of motion for the system operator, 
we have 

a = -iuj a - -a + ^/jb in (t) (A.55) 

where 

Mt) = / b ( k > to)e-^ {t - to) dfc , 7 = K\ (A.56) 

V27TJ-00 e hL z 

which is the negative frequency part of the input vector potential. In particular, continuing to make 
the narrow bandwidth approximation, we have 



A ta (t) = \j 2^Z*(M*) + &L(*))- (A.57) 

It is clear now that the equation of motion for the cavity mode annihilation operator has two terms 
due to the interaction with the external field. The first of these is a damping term, and the second 
is a driving term, being the negative frequency part of the input field. 

Defining b ou t(t) in a similar manner to bi n , that is, as the negative frequency part of the output 
vector potential, and using the same scaling, we find that the input-output relations become 

&out(<) = hn(t) - ^ya(t). (A.58) 
Note that to obtain this relation we have used the 'convention' that there is no phase change upon 



reflection at x = (see subsection A.2.2 ). To change this convention, all we need do is to multiply 
the output field by a phase factor. In particular, if we want a 7r phase change on reflection, we 
multiply the output field by minus one, which is why one will often see the relations written 

6out(<) = V7<»(*) - M*)- ( A - 59 ) 
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It is this latter convention that we use for our analysis in Chapter ||[ 

We now wish finally to make a connection with photo-detection theory to see the physical signif- 
icance of bi n (t) and b out (t). Glauber's theory of photo-detection gives the rate of photon detections 
for a detector placed in the field E(x,t) as proportional to (E + (x,t)E~(x,t)), where E + (x,t) is the 
positive frequency part of the Electric field operator p6| . We are interested here in detecting the 
light output from the cavity, so the photon detection rate is given by (E^ ut (x,t)E~ ut (x,t)} . Using 
E ut(x,t) = -dtA out (x,t) we have 



The rate of photo-detection is therefore proportional to (& ut (t)b 1 out (t)} , and we just need to determine 
the constant of proportionality. Since 7 is the rate at which photons leave the cavity, for an detector of 
unit efficiency, and for a vacuum input field, the average rate of photo-detection must be 7(6* (t)b(t)), 
being the average number of photons in the cavity field at time t. Calculating (bl, ut (t)b ovl t(t)) gives 
exactly that result, and hence the constant of proportionality is unity. We see now that b out (t) is 
simply the positive frequency part of the electric field scaled such that the average photo-counting 
rate is (b ont (t)b QU t(t)) . The operators b out (t)b ou t(t) and b in (t)b ln (t) may therefore be interpreted as 
photon flux operators, so that (b hl (t)bi n (t)) is the rate at which photons are incident on the cavity 
from the input field. This concludes our discussion of the input-output relations. 



Appendix B 

The Action of Exponentials Linear 
and Quadratic in P and Q 



We first calculate the effect of an operator of the form 

tfP+nQ ( B1 ) 

on a coherent state \a). The coherent state is defined as the eigenstate of the annihilation operator 
a, such that 

a\a)=a\a), (B.2) 

and 



/ muJ /I 

a = y^n x + l2n^ p - (R3) 

Here m and ui are the mass and frequency of a harmonic oscillator which serves for the purposes 
of defining the coherent state. In particular we are interested in the position wave-function of the 
result. We therefore wish to calculate 

(a#) = {x\e vP+ " Q \a), (B.4) 

where \x) is an eigenstate of the position operator Q such that Q\x) = x\x). Note that in general 
will not be normalised. To perform this calculation we will need the BCH formula given in 



Eq.(3.41), and the position wavefunction for a coherent state, 



„2\l/4 /o„2\l/4 



2' S I e -s 2 x 2 +2sxa-±(\a\ 2 +a 2 ) _ f ^ S \ £ -s 2 x 2 +2sxa-a 2 -ia r ai ( B *) ) 



where 



3 = (R6) 

a — a r + 1014. 
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Note that this expression contains the phase factor — ia r ai. This is left out in many texts, but is 
essential for consistency with the completeness relations for the position states. We also require the 
inner product of two coherent states, 

(a\(3) = e -i(M a +l/»l a )+°V (B.8) 

and the well known integral formula 

J e - a x 2 ~f3x dx = y/Z e /3 2 /(4 a ) ; Re [ a ] > Q ( B 9 ) 

We proceed first by rewriting the exponential in terms of annihilation and creation operators, so 
that we have 

(B.10) 



p vP+nQ _ e ea+<t,a'' _ e <j,a^ e 9a e 8<f>/2 



in which 



W (B.ll) 



2muj V 2 

We may now use the completeness relation for the coherent states to obtain 



(x|e^V»e^/ 2 = i ff(x\(3)(P\e^e ea \a)e e ^ 2 d 2 (3 



= (x\a + <f>)ei^ 2+Rc ^ +ea+e ^ 2 . (B.12) 

We see that the state remains coherent, although it is no longer normalised, and is shifted in phase 
space by <f>. 

We now wish to calculate the effect of an operator of the form 

e vP 2 +CQ 2 +iQP (B.13) 

on a coherent state. This time we require to calculate 

(a#) = (x\e^ p2+<:Q2+ ^ QP \a). (B.14) 
For this calculation we will need the disentangling theorem for the exponential of a general quadratic 



form of the annihilation and creation operators, which is given by [107] 

e ua 2 +va u +wa ] a _ e (w+x) /2 e la f2 e xa f a e ma 2 -^g-j 

in which 



1 = J^Hf)^ (BJ6) 

X = M l ^ tt J r^rjr) (B.17) 



/coth(/) - uisinh(/) / 

™ - /cot,,"/)-. C- 18 ' 
/ = x/(w 2 - Auv). (B.19) 
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First of all rewriting the exponential containing P and Q as an exponential in the annihilation and 
creation operators, we have 



(x|e" p2+ ^ 2+ ^ p |a) = (x\e ua2+va ' 2+wa ' a+u \a) 



(B.20) 



in which 



u = 



w — 



( Cft ymhw ^ 

2mw 2 2 ' 



+ rjmhui 



(B.21) 
(B.22) 



We now proceed by using the disentangling theorem, and employing the completeness relation for 
the coherent states. 



(a#) = (x\e^ w+ ^/ 2 e la ' 2 e^ a e ma2 \a} 



(">+x)/2„ia t p xa'a ma 



a) (f(3 



= i jf{x\P){fi\t 

= l e kH 2 (\e 2 -\-l)+rn^ e v + ( W +k)/2 ff^^X) ^ ^ (R23) 



Performing the integral over the real and imaginary parts of a, we obtain 



(x\ 



2s 



VTT27 V 77 



2\ 4 



e i|Q| 2 -ma 2 e t/+(«,+fc)/2 



x exp ^— s a; 
x exp ^2sxae x 
x exp |o; 2 e 2x 



1-21 
3-21 
1 



1 + 2 



1 - 21 



3-21 
1 



1 + 2- 



1 + 2/ 
1 - 21 



3-21 



1 + 2/ 
2 



1 

2 + 1 + 21 



(B.24) 



It is easily verified that this reduces to (x\a) as required when we set I = \ = m = 0. 
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Appendix C 

Detection Probabilities Following 
Two Consecutive Atom-Cavity 
Interactions 



We consider a composite system, consisting of a cavity and two atoms, the initial state for which is 
given by Eq. (5.32). We allow the first atom to interact with the cavity for a time t, and at some 
time after this allow the second atom to interact with the cavity for a time s. Using the interaction 
picture throughout the calculation, so that the Hamiltonian for the interaction of each atom with 
the field is given by the term proportional to k in Eq. ( 5.17 ), the state of the system after both 
interactions is given by 



|*)= r e ie °\0,g,g) + (r*r e"°+* - iT^e* 1 ) |0, e, g) 
+ (T s r e ie ° - iTgrfoe* 1 - T^T*r oe ^) |0,.g,e) 
+ (T^ine^ - iT» o rlr e ie ° + * - iT s r e ie °) \l,g,g) 

oo 

+ J2 [r s (rV n e ie "+* - l TV„ +1 e ie "+ 1 ) - iT s (r + V n+1 e i < 9 »+ 1+, » - *T + t r„ +2 e ie "+ 2 )] \n,e,e) 



(C.l) 



n=0 
oo 



]T [r s (r + V n+1 e^ e »+ 1+ « - iTlr n+2 e w ^) - *T s (rV„e ie "+^ - i T t r rl+1 e je "+ 1 )] \n+l,e,g) 

71=0 
OO 

[r s (Tlr n e t0 " - *T + t r„_ 1 e^"- 1 +«) - iT s (T' r n+1 e ie ^ - t T t r n e^+^)] \n,g,e) 



71=1 
OO 



In this expression, the labels on the kets refer respectively to the cavity photon number, the state 
of the first atom, and the state of the second atom. The joint probabilities for detecting n photons 
in the cavity, and the atoms in either their ground or excited states, are given by the following 
expressions: 



Ql 3 



Pa 



<)<) 
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Qf = P (TS + T t fg + 2TCT*cos(</»)) + P 1 fgf* 
-2^I\r s T s T t sm(Ae 1 ) 
+2/fV\fgr*T* sin(A0 1 - <P) 

Qn>2 = Pnf_*rf+P n _ 2 fff S 

+P„_ 1 (f_ s T_* + fj + 2r s T s r = t T t cos(<£)) 
-2 v /P n P„_ 1 f s r t T t sin(A#„ - <j>) 
-2 v /P„_iP„_ 2 r = t T t sin(A0„_i - 0) 

-2VPnP„-ir s T s r t ri sin(A# ri ) 

-2VPn-iP^ 2 r s T s T t T t sin(A0„) 
-2 v /P n P„_ 2 r s T_ s r t T t cos(A0„ - 0) 

Of = Po(fg + fgf*-2TCT*cos(</»))+P 1 f*fg 
+2/fVPrgT* sin(Afli) 
-2v^Plr*T^ sin(Afli - 0) 
0r>! = P„_if'T* + P n+1 f'f* 

+P„(f s f t + f s f ' - 2r s T s r_*T t cos(^)) 
-2 v /P„P l -if s r t T + t sin(A#„ - <t>) 
-2 v /PnP„+iT s r t T t sin(A0„ +1 - 0) 

+2VPnPn+ir s T s r t r t sin(A0„ +1 ) 

+2VPnP„-ir s T s T + 'T* sin(A0„) 
+2 v /Pn + iP«-ir s T s T + *r t cos(A0„ +1 - 0) 

= p f* +p 1 f^ + v ^p7r*T^ 

q: 9 >! = p„_itff_ t +p n+1 f s f t 

+P„(fff« + Tff « - 2rfT_ s r t T t cos(0)) 
+2 v /P n P Il+1 f s r t T t sin(A0 n+1 - 0) 
+2 v /PnP„-iX s r t T* sin(A#„ - <P) 
-2 v /P„P„-ir s T s f t sin(A0„) 
-2 v /PnP l+ ir_ s TfT*T_ i sin(A0„ +1 ) 
+2 v /P Il _ 1 P„ +1 r s T_T_ t T t sin(A0„ +1 - 0) 

Qn>0 = P„f S f + P„ +2 f S f + * 

+P n+ i(P s T* + T s f + * + 2r s T s r + *T* cos(») 
+2 v /P„P„+if s r t T t sin(A0 n+1 - $ 
+2 v /Pn+iP n + 2 T s r + t T + t sin(A0„ +2 - 0) 

+2 v /P„P„+ir s T s r t r + t sin(A0„ +1 ) 



101 



+2 v APn+i J p« +2 r s r s r t r + t sm(A9 n+2 ) 

-2 v /P„ J P„+2r s T s r t T + t cos(A0 n+2 - <f>) 

In the above expressions Q^ e , for example, denotes the joint probability for detecting n photons in 
the cavity and the first and second atoms in their ground and excited states respectively. 
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